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Preface 


The past two decades have been characterized by the dynamic devel- 
opment of many technical and scientific domains, and the development 
of machine tools has been on this same trajectory. The complexity of 
machine tools has increased, as well as their intersection with other 
fields, primarily electronics. The new generation of machine tools is very 
different from the previous generation, from which a large amount of 
experience has been gained. The modern systems of machine tools have 
new rules of design, research, and fabrication. 

This book presents the new research in this large domain of ma- 
chine tools and tools. The book's main contribution is the establishment 
of technical diagnosis by the aid of the evolution of machine parameters, 
which can be correlated with the functionality of the machine tool as 
well as that with the tool. 

Two factors have contributed to these developments: one is the in- 
crease in the complexity of machine tools and their integration in flexible 
manufacturing systems, and the second is the success of the monitoring 
and diagnosis of other fields. By correlation of these aspects with the 
development of computer hardware, the book details the replacement 
of conventional systems with virtual instrumentation. This represents 
the step from an “art” in the field of machine tools to the science of 
the diagnosis of these systems. The specialized software systems create 



the possibility of scientific diagnosis and monitoring of the technological 
systems. 

Together with the traditional instruments used in the monitoring 
and diagnosis of machine tools (investigating such aspects as forces, vi- 
brations, and acoustic emission), the expert system is presented together 
with neural processing of the information. 

The first chapter reviews general notions about technical diagnos- 
tics of machine tools. The second chapter is dedicated to the state of 
the art of research concerning vibroacoustic diagnosis of machine tools. 
Chapter 3 presents an original systematization and application of a pack- 
age of vibroacoustic methods to establish the technical diagnostics of 
machine tools. Chapter 4 describes the experimental research in diag- 
nostic analysis of the mechanical system of the feed kinematics chain. 
Virtual instrument packages for diagnosis by vibroacoustic methods are 
described in Chapter 5. Specific elements from the feed kinematics chain 
are diagnosed by virtual instrument packages in Chapter 6. Chapter 7 
deals with a neural approach to the problem of establishing technical di- 
agnosis for machine tools. In Chapter 8, we conclude by presenting the 
original research on which this book is based. In order to prove that the 
method works, a few case studies are presented in the three appendixes. 

This book is a joint effort with my former professor, Constantin 
Ispas, a well-known specialist in machine tools in Romania, and one of 
his graduate students, now Dr. Dan Boboc. The topic of the book was 
the main theme of Dr. Boboc’s dissertation a few years ago. 

The authors would like to express their gratitude to two of their 
graduate students, Dr. Radu Pavel and Liviu Luca of the Precision 
Micro-Machining Center at The University of Toledo, Ohio, for assisting 
with the manuscript. I give special thanks to my wife, Jocelyn, for her pa- 
tience in polishing the English of this book, which was written by native 
Romanian authors. I also thank Rita Lazazzaro and Barbara Mathieu 
of Marcel Dekker, Inc., for their help and for keeping us on track. 


loan D. Marinescu 
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General Notions About the 
Introduction of Technical 
Diagnosis for Machine Tools 


1 . 1 INTRODUCTION: THE PLACE AND THE 
ROLE OF DIAGNOSIS IN MODERN 
TECHNICAL SYSTEMS 

The diagnosis of technical systems can be defined as a process of func- 
tional faults and their causes, on the basis of data obtained by control, 
supervising, or monitoring. The rudiments of establishing technical di- 
agnosis have been present for a long time regarding the estimation of the 
functioning state of a machine or equipment. Thus, it was natural that 
the functioning anomalies be noticed by the operator of the machine, 
such as: modifications of level or type of emitted noise, too much energy 
consumed, low output, the vibration level increasing, or local or global 
excessive heating. On the basis of these observations the operator could 
intervene to correct the deficiencies or to decide to stop the machine. 

Even by this subjective elementary diagnosis method, some impor- 
tant advantages could be identified: 

Increased sensorial accuracy of the operators 

Enhanced “database” of the operator’s memory, built upon accu- 
mulated experiences 




The possibility to control, analyze, diagnose, and intervene during 
the machine’s functioning, before reaching the faults 

Objective diagnosis, on the basis of some precise techniques and 
methods, has been propelled during the last decade by two important 
factors: the large scale introduction of automated production processes, 
and the development of data acquisition, measuring, and signal process- 
ing equipment. In these conditions, it becomes necessary to use some 
supervisory/control systems to ensure machine reliability. 

What is the role of such systems? In short it can be described as: 

To supervise the functioning state of the production system to 
prescribed performances 

To reduce and even avoid accidental breaks in functioning by super- 
vising some significant factors during the technological process 

To avoid the limit situations, which can generate faults and dam- 
ages 

To reduce, even eliminate, fabrication waste by supervising every 
step in the evolution of the product, parallel with the reduction 
of materials consumption 

To analyze functioning state tendencies, locally or globally, in order 
to plan interventions 

The supervising/diagnosing system for the technological process may 
be identified as an auxiliary system that has the possibility of being 
integrated in a reaction buckle. 

The signal, from a mechanical, pneumatic, acoustic, thermic, or 
such nature, is captured by adequate sensors and then converted to 
electrical signals, and sent to a postprocessing block (Figure 1.1). The 
white noise and the parasite signals are eliminated here; and the charac- 
teristic parameters of the supervised signal are extracted. Interpreting 
these signals on the basis of adequate software, the supervising/diagnosis 
block, which is the ensemble of electronic devices, offers data for fault 
recognition and their causes, and identification of where and when the 
fault arises. Depending upon the integrating degree of the supervising/ 
diagnosing system, the reaction buckle can be: 

Half-closed (warning/alarming reaction; stopping of the technolog- 
ical process) or 

Closed (continuing the data acquisition, concentrating on the el- 
ement predisposed to fault, with the beginning of some testing 
procedures; the modification of some parameters of the tech- 



-pre-processing - 



Figure 1.1 Preprocessing block. 


nological process to improve the functioning state; warning the 
experts and stopping the process to the limit) 

1.2 MONITORING, THE FIRST 
STEP IN ESTABLISHING THE 
TECHNICAL DIAGNOSTIC 


Monitoring is defined as the activity of gathering information about the 
functioning state in a given system, by means of adequate observation 





of instruments and measuring apparatus in order to supervise and inter- 
vene for correction purposes. The axiom of maintenance by monitoring 
is that the intervention maintenance/repairing is done only when the 
measurements show its necessity (Figure 1.2). 

The monitoring technical system can have two functions: 

The protection or preventive function, which imposes a warning or 
automatically stops the equipment if it senses the possibility of 
a fault 

The analysis and prediction function, which selects the main state 
modifications and surveys their evolution; the system expert will 
provide efficient solutions for remedy before any final damage 

Preventive monitoring is recommended for machines that are not 
integrally doubled or where unplanned interruption would lead to large 
production loss. In this case, the monitoring is directed to those com- 
ponents whose failure (in a shorter time than estimated) would have 
serious consequences for the entire assemblage. 

Predictive monitoring is mainly oriented to those machines and 
equipment with a continuous technological flux, and which might have 
accidental stops or revision/repairing periods implying considerable pro- 
duction loss; it is also oriented to complex production systems where 
the damage of a subsystem would block the entire flux. In this case, the 



Figure 1.2 Levels of necessary intervention. 



functioning state is intermittently or continuously supervised until the 
modification of the supervised parameters is able to signal imminent or 
predictable spoiling, and then take the imposed measures. 

The stages that have to be crossed to build the monitoring schema 
are as follows. 

1. Make schematic descriptions of the machine, equipment, or in- 
stallation, divided on technological subassemblies and groups, 
from the functional and constructive point of view. 

2. Make schematic descriptions of each component, marking the 
controlled elements and the monitored control points. An ade- 
quate monitoring for the fixed control points and for arranging 
and retaining the recordings should be used. The number of 
control points depends on the system complexity and on the 
spoilage consequences. This number will determine the com- 
plexity of the monitoring schema and also the number of qual- 
ified personnel required for analyzing and supervising. 

3. The conceived monitoring schema will be implemented into 
the system and will provide information concerning the sys- 
tem functioning state. Based upon this information, graph- 
ics can be drawn manually or automatically which will show 
the evolution of the studied parameters, and indicate the val- 
ues considered normal. Preliminary diagnostics can be estab- 
lished, and preventive repairing can be planned. Previous spe- 
cial events, faults, and other evaluation criteria, in short, “the 
machine history,” can be added to these data by adequate 
description. 


1.3 APPLICATION OF VIBRO ACOUSTIC 
DIAGNOSIS TO MACHINE TOOLS 

Modern machine tools are dynamic complex assemblies composed of me- 
chanical subsystems (gears, bearings, cam mechanisms, couplings, belt 
transmissions), hydraulic equipment (pumps and hydraulic engines, dis- 
tribution equipment, etc.), and electric equipment (motors, contactors, 
etc.). Because of this structural complexity, it is very difficult to establish 
technical diagnostics for machine tools. From the same point of view, the 
diagnosis of the cutting tools is much simpler and much of the present 
research is oriented in this domain. 



1.3.1 Sources of Machine Tool Vibration 
and Noise 

Not only the functioning of subsystems and subsystems of the machine 
tool but also the cutting process generates vibration and noise. A rapid 
way to identify the sources of this phenomenon, which might generate 
damage, is represented by the connection (gradual, separate, or in di- 
verse combinations) of the functional elements of the machine (electrical 
motors, transmissions, main kinematic chain, feed kinematic chains, lu- 
bricating system, etc.), by analyzing each time the vibrograms result in 
the measurement points. Figure 1.3 presents the internal and external 
vibration and noise sources that determine the vibroacoustic behavior 
of the machine tools [76]. 

1 .3 . 1 . 1 Free Vibrations 

The free vibrations appear in the absence of some perturbing forces, as 
components of some transitory processes of the machine tool. These pro- 
cesses and vibrations have a relatively short duration, as the vibrations 
imply. 

Free vibrations develop with the self- frequency of the elastic sys- 
tems where they were born, and this phenomenon is interesting from the 
point of view that the generating transitory systems can determine or 
influence the dynamic behavior of the machine tool. Considering this, 
some important transitory processes can be mentioned. 

The on/off chip contact of the tool leads to deformation variations 
of the elastic fixing system in a transitory process whose duration 
surpasses the time of a complete rotation of the main shaft; this 
sometimes happens during the entire passage. Analysis of this 
process has established an exponential law of variation of the y 
deformation of the elastic system: 

s = “- eip (l(n!«)) (L1) 

where a Q is the nominal thickness of the cutting layer, t is the 
time, t is the duration of a rotation or technological cycle, and 
Rse , Ka are the static characteristics of the elastic system and 
cutting process. It must be noted that the numerator of the expo- 
nent represents the time constant T of the technological process. 
Generally its duration is appreciated to be t = (3 ... 5) T. 
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Figure 1.3 Internal and external vibration and noise sources that determine 
the vibroacoustic behavior of the machine tools. 











Accelerating/breaking of the mobile elements appears when the 
still state or the uniform movement of some mobile subassem- 
blages is changed, such as tables, sledges, supports, and tra- 
verses. The developing law of these processes is also exponential, 
similar to that previously presented. The experimental results 
have established the value of the time constant of this process 
at 0.3 to 1 sec, so the transitory process could take maximum 
5 sec. 

Inverting the direction of movement is a more complicated process, 
consisting of three phases: breaking the movement, stopping, and 
restarting in a new direction. This process is accompanied by 
a series of complex dynamic phenomena: shocks, changing the 
direction of the friction forces, the modification of the pressure 
repartition on the guiding, abrupt temperature variations in the 
kinematic couple, and so on. Also taking into consideration the 
presence of the plays from the kinematic couplings, it can be 
shown that changing the direction of movement is a random 
process with unfavorable consequences for the dynamic of the 
machine tools. The greater the number of masses in movement, 
the greater are the consequences. 

The study of machine tool free vibrations is important because 
the dynamic parameters of the elastic systems that compose them are 
determined with the help of linear differentially homogeneous equations 
which describe the free vibrations of these systems. 

1 . 3 . 1 . 2 Auto vibrations 

The autovibrations are kept-up vibrations, caused by excitation factors 
generated from the vibration movement. These can be as follows. 

1. The interdependence between the magnitude of the cutting force 
and the relative displacement between the tool and semiproduct: due to 
an external cause, the cutting force varies, making a resultant mechanical 
work positive per variation cycle, and the energy resulting from this 
keeps up the vibration process. The rise of autovibrations depends upon 
the values and the ratio of the rigidities from the elastic system, and also 
upon the position of the main rigidity directions, while the frequency 
of the autovibrator process is very close to the self-frequency of the 
elastic system, without being equal to it. Figures 1.4a, b, c illustrate the 
influence of the cutting process parameters (respectively u, s, t) over the 
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Figure 1.4 Influence of cutting process parameters (respectively, v, s , t) over 
the autovibrations amplitude. 


autovibration amplitude. It can be observed that the cutting depth has 
a very powerful influence on the autovibrational process and is able to 
raise the vibration amplitude to 140 /im. The depth of cut (when the 
stability of the elastic system to vibrations is surpassed), is called critical 
depth, and it is adopted as a parameter for evaluation of the vibration 
stability of the elastic system of the machine tools. 



2. The dependence between the friction force and the sliding speed 
from the kinematic couples: the autovibrations of this process manifest as 
jerky motions that appear with the movement of the mobile assemblies 
with speeds under a specific critical value. This phenomena is known 
as stick-slip. These autovibrations, also called “of relaxation,” can be 
attenuated by an adequate selection of the couple of materials, rational 
lubrication, and discharging of the contact surfaces. 

3. The phase difference between the variation of the cutting force 
and move away: this autovibration type highlights the copying kinematic 
chains. Their occurrence is favored by the presence in the elastic system 
of the nonlinearity of playing type, hysteresis, and instabilities. The fre- 
quency of these vibrations is approximately equal to the self-frequency 
of the excited elastic system. 

4. The regenerative vibrations appear due to the existence of sur- 
face waves that were manufactured at the previous pass and which, at 
a new pass, generate dynamic forces conditioned by the intensity of the 
dependence between the variation of the chip depth and the cutting 
force. The effect of regenerating vibrations then has the proportional 
variation of the autovibration frequency function of the main shaft or 
tool rotation. 

A characteristic aspect of autovibrations consists of the possibility 
of raising the amplitudes in time that lead to the trepidation occurrence, 
which is a very dangerous phenomenon. 

Concerning noise and acoustical emission, the machine tools cover 
a very large domain: 

Infrasounds (vibrations with a frequency lower then 16 Hz) due to 
structural dislocations, phase transforming, and the like 

Sounds (the audible domain: 16-20,000 Hz) due to the function- 
ing of diverse types of mechanisms that compose the kinematic 
chains and are grouped in low sounds (16-360 Hz), medium 
sounds (360-1400 Hz), and high sounds (1400-20,000 Hz) 

Ultrasounds (vibrations having their frequency higher then 20 kHz) 
due to the increasing occurrence of crackings, microfrictions on 
the level of cracked surfaces, microcollisions with and between 
the microscopic particles, and so on 

Noise measurement allows us to estimate the “silence” of the ma- 
chine tools when running out of job, from the point of view of the tech- 
nical criteria of execution and operator protection. It permits us to es- 
tablish if the level of maximum noise produced by a machine tool is 



within the admitted limits imposed by the standards. This also allows 
us to locate the machine elements that produce high-level noise and to 
make a comparison between the noise levels produced by the same type 
of machines. 

The machine tool whose noise level is measured has to be com- 
pletely equipped (with all the covers and guards), adjusted for proper 
functioning, and run in and installed on a basis plate under the same 
conditions that it is installed for exploitation. The noise level is mea- 
sured when the machine is idle-running with the usual rotations and 
advances. 

The measurement points are disposed on a measurement line whose 
trace in the horizontal plane is situated at a height of 1.5 m from the 
basis plane, and at a distance of 1 m from the contour line of the ma- 
chine tool. The nonnoisy prominences are not taken into consideration. 
For each type of machine tool, the location of the measuring points 
is established in a concrete way. Their minimum number is four, and 
other supplementary points can be chosen as a function of the machine 
tool sizes. The placement of the measuring points for different types of 
machine tools is shown in Figures 1.5a to k, and the admissible val- 
ues of the noise levels for these machines are presented in Table 1.1 
[64, 149, 150]. 

Acoustical emission is a succession of high frequency elastic waves 
(>100 kHz), generated by freeing the internal energy stored in a struc- 
ture. It allows detection of fissures and/or ruptures. When the cutting 
process is running on the machine tools, the sources of acoustical emis- 
sion are numerous: continuous or discontinuous chip forming, processed 
material deformation, Assuring of the semiproduct or the tool, friction 
among the tool, semiproduct, chip breaker, fracture, and collision of the 
chip. In addition to all of the above, other sources of acoustical emis- 
sion come from the functioning of some mechanical subassemblies (gear 
wheels, bearings) and high-frequency electrical sources. 

1.3.2 Requirements for Machine Tool 
Diagnostic Systems 

Utilization of the diagnostic systems should correspond with the impor- 
tance of the supervised system of the production system, with its com- 
plexity, and with its performance. Usually the monitoring is oriented 
to numerical command machine tools, processing centers, cells, flexible 
processing lines, and also to cutting tools. 




















Figure 1.5 Placement of the measuring points on different types of machine tools: (a) univer- 
sal grinding machine; (b) surface grinding machine; (c) centerless grinding machine; (d) shaping 
machine; (e) portal milling machine; (f) reaming and milling machine; (g) vertical milling ma- 
chine; (h) radial boring machine; (i) boring machine with column; (k) normal and revolver 
lathes. 




Table 1.1 Admissible Values of the Noise Level for Different Types of Machine Tools 


Machine 

type 

Turn 

rpm 

power 

kW 

0- 

-500 

500- 

-1000 

1000- 

-2000 

2000- 

-4000 

4000- 

8000 

8000-16000 

L 

dB 

(A) 

Cz 

L 

dB 

(A) 

Cz 

L 

dB 

(A) 

Cz 

L 

dB 

(A) 

Cz 

L 

dB 

(A) 

Cz 

L 

dB 

(A) Cz 

Grinding 

1.0-1. 6 









80 

75 

83 80 

machines 

2. 5-4.0 



75 

70 

78 

75 

81 

75 





4. 0-6. 3 



77 

70 

80 

75 

83 

80 





6.3-10.0 



77 

70 

80 

75 

83 

80 




Shaping 

1.6-2. 5 

75 

70 










machines 

2. 5-4.0 

80 

75 











4. 0-6. 3 

82 

80 











6.3-10.0 

85 

80 










Portal 

10.0-16.0 

79 

75 

82 

80 

85 

80 






milling 

16.0-25.0 

81 

80 

84 

80 

87 

80 






machines 

25.0-40.0 

83 

80 

86 

80 

89 

85 










Reaming and 

6.3-10.0 

milling 

10.0-16.0 

machines 

16.0-25.0 

Vertical and 

1.6-2. 5 

universal 

2.5-4. 1 

milling 

4. 0-6. 3 

machines 

6.3-10.0 

10.0-16.0 

Boring 

1.0-1. 6 

machines 

1.6-2. 5 
2. 5-4.0 
4. 0-6. 3 

Normal and 

1.6-2. 5 

revolver 

2. 5-4.0 

lathes 

4. 0-6. 3 
6.3-10.0 
10.0-16.0 



79 

75 


82 

80 


84 

80 


75 

70 


77 

70 


79 

75 


82 

80 


75 

70 


77 

70 





82 

80 

85 

80 





85 

80 

88 

85 





87 

80 

90 

85 





78 

75 

81 

80 

75 




80 

75 

83 

80 





82 

80 

85 

80 





85 

87 

80 

80 

87 

80 





75 

70 

78 

75 

80 

75 



76 

70 

79 

75 

80 

80 



78 

75 

81 

75 





80 

75 

83 

80 





77 

70 

70 

75 

85 

80 



79 

75 

82 

80 

85 

80 



81 

75 

84 

80 

90 

80 



83 

80 

86 

80 

90 

85 



85 

80 

88 

85 







In order to efficiently detect and supervise the spoilage that can 
appear, the supervising and diagnosis systems must follow a minimum 
set of requirements: 

Adequate match to the supervised machine tool or tool, having the 
necessary precision and ability to detect the faults. 

Fault warning must be made on time, and false alarms must be 
avoided. 

Location of the faults must be determined in order to minimize the 
intervention time for repairs. 

Correlation of many parameters that accompany the machine func- 
tioning (vibration, noise, temperature, pressure) must be en- 
sured in order to have as much complete and correct information 
as possible. 

The reliability of the supervising/diagnosis system must be supe- 
rior to the reliability of the measured system. 

The system has to be easy to use and maintain, and be resistant to 
dust, dump, and industrial liquids, low and high temperatures, 
and sometimes to radiation. 

The system must not be connected to the same energy sources as 
the supervised system but to special stabilized and protected 
sources. 

The cost of the supervising system, including its installation, should 
have reduced weight, usually under 10% of the cost of the su- 
pervised system. 

1.3.3 Stages of Technical Diagnosis 

Implementation on Machine Tools 

It is possible to establish technical diagnostics for complex systems such 
as machine tools only after knowing their dynamics and kinematics. By 
comprehending these elements, the parameters that have to be super- 
vised, the types of transducers able to be used, and also the points where 
they will be installed can be established. When processing the signals, it 
is necessary to keep in mind the normal and maximum admitted levels of 
the supervised parameters. Figure 1.6 presents a pyramid of the stages 
that must be crossed in order to conceive a diagnostic system for a given 
machine tool or manufacturing process. 

The sensible points of the supervised systems are: the signal pro- 
cessor , which has to extract from the raw signal the necessary data, and 
the diagnosis processor , which has to use the data in order to identify the 
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Figure 1.6 Pyramid of a proposed diagnostic system. 


state of the supervised system and to locate and isolate the fault. This 
should be done according to the specific algorithms of the diagnostic 
proceeding used. 

1.3.4 Efficiency of the Diagnostic Systems 

The term “efficiency” designates the ratio between the ensemble of the 
effects obtained as the result of an activity and the totality of effort 
(expenses) that this activity implies. In the case of technical systems, 
we can speak of the introduction efficiency of technical diagnosis. 

Concerning the expenses implied by the introduction of the diag- 
nosis system, the following must be taken into consideration. 

The supervising/diagnosing equipment has to be adjusted to the 
supervised parameters. 

The sizes and the complexity of the supervising system. 

The cost of the supervising system, installation, and maintenance. 

The cost of training personnel. 

Concerning the useful effects, these are manifested by the following. 

Significant increase of effective production time. It was noticed that 
for machine tools without monitoring only 10% of the theoretical 
production time could be reached, and for machine tools having 
supervising systems and diagnosis, the production time increased 
up to 60% from the theoretical production time. 

Reduction of accidental interruption times, and also decreasing of 
the time necessary for repairing. 


Decrease in the consumption of cutting tools and elimination of 
waste because of the supervising during the time of the cutting 
process, including wear and breakage of the tool. 

Sustaining the adaptive command for using the intensive cutting 
regimes in order to increase productivity. 

Recent studies have shown the economic efficiency of introducing 
diagnostic systems, estimated on the basis of the above criteria, increased 
on average from 2 to 3% of the value of the net production. Of this 
percentage, 65% is the effect of avoiding production stoppage, and 34% 
is the effect of reduction of maintenance and repairing costs [54, 108]. 


1.4 VIRTUAL INSTRUMENTATION FOR 

ESTABLISHING TECHNICAL DIAGNOSIS 

The analysis and measurement branch of industry is undergoing spec- 
tacular changes because of rapid developments in hardware and software 
technology. In research, designing, testing, measuring, and control activ- 
ities, PC computers are being used more and more, so companies pro- 
ducing instrumentation have reoriented methodologies and measurement 
equipment in order to better exploit the PC's hardware and software re- 
sources. On the other hand, the limits imposed by the rigid architecture 
of traditional instruments has generated, during this time, nonconcor- 
dances between the offer and the functionality request, namely, between 
what the instrument producers offer and what the user wants. 

In the last decade, the idea of combining a new measuring instru- 
ment programmable by standard PC computer has led to the creation 
of a new concept, that of “virtual instrumentation,” whose functions are 
defined by the user and not by the producer. This innovation became 
possible only after the appearance of digital instruments and the commu- 
nicating interface type GPIB (general purpose interface bus). Through 
this interface, these digital instruments can now be controlled by a pro- 
gram. The current generation of technical measurement offers more flexi- 
bility and performance since the instrument itself is built as a component 
of the PC. 

In 1986 National Instruments (Texas) Company launched the first 
release of virtual instrumentation software for engineering purposes, 
named Lab VIEW (Laboratory Virtual Instrument Engineering Work- 
bench) 1.0. This software uses a graphical programming language G in 
order to create programs shaped as block diagrams, without sacrific- 



ing any of the power of a traditional programming language. Lab VIEW 
uses terminology, symbols, and ideas familiar to researchers and engi- 
neers, based on graphical symbols rather than on textual language for 
describing programming activities. Lab VIEW combines the most recent 
operating system technology with specialized programming techniques 
(OOPT, object-oriented programming techniques) in order to obtain a 
simple and flexible operating environment. 

Under the slogan, “The software is the instrument,” the National 
Instruments Company identifies and structures the steps for use of the 
virtual instrument (Fig. 1.7): 

At the first level, the virtual instrument uses hardware controllers 
and interface type GPIB for connection with the programable 
instruments. 

At the medium level, the virtual instrument uses standard hard- 
ware architectures and corresponding drivers (specific driver 
packet containing libraries for the functions that operate the 
hardware). 

The superior level is represented by the software applications Lab- 
VIEW or LabWINDOWS. 

The Danish firm Bruel & Kjaer, specializing in measurement equip- 
ment for the vibroacoustic domain, in 1992 culminated their efforts of 
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Figure 1.7 Virtual instrument usage steps. 









integrating the computer in the measuring chains by launching the soft- 
ware Modular Test System Type 3538. The company, which has had 
exceptional experience in the field of signal processing with dedicated 
traditional equipment, created a virtual instrumentation that couples 
widely used classic instruments with a powerful software capable of pro- 
viding the instrumental functions that the user needs. Moreover, using 
the X-Window System and Motif application, the user can create and 
use his own instrument, front panels, and graphical interfaces. The per- 
formance of the virtual instrumentation launched by B & K is based 
on two special qualities: the precision of the digital/analogue and ana- 
logue/digital conversion and the processing power of the digital signal. 

Figure 1.8 presents the front panel of one of the most powerful vir- 
tual instruments offered by the B & K company, a double-channel spec- 
tral analyzer, Spectrum and System Analyzer Type 7627. It includes not 
only the measuring instrument but also the signal generator and real- 
izes a large range of procedures such as spectral analysis, measurement of 
frequency response, measurement of harmonic distortions, simultaneous 
measurements in time and frequency (wave shape, magnitude spectrum, 


• ; 3 '•! - • >y'- ' - f Brug* &. Kjaa r . Typft T6Z ? •- ’ Spectrum and 


11 id 


uose 


Configuration 


Setups 


Information 


Results 


Opti 




Help 


Dteptey 





and the phase, coherence, and correlation functions), and finally, testing 
and finding the measurement schema faults. 

The flow of information in this measurement and analysis system is 
indicated in Figure 1.9. The physical signals enter the instrument by the 
hardware module, which is controlled by a dedicated driver. The server 
for measurements controls this driver, conceiving an interface between 
the hardware and the virtual instrument. The virtual instrument takes 
over and processes the data on the basis of its functions, and the results 
are communicated to the operator by the front panel. 

The well-known German firm, Hottinger Baldwin Messtechnik, has 
merged with the Group Spectris concern, which produces control and 
measurement equipment and has introduced in its turn virtual instru- 
mentation as a working environment. The new generation of Hottinger 
equipment is adapted to acquisition, computerized processing, and anal- 
ysis of signals. Moreover, one of the fields of interest of the company is 
monitoring industrial processes. 

A virtual instrument is defined as being a software/hardware in- 
terface that is added to the computer in such a way that the user can 
interact with it as she interacted with the traditional instrument. This 
instrument can accomplish a compulsory family of functions. 

1. The function of data gathering is executed through a data ac- 
quisition board, connected straight to the processor bus; the memory 
registers from this board are accessible to addresses from the I/O space 
of the computer memory. The computer controls the data acquisition 
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Figure 1.9 Flow of information in measurement and analysis system. 















board, and the data transfer between the board and processor is made 
using hardware interruptions on DMA dedicated channels. The problem 
that arises when performing this function is that of the speed of data 
acquisition: the computer can not compete with the dedicated proces- 
sor of traditional instruments, even though the computer is completely 
occupied with solving the operation executed by the virtual instrument. 
The computer interruptions have priorities, and they are stored in a 
waiting file that results in an increase of the working time. Traditional 
instruments can usually gather data with a speed of 10 GHz, but an 
acquisition board does not surpass the speed of 1 MHz. The solution for 
this problem is to endow the acquisition board with counter /timer chips 
that use internal frequencies up to 10 MHz. In this way, the data acqui- 
sition will be made correctly by the virtual instrument, which means to 
the speed specified by the user. The concurrence of tasks that have to be 
solved by the computer’s processor eventually may delay data storage 
and/or the presentation of the results. 

2. The data analysis and control functions are completely carried 
out by the hardware already existent in the computer and by the soft- 
ware that, to a large measure, is already familiar to users. The virtual 
instrument uses, as do traditional ones, software modules from a large 
package, but the difference is that as traditional instruments close this 
software on its RAM memory, the virtual instrument holds the functions 
on the computer’s hard drive or on a floppy disk that can be installed 
on every computer. Moreover, many virtual instruments can coexist on 
the same computer, using the same display, independently or in direct 
relation with one another. 

3. The function of results display is another compulsory function. 
The existence of a driver with a graphical interface enormously dimin- 
ishes the handling and control of the application. The instrument is 
presented on the virtual panel display on the computer’s monitor, and 
may look like the front panel of traditional instruments. The virtual 
panel has in its background the software program, which means the 
instrument commands that are blended into that application, together 
with the acquisition routines, data analysis, graphical presentation, and 
finally the capacity to eventually store data/results in files. This con- 
struction of the virtual instrument makes possible, for the first time, for 
its functionality to be defined entirely by the user. The virtual instru- 
ment uses standard hardware architectures and corresponding drivers, 
which is a specific package with libraries that operate the hardware. 



In principle two programming methodologies exist for the use of 
virtual instrumentation: 

Graphic programming, such as Lab VIEW software 

Programming in a traditional language, which is the case of 
Lab Windows/ C VI (C for virtual instruments) software [160] 

Launched on the market during the second half of the 1990s, vir- 
tual instrumentation has seen great development and, nowadays, is a 
standard factor in testing and control. 



2 


The State of the Art of Research 
Concerning Vibroacoustic 
Diagnosis of Machine Tools 


Monitoring and diagnosis of machine tools and technological equipment 
on the basis of vibration and noise analysis has greatly expanded in 
the last twenty years. Although a relatively young field of research, this 
domain is marked by important theoretical and experimental findings 
that are selectively presented in this chapter. 


2.1 DIAGNOSIS OF ROLLING BEARINGS 

Bearings are a major source of vibration and noise in machine tool manu- 
facturing. Their vibration and noise level correlate with their functioning 
state, indicating the presence and eventually the evolution of some dam- 
age. Thus, interventions must be made in order to avoid serious damage. 

For bearings in a nonmounted state a sum of characteristic fre- 
quencies is highlighted: the rotation frequency of the bearing cage / c , 
the frequency of the ball running on the internal inner //, and the ball 
rotation frequency /&. These frequencies depend on the rotation fre- 
quency of the spindle f n . Table 2.1 presents the values of these frequen- 
cies for several types and sizes of bearings. It can be observed that the 
fc<fn<fe<fe< Si relation always exists. Maintaining the relation 


Table 2.1 Characteristic Frequencies for Bearings in 
Nonmounted State 


Frequency 

(Hz) 

Bearing 6206 

Bearing 6209 

Bearing 7309 

fc 

0.391 f n 

0.402 /„ 

0.410 f n 

u 

3.132 f n 

3.621 /„ 

4.506 f n 

fi 

4.868 f n 

5.379 f n 

6.497 fn 

h 

2.191 /„ 

2.460 /„ 

4.022 fn 


between the characteristic frequencies shows that the shape and sizes of 
the rolling bodies influence these frequencies only in small measure [46]. 

The vibrations and noise produced by the bearings in a mounted 
state manifest themselves not only by direct effects (acoustic radiation, 
vibration), but also indirectly (on the level of other elements, on the 
route source-receiver). The main causes that lead to an increase of the 
noise and vibrations on the mounted bearings are: 

Changing the position of the rolling bodies in the charged (respec- 
tively, uncharged) zone of the bearing, taking into consideration 
the elastic contact deformations and the working play 

The nonuniform movement of the rolling bodies, which leads to 
friction and collisions of the rolling bodies with the rings and 
the cage 

The contact that appears on the rolling movement on the surfaces 
with dimensional, shape, and position deviations (deviations of 
diameter, axial or radial runout of the rolling paths, eccentricity, 
ovalization, polygonality, waviness) 

Passing of the rolling bodies over the impurities from the contact 
surfaces or over damages or local deterioration (pitting or peel- 
ing, abrasive wear, prints) 

Table 2.2 synthetically presents relations concerning the position 
of the peak in the frequency spectrum for these damages; the data have 
been considered for the bearing with turning internal ring ( fi = f n = 
n/60) and fixed external ring. The notations represent: / c *, turning fre- 
quency of the cage function of the internal ring; / ce , turning frequency 
of the cage function of the fixed external ring; and /^ e , the turning fre- 
quency of the rolling bodies function of the fixed external ring. 


Table 2.2 Relations Concerning Position of the Peak in the Frequency 
Spectrum for Damages at Mounted Bearings 


Element 

Deviation 

Position of 
peak in 
spectrum, Hz 
(linear effects) 

Position of 
peak in 
spectrum, Hz 
(nonlinear effects) 

Internal 

ring 

Radial runout 
Waviness 

fi 

mzfci ± pfimz 

'ffl'fci T fee 

External 

Singular damage 
(pitting) 

Waviness / singular 

fa ± pfi 
mzfce 

(m± 1 )fce 

ring 

Rolling 

body 

damage 

Diameter variation 
Waviness 
Singular damage 

mfce 

2 mfbe ± pfce 

2 mf be ± pfce 



The frequency that corresponds to the succession in-out of the 
rolling bodies from the singular damages can be estimated with the 
equation: 


ft = 7 rD m n 




where b is the damage length on the circumference of the rolling path. 

Monitoring of the vibration level can indicate nearly all types of 
mechanical faults of the dynamical equipment. An increased vibration 
level signals imbalance, nonalignment, general wear, or play between 
subassemblies; a high shock impulse level signals a nonalignment, an 
incorrect lubrication, or the existence of an overload on the bearing. 

A rapid and efficient method of monitoring the functional state 
of the bearings, named the shock impulses method, has been patented 
by the Swedish company SPM (Shock Pulse Meter). The principle and 
mathematical apparatus of the method are presented in great detail in 
this paper. The level of shock impulses of a bearing depends on the size 
of the bearing, the rotation speed, and the working conditions. 

The SPM equipment can indicate the limit values for good, re- 
duced, and unsatisfactory functioning conditions for any kind of bearing. 
Figure 2.1 presents a monitoring system for simultaneous supervision of 
16 measuring points. The scavenged frequency domain is between 10 and 
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Figure 2.1 Monitoring system. 


10,000 Hz, corresponding to the ISO 2372 recommendations (equivalent 
to BS4675 and VDI2056). The measurement is made on (a) channel 
ALERT, which is mainly used to monitor significant increase of the vi- 
bration level due to wear, a channel with a discrete adjustable level from 
25 to 75 dB sv ; and (b) channel ALARM, which is used to warn and to 
stop voltage feeding of the monitored machine. Figure 2.2 illustrates the 
estimation graphics of the functioning state of the bearings. 

Frequency analysis of the vibration signal is the most efficient diag- 
nostic method for bearings. During the first investigation, the vibration 
level in the large band (third octave) or narrow band (10%, 3%, 1%) 
can be analyzed. This way the power of the signal can be appreciated, 
and the frequency peaks can be located in order to diagnose them (see 
Table 2.2). This process is complicated by the fact that the frequency 
peaks do not manifest themselves in a uniform way on the whole spec- 
trum: at low frequencies they can be masked by white noise, and at high 
frequencies the danger of resonance appears. It must be noted that the 
position in the frequency spectrum presented previously depends on the 
bearing turn, because the position in the frequency spectrum previously 
presented is independent of turn (Figure 2.3). 

Because of the presence in the spectrum of numerous harmonics 
that disturb the signal interpretation (some of them hide significant 
frequencies), Cepstrum analysis is recommended in bearings diagnosis 
[19]. This way it is possible to detect the periodicity in the spectrum 
(families of harmonics and lateral bands), the analysis also being in- 
sensible to perturbations that appear on signal transmission. Detection 
by Cepstrum analysis of damage in a large-sized bearing (Figure 2.4) 
highlights the presence in the spectrum of the frequency of the pass of the 



Figure 2.2 Estimation of the functioning state of the bearings. 

rolling bodies over a pit from the external ring (FTBO), and also its first 
harmonics. 

Frequency analysis has given birth to a spectacular method of in- 
vestigation of the vibration signal, a method of comparing the spectra. 
The Danish company, Briiel & Kjaer, launched the method and made it 
operational. The differences between the measured current spectrum and 
the spectrum templates come from the boarding of the reference signal. 
If the current spectrum surpasses the tolerance template, a warning sig- 
nal is emitted and, if the limit template is surpassed, the emitted signal 
is one of damage and in parallel the machine is stopped. The method 
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Figure 2.3 Analysis in frequency — spectrum. 


can be completed by an analysis of the tendency of the frequencies that 
bypass the tolerance templates, the analysis that approximates the time 
until the final damage. Figure 2.5 illustrates the use of this method as a 
diagnostic of large-sized bearings: (a) comparing the current spectrum 
with the tolerance template; (b) graphic highlighting of the spectrum 
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Figure 2.4 Detection of damage on bearings by Cepstrum analysis. 
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Figure 2.5 The method of spectrum comparison: (a) comparing the current spectrum with the tolerance 
template; (b) graphic highlighting of the spectrum displacements on narrow frequency bands; (c) selection of 
the alarming bypasses of over 6 dB; (d) drawing the tendency graphic for one of the selected frequencies. 





displacements on narrow frequency bands; (c) selection of the alarming 
bypasses of over 6 dB; and (d) drawing the tendency graphic for one of 
the selected frequencies. 

Another powerful method of frequency analysis, useful in diagnosis 
of bearings and also of gears and hydraulic elements is frequency analysis 
of the evolute of an intermediate signal from the time domain. The sure 
advantages of this method consist of highlighting the periodicity of the 
peaks corresponding to the damages from the bearing, their separation 
from the signals with close frequencies, or from the frequencies excited 
by resonance. Dr. D. Pupaza used this method with good results in the 
diagnosis of roller bearings, in the frame of a specialization at the Univer- 
sity of Leicester (United Kingdom). He applied the method to diagnose 
the bearings of the main shaft of a NC milling machine (Figure 2.6). 
The tests were made at idle running (2400 rpm) to minimize slippage 
and contact angle variation. The results obtained at the introduction of 
damage on the internal ring confirmed the increase in amplitude at the 
frequency estimated by calculus [126]. 

Acoustic signals are used for the diagnosis of bearings in the non- 
mounted state. Because of the reduced noise general level, the diagnosis 
of the bearings mounted in an assembly that contains other sources 
presents serious difficulties. The bearing under research is placed on a 
stand having very quiet characteristics, located in an anechoic chamber 
having inferior minimum frequency under 250 Hz. The placement of the 
free field microphone is indicated in Figure 2.7. 

In his doctoral thesis, D. Radauceanu [128] analyzed the problem 
of noise dependence of the radial bearings upon the manufacturing and 
functioning conditions, and highlighted some interesting general aspects: 

A correlation exists between the noise and vibration levels of the 
same bearing under the same measuring conditions, and is evi- 
dent in the frequency analysis of large band and less perceptible 
to frequency analysis of narrow band. 

The shape deviations of the rings and rolling bodies (microgeom- 
etry, waviness, roughness) negatively influence the noise and vi- 
bration levels. 

Ball bearings are quieter than roller bearings (5 to lOdB). 

A watertight bearing, lubricated with low viscosity grease has a 
noise level smaller than an open one (3 to 4dB). 

Mounting a roller in a soft gearing box (i.e., silicone rubber) or in 
a fixed position, reduces the noise with 6-10dB. 
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Figure 2.6 Diagnosis of main spindle bearings of a milling machine using frequency analysis of the 
evolute of an intermediate signal. 





Figure 2.7 Diagnosis of main spindle bearings using spectrum comparison. 


Load influences the noise level of the bearings very little, but the 
rpm considerably influences it (the acoustic pressure is approxi- 
mately proportional to the square of the turn). 

No relevant differences of the noise level were observed between 
lubrication with low viscosity level and oil lubrication [128]. 

2.2 DIAGNOSIS OF SLIPPING BEARINGS 

Slipping bearings are less frequently used in machine tool manufactur- 
ing. However, they may determine supplementary perturbations in the 


functioning of kinematic chains, usually with minor effects. The axle 
bearing in this case constitutes the main cause of the noise and vibra- 
tions. If radial plays in this kinematic couple also exist, they generate in 
spectrum peaks of the frequency of the axle f Q and of the subharmonics 
of this frequency (/ c /2, f a / 3). 

If the following phenomena are present: insufficient feeding with 
lubricant of the bearing, overloading of the axle (forces and moments), 
changing the rotation direction of the axle, and instability of the lubri- 
cant film at high rpm, then these lead to a thinner lubricant film and 
the occurrence of the “fiddlestick” phenomenon, which produces a noise 
with spectrum peaks at the rotation frequency of the axle f Q and at the 
harmonics of this frequency (2/ 0 , 3/ 0 , . . . ). 

The diagnosis of these bearing types is chiefly carried out by su- 
pervising the noise and the acoustic emission. The acoustic signal is 
detected, amplified, and filtrated by filter pass-up and then analyzed in 
the time domain. The averaging of the signal in amplitude (RMS) and 
the recording of the impulses that pass over a limit value, adjustable, of 
the amplitude are executed [39]. 

The change from the hydrodynamic friction to the mixed frequency 
is indicated by the occurrence of some peaks in the domain of the high 
frequencies, emitted as a result of going beyond the carrying capacity of 
the bearing. The frequency analysis puts in evidence the phenomenon 
through a dominant in front of the rotation frequency of the shaft. 

The spectral density of the impulse function can correlate very well 
with the functioning state of the bearing, as highlighted in Figure 2.8. 
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Figure 2.8 Spectrum of the impulse density for a hydrodynamic bearing. 



The upper part of the figure presents the spectrum of the impulse density 
for the hydrodynamic bearing in perfect condition, and the lower part of 
the figure presents the spectrum of the same bearing, this time damaged. 

In order to reduce the noise and vibrations generated by the slip- 
ping bearings, some measures can be taken: 

Choosing an eccentricity relatively bigger than an admissible value 
e = 0.2-0. 3, ensuring the stability of the axle 

Ensuring a correct lubrication 

Modifying the shaft rpm in order to avoid functioning at the critical 
rpm or double of this critical rpm 

Equilibrating the shaft in bearings, by the mean of the marks that 
are mounted on the axle 


2.3 DIAGNOSIS OF GEAR WHEEL 
TRANSMISSIONS 

The gears constitute one of the most important sources of noise and vi- 
bration in the machine tool structure. The gears 5 vibroacoustic behavior 
becomes more defective as the transmitted power and the rpm increase, 
and with the decrease in the size of gears. The parameters that deter- 
mine the vibroacoustic behavior of the gears can be grouped depending 
upon their nature: 

Constructive parameters (module, inclination angle, covering ex- 
tent, number of teeth, material, etc.) 

Technological parameters (the class of precision, execution errors, 
roughness, etc.) 

Functional parameters (rotation speed, load, moment, existence of 
lubricant) 

Table 2.3 synthetically presents the influence of the constructive 
parameters on the vibroacoustic level of the gear wheel transmissions. 
It should be noted that the main factor influencing the vibroacoustic 
level of the gear wheel transmissions is the teeth inclination followed 
by the stiffness of the gear wheel material. The gear's behavior on the 
dynamic regime is especially influenced by the parameters belonging 
to the two groups. The execution and mounting errors determine sup- 
plemental relative movements of the wheels, which superimpose those 
determined by the variable stiffness of the gear. The most important in- 
crease of the noise level (up to 20 dB) is determined by the shape errors 


Table 2.3 Influence of Constructive Parameters on Vibroacoustic Level of 
Gear Wheel Transmissions 


Analyzed parameter 

Analyzed 

values 

domain 

Influence of 
parameter value 
increase on noise level 
(dB) 

Tooth inclination 

0-40 

420 

Material stiffness 


412 

Tooth flank 

0 - 0.02 mm 

46 

Covering extent 

1 - 2.2 

43 

Module 

1-12 mm 

42 

Number of teeth 

z-2z 

43 

Specific deviation of profile 

0.092-0.883 

41 

Width of tooth 


Low effect 

Tooth bulging, stiffness, 


To large direction 

mass of frame box 


errors of tooth 


of the profile. They imply increased frequencies that can coincide with 
the natural frequency of the transmission. The decrease of the working 
surfaces' roughness (by shaving, grinding) leads to a decrease of 3 to 
4 dB of the noise level. If the relative error of the base step is positive, 
the gearing begins and ends beyond the theoretical gearing line; thus, 
the entrance and exit of the tooth from the gearing is made in shock 
conditions. The rotation error cumulates the effect of deformations and 
the errors of the teeth, providing complete indications on the mechanism 
of the excitation. 

Concerning the exploitation parameters, experimental research has 
shown that the noise and vibration levels increase at the same time as the 
load and angular speed. Doubling the loading moment has the effect of 
increasing the noise level by 3 to 4 dB. On the other hand, the doubling 
of the rpm produces an increase of 6 to 7 dB, at the same loading. The 
influence of the oil viscosity on the vibroacoustic behavior of gears is 
negligible [106]. 

The analysis in the frequency domain is the most adequate tech- 
nique for processing the vibroacoustic signal coming from the gearing, 
because the modifications of some spectral components correlated to the 
evolutionary damages that do not sensibly affect the global vibration 
level are detectable in the frequency spectrum, and the modifications 
of the spectrum provide important information concerning the causes 






of the damage. Usually, the frequency spectrum of the noise and vibra- 
tion from gearing presents peaks at: (a) the gearing frequency and its 
harmonics, (b) secondary components, (c) lateral bands, and (d) low 
harmonics of the rotary frequency. 

(a) The harmonics of the gearing frequency: the gearing frequency 
can be calculated as the multiplication between the rotation frequency 
of the shaft (/ = n/60) and the number of teeth of the wheel 2 , so f a = 
/ * z. The spectrum peaks corresponding to the gearing frequency and 
its harmonics are caused by the shocks, which arise in gearing because 
of the deviations from the ideal profile (Fig. 2.9). These deviations occur 
either during manufacturing or by the nonuniform wear of the flanks, 
which is larger on both sides of the rolling circle. The wear presence is 
more evident to the superior harmonics. 

(b) The secondary components (accidental) occur in the same con- 
ditions as the gearing frequency and its harmonics, but correspond to 
a different number of teeth, and these secondary components are due 
to kinematic errors of the generating kinematic chains of the gear- tooth 
cutting machine. The most important of these errors is the gear-tooth 
cutting error on the rotation of the gear-tooth cutting table. If the wheel 
has a number of N 0 teeth on the functioning of the gear wheel on the 
table activated by it, a frequency will appear corresponding to N 0 , as 
the harmonic integer multiple of the rotation frequency. This component 
is easily recognized because of its property of being independent by the 
increase of the load. 
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Figure 2.9 Gearing frequency harmonics. 



(c) The lateral bands around the gearing frequency and its har- 
monics can be explained by the amplitude modulation or the nonuniform 
vibration frequency of the gearing, caused by some local damage. For 
example, the radial runout of the gear or the eccentricity of the gear 
wheels’ sustaining shafts make the division circles of the two gearing 
wheels not remain tangent during the whole gearing time. Because of 
this radial runout, the low frequency corresponding to the rotational 
movement with A frequency superimposes a relatively high frequency 
vibration corresponding to the gearing process / a , resulting in a vibra- 
tion modulated in amplitude. In the limit case of a 100% modulation, 
the fi component disappears, and the lateral frequencies characteristic 
to modulation remain: /ii n f = f a + /i and /i sup = fa + fi- Because 
the modulation process is distorted by manufacturing errors and 100% 
modulation is not realized, peaks arise in the spectrum for frequencies 
fj inf = fa - j * fi and fj sup = fa + j * /l, with j = 2, 3, 4, . . . . Also, 
every fluctuation in the tooth load will have the tendency to generate 
a variation of the vibration amplitude, so an amplitude modulation ap- 
pears. At the same time, those fluctuations also determine fluctuations 
of angular speed, so a frequency modulation appears. The two modula- 
tions determine an increase of the lateral bands of frequency amplitude, 
their distance in the spectrum being equal to the modulation frequency. 
This contains important information on the source of the modulation 
effect. Each modulation with a pure frequency tends to give a family of 
lateral bands, and in the case where the modulation is not pronounced, 
this can be represented by two or three groups of lateral bands. Under 
these conditions, an additional effect is to increase the number of lateral 
bands. 

(d) The low harmonics of the rotation frequency correspond to 
the cumulated impulses that repeat with each rotation of the wheel. 
The modulation in amplitude or frequency has the tendency to give 
symmetric signals. Any asymmetry of the signal can be interpreted as 
an additive signal that also gives a number of harmonics of the respective 
frequency. 

In addition to the mentioned frequencies, other significant compo- 
nents may appear in the spectrum, independent of rotation, correspond- 
ing to some self-modes of vibrations for the gearing, shafts, gear box, and 
the like. Also, the variation function of time of the rigidity of the tooth 
gearing has a specific action. The parametric vibrations of the gearing are 
characterized in this case by the existence of many domains of instability 
to frequencies equal to rational fractions, multiples, or submultiples of 
the self-frequencies of the system. A resonance phenomenon can appear 



for the same reason, when the excitation frequency f e corresponds to 
one of the first subharmonics of the fundamental self- frequency fa. 

In conclusion, the frequency spectrum of the gearing transmission 
permits an evaluation of the mechanical state of the transmission to be 
made. In order to identify the sources of noise and vibration, a correlation 
of the frequencies corresponding to the maximums of the spectrograms 
with the frequencies of the possible sources from transmission is made, 
calculated on the basis of their constructive and functional parameters. 

Cepstrum analysis is a powerful and useful instrument of investi- 
gation for gear transmissions. It allows detection of periodicity of the 
spectrum and the precise determination of the frequency of the signal 
modulation. In the gear diagnosis domain, the Danish company Briiel & 
Kjaer has had notable results beginning in the 1980s. An excellent ex- 
ample is the diagnosis of some problems that occur in the differential 
mechanisms in trucks produced by Raba (Hungary) in collaboration with 
Briiel & Kjaer. 

The test consists of rolling the differential mechanism on a stand 
until the damage occurs, simulating diverse loading (resistant moments). 
The gear vibrations are permanently monitored, making it possible to 
stop the test at any moment in order to study the damaging mechanisms 
before the progress of the damage makes it impossible to recognize the 
source. The initiation of the damage is recognized by the increase of the 
characteristic frequencies and of the lateral bands around these charac- 
teristic frequencies (the Cepstrum analysis method) [14]. 

The gearing frequency is sufficiently low, so many harmonics of 
this frequency are included in a spectrum of 2 kHz. The monitoring of 
this zone has indicated the increase of the second and third harmonics’ 
amplitude because of the debut of the pitting phenomenon to the gear 
teeth. As the phenomenon progresses, the lateral bands around these 
harmonics increase. Figure 2.10 presents recordings obtained after the 
damage due to pitting has put gear functioning into danger. In Fig- 
ure 2.10a, the gearing frequency and its harmonics are evident, and a 
magnification in the zone of interest (the second harmonic, Fig. 2.10b) 
shows an excessive increase in the lateral bands, at a distance of 12.05 
Hz. The point of occurrence of the lateral bands can not be exactly pre- 
dicted but mathematical modeling of the damages and the perturbations 
introduced by these damages allows prediction. 

By processing the signal by Cepstrum analysis (Fig. 2.10c), the pe- 
riodicities become evident, cumulated in a single band that corresponds 
to the frequency of 12.05 Hz. The analysis also indicates the occurrence 



Figure 2.10 Recordings obtained after damage due to pitting that puts a gear function into danger: (a) gear- 
ing frequency and its harmonics; (b) magnification of interest zone; (c) Cepstrum analysis; (d) occurrence of 
wear and periodicities of the lateral bands. 









of wear for the gear crown, and, in this case, the periodicities of the lat- 
eral bands are grouped to 1.7 Hz (Fig. 2.10d). This highlights the fact 
that diagnosis of gear damage is possible by correctly processing the 
vibration signal, using methods from the domain of frequency analysis. 

The Department of Machine Parts and Mechanisms of the Poly- 
technic Institute of Iasi has developed some original research on the 
noise and diagnosis of gearing. V. Merticaru (1971) tackled the problem 
of dependence of the cylindrical gear wheels on manufacturing and work- 
ing conditions [106]. In his doctoral thesis (1987), B. Dragan continued 
this research with his contribution to noise and vibration reduction in 
the gearbox of machine tools. His research focused on identification of 
critical natural frequencies and of different kinds of damage in the gear- 
box of a milling machine, the kinematic scheme of which is presented in 
Figure 2.11 [39]. 

To diagnose damage, the determination of natural frequencies, 
modulation frequencies, and periodicities from the frequency spectrum 
has necessitated the building of a complex schema for analysis of the 
vibroacoustic signal, an analogue type schema having compression pos- 
sibilities, Zoom analysis, and Cepstrum analysis (Figure 2.12). 
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Figure 2.11 Schematic of gearbox of a milling machine. 



Analyzer 

B&K201 


Figure 2.12 Zoom analysis and Cepstrum analysis. 

Initially, reference spectra and “Cepstrums” for the gearbox, ac- 
cepted as functioning in good condition, have been recorded. Modifica- 
tions of the spectrum on the gearing frequency and its harmonics have 
been observed by introducing a damaged gear wheel into the gearbox 
(position 4 in the kinematic schema); in addition, the global noise level 
increased by 3 dB. Using Cepstrum analysis, an increase of approxi- 
mately three times over the reference level of the first harmonic, corre- 
sponding to the rotational frequency of the fault wheel (radial runout) 
was observed. Thus the advantages of Cepstrum analysis were high- 



lighted, advantages related to a special sensibility and a simplification 
of data interpretation, because in cepstru the periodical peaks from the 
frequency spectrum are cumulated. 

2.4 DIAGNOSIS OF BELT TRANSMISSIONS 

Belt transmissions generally produce a reduced noise and they even have 
the quality of attenuating the vibrations produced by other sources from 
the kinematic chains. The three main causes that lead to the occurrence 
of vibroacoustic phenomena in this case are: the variation of the friction 
forces on the relative slipping of the belt on the wheel, some aerodynamic 
phenomena, or imperfections on the belt joint. In order to decrease the 
intensity of such phenomena, the following are recommended: 

Use of plastic or rubber plated belts. 

Belt stretching has to be optimum. 

Belt pulleys have to be made from nonmetallic materials with high 
damping capacity. 

Use of trapezoidal belts instead of flat ones. 

The damaged driving wheels produce in spectrum components to the 
rotation frequency and its harmonics, and damaged belts generate spec- 
tral peaks on the f c = vxj L frequencies, where v is the speed of the belt, 
x is the number of belt wheels, and L is the length of the belt. The belt 
damage can be located with the aid of a stroboscope, fixing the image 
by synchronizing at the f c frequency. 

2.5 DIAGNOSIS OF DRIVING 
ELECTRICAL MOTORS 

The driving of the kinematic chains of machine tools is mainly realized 
by asynchrony or direct current electric motors. This section deals with 
driving by asynchronous triphased motors. Their functioning in a sta- 
tionary regime is accompanied by phenomena that are constituted of 
forced vibrations and noise, influencing the elastic system of the driven 
kinematic chain. 

2.5.1 Vibroacoustic Phenomena of 
Mechanical Nature 

These phenomena are determined by the imbalance of the rotational 
parts, by the dynamic efforts variable in time (raised due to an incorrect 



execution and/or mounting), by the vibration of the brushes on the 
alternator, and so on. For example, the brush noise depends on their 
quality and on the quality of the sliding surfaces, on the holders’ guiding 
way, and on the contact pressure. The noise increases at the same time 
as the rpm n of the rotor and as the number z of the collecting lamella 
increase. In spectra, the brushes’ noise occurs to the fi = n * z/6 0 
frequency and superior harmonics. 

In order to identify mechanical noise, which is always combined 
with magnetic noise, the motor is disconnected from the power source; 
the electrical components of the noise disappear immediately, allowing 
study of the mechanical ones. The mechanical noise can be reduced by a 
static and dynamic equilibration of the rotor, using the slipping bearings 
instead of rolling bearings by the optimization of the brushes/collector 
pressure [4, 11]. 

2.5.2 Vibroacoustic Phenomena of 
Magnetic Nature 

Magnetic noise of electrical engines is due to the periodic magnetic forces 
from the air gap of the dynamo, the magneto-driving forces, or parasite 
magnetic fields. In the case of asynchronous motors, the angular fre- 
quency of the parasite magnetic fields is superimposed on the harmonics 
of the main field, which leads to an important increase of the magnetic 
noise level. The parasite alternative magnetic forces arise when the sta- 
tor is not circular, the rotor has no axial symmetry, or is not centered. 
The notches of the rotor and stator distort the magnetic field by con- 
centrating the lines of the flux density in the air gaps from the teeth 
(Fig. 2.13). Under these conditions, a complex vibration is born, whose 
frequency is given by the equation: 

1 — s 

( nR s ± k e ) =t hi 

V 

where cj is the frequency of the electrical network (50 Hz), R s is the rotor 
notch number, k e represents electricity coefficient (0 for static eccentric- 
ity, and 1, 2, 3 for dynamic eccentricity), s represents the slip, and p rep- 
resents the number of pairs of poles, n = 1, 2, 3, ... , k\ = 0, 2, 4, 6, . . . . 

Reduction of magnetic noise can be obtained by the inclination of 
the rotor notches or even by using rotors with double-inclined notches, 
executing some symmetrical electric reeling, and by reducing the air gap. 
The reduction of the magnetic field intensity by decreasing the excitation 




Figure 2.13 


Magnetic flux 


in an electric motor. 


current leads to substantial reduction of noise; however, this measure is 
disadvantageous because it implies a substantial increase of motor size, 
for the same power. 

Analysis in the frequency domain and also Cepstrum analysis allow 
for the detection of mechanical and magnetic damage, which occur dur- 
ing the functioning of the asynchronous motor. Figure 2.14 (left) presents 
a recorded example of the base spectrum (low frequency) of the function- 
ing of a mechanical damaged motor (rupture in the motor). The effect 
becomes evident by the amplitude increase to the rotation frequency; 
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Figure 2.14 Detection of mechanic and magnetic damages: (left) analysis in frequency domain; (right) Cep 
strum analysis. 






Table 2.4 Synthesis of Direct Current and Asynchronous Motor Diagnosis Research 


Vibrations cause: 

Arise 

frequency 21 

Dominant 

plane 

Observations 

Lack of rotor axle poise 

1 X fr a 

Radial 

Unbalanced type can be established analyz- 
ing signal phase (static 0°; couple 180°; 
dynamic 0-180°) 

Bent shaft or angular 
unaligned 

1 X, 2 X f r 

Axial 

180° axial 
0° radial 

Parallelism error 

1 X, 2 X fr 

Radial 

180° axial 
180° radial 

Mechanical plays 

1 x, 2 x, 3 x, 
4 X f r , . . . 
and 0.5 x, 
1.5 x f r , . . . 

Radial 

Large number of harmonics and half 
harmonics 

Damages of rolling bearings 

Wide frequency 
domain (1 to 
20 kHz) 


Resonance frequencies are excited by impact 
with local damages from rolling paths 

Incorrect lubrication (shocks) 
of slipping shafts 

0.43 T 0.48 x f r 

Radial 

Slipping bearings are used on very high power 
motors 

Nonrigidity of stator holder, 
nonequilibration of phase 
reels’ resistance, stator 
exfoliation 

2xf R 

Radial 

These damages are hard to differentiate by 
vibration analysis, but it is significant that 
they appear not only without, but also 
under load 

Dynamical eccentricity, 
fissured or broken 
stator bars 

1 x f r and 2 x f a 

Radial 

Bent, worn-out, or deformed rotor by local 
heating 

Slipping frequency is usually low 


a / r represents the rotor rotation frequency; /# represents the electrical network frequency; f a represents the slipping frequency. 






at the same time, lateral bands arise around this amplitude, produced 
by the increase of magnetic noise due to the slipping increase. Cepstrum 
analysis [Fig. 2.14 (right)] confirms the vibroacoustic phenomenon ob- 
served in the frequency domain. 

2.5.3 Vibroacoustic Phenomena of 
Aerodynamic Nature 

The aerodynamic noise is caused by the forced circulation of the cooling 
air from inside the electric machine, and also by the resonance of the 
air from the gaps and orifices of the rotor and stator. These phenomena 
have a harmonic character and can be identified in the spectrum on the 
frequencies f a = k * n * z/6 0, where z is the number of pallets of the fan 
(for the noise of the air of forced cooling) or the number of rotor notches 
(for the rotor noise), and n is the motor rpm, k = 1, 2, 3, . . . . 

Decreasing the aerodynamic resistance of the ventilation network 
can reduce the aerodynamic noise. Thus, aerodynamic shapes must be 
designed for the ventilation radial orifices of the stator and rotor; the 
fans must be built with aerodynamic pallets; and parts to direct the 
air must be designed in order to avoid the bends and narrowing of the 
cooling air jet. In order to avoid formation of the acoustic resonators, 
the existing gaps and orifices in the rotor and stator of the electric 
machine should be filled with phonoabsorbing material. For high rpm 
motors, noise attenuators for air admission and evacuation channels can 
be used. 

Table 2.4 represents a synthesis of research in the domain of direct 
current and asynchronous motor diagnosis currently used in industrial 
drives. The superimposing of all the excitation sources presented leads 
to the occurrence of numerous spectral peaks. The experimental research 
on the vibroacoustic behavior of the usual triphase asynchronous motors 
(0.6 to 7.5 kW) has shown that the zone of interest of the frequency 
spectrum is between 16 and 530 Hz. The same domain is significant also 
in the case of direct current motors that are used to drive the machines 
and equipment. 


2.6 DIAGNOSIS OF SOME DEVIATIONS OF 
THE ROTATIONAL PARTS 

The dimensional, shape, and position deviations of the motion rotational 
parts constitute an important source of vibrations frequently encoun- 


tered on machine tools and equipment. They become more and more 
destructive with the increase of the machining rpm. The causes of these 
deviations can be execution and/or mounting errors, nonhomogeneity of 
the materials used (inclusions on the cast parts or density variations on 
the cast or forged parts), the asymmetrical configuration of some pieces, 
hydro- or aerodynamic phenomena, wear phenomena, and so on. 

2.6.1 Imbalance (Lack of Poise) 

Imbalance, or lack of poise, arises when the rotational and mass centers 
of a rotational part are not coincidental. If imbalance can be separated 
into a single plan, this is denominated static imbalance, and if this arises 
in many more planes, dynamical imbalance results. 

In the case of static equilibration the two centers are made to coin- 
cide by adding or eliminating a mass on the rotational plane [the case of 
the grinding wheels, flying, and belt pulleys as presented in Figure 2.15 
(left)]. In the case of dynamic equilibration, not only the centrifugal 
forces but also the coupling of these forces must be taken into consider- 
ation; this coupling produces an angular movement of the main inertial 
axis compared to the rotational axis. By equilibration, these two axes 
must be brought into coincidence [the case of rotors, Fig. 2.15 (right)]. 
Because of the centrifugal forces created by the unbalanced masses the 
following occur: dynamic forces appear in bearings formed from a pri- 
mary component due to the nonequilibrated mass of the bearing plane 
and a secondary component due to the unbalanced coupling. The im- 
balance on the bearing level can be defined with the help of a turning 
vector whose magnitude is given by the magnitude of the resultant force 



rotation axis 





inertial axis 




Figure 2.15 Equilibration: (left) grinding wheel; (right) rotors. 



and whose angle is given by its action direction. This direction may be 
determined precisely by comparing the phase of the signal periodically 
created by the vibration of the bearing with a reference periodical signal, 
obtained in a known position of the rotational part. 

The imbalance of the rotational parts produces vibrations on the 
rotation frequency fi = n/60 and on the superior harmonics of this fre- 
quency, having high amplitudes especially in the radial direction. The 
amplitude of these vibrations increases at the same time with the in- 
crease of rpm. The imbalance also can be detected by measuring the 
vibration movement phase: the perturbation force is a rotational vec- 
tor, so the phase will be dependent on the transducer’s location, as the 
amplitude of the movement will practically remain constant. 

2.6.2 Axial Imbalances 

The vibrations produced by axial imbalance are characterized in the 
spectrum by the existence of some significant peaks of the rotational 
frequency and especially of its second harmonic (2x). The ratio be- 
tween the amplitudes of the two components gives information on the 
degree of axial imbalance. The axial imbalance produces some high-level 
axial vibrations that can be perceived at both ends of the shaft with a 
dephasing of 180°, making them different from those produced by lack 
of poise, also located at the ends of the shaft. 

2.6.3 Dimensional and Shape Deviations 

These errors become significant on machine parts that have rolling con- 
tact (gears, bearings, ball screws, roller guidings, etc.). 

In the case of gear transmission, when the leading wheel has con- 
stant angular speed, it can be observed that the angular speed of the 
driven wheel is influenced by the transmission errors cumulated from 
the dimensional and shape errors of the gear geometry. As the transmis- 
sion error is a function dependent on time, each spectral component of 
this error will be determined by the force of interaction between the two 
wheels. 

The shape errors of the profile (Ef p ) lead to an increase of up to 
20 dB in the noise level of the gearing (Fig. 2.16), because their frequen- 
cies are high and can coincide with the natural transmission frequencies. 
These errors are signaled in the spectrum as harmonics of the gearing 
frequency and low harmonics of the rotation frequency [46]. 



2.6.4 Mechanical Plays 

The existence of the mechanical plays is highlighted by the presence 
of some strong directional vibrations, a characteristic that permits their 
separation from those produced by lack of poise and axial imbalance. If a 
harmonic load stresses an assemblage with plays, the answer of the struc- 
ture is no longer harmonic, even though the frequency remains periodic 
(Fig. 2.17). The detection of the mechanical plays from the joints and 
kinematic couplings of the machine tools, even in the incipient phase, 
is very important because their influence on the shape and dimensional 
precision and quality of the part surfaces is of first rank importance. 


2.7 DIAGNOSIS OF THE TOOL AND 
CUTTING PROCESS 

2.7.1 Diagnosis of the Cutting Tool 

Classifying criteria for evaluation of the integrity and wear of cutting 
tools have been chosen: (a) the moment of control, divided into offline 



Figure 2.17 Structure response when mechanical plays exist. 


methods or evaluation methods outside the cutting process, and online 
methods or evaluation methods during the cutting process; and (b) the 
manner of measurement, also divided into groups: direct methods, in- 
cluding measurement of the cutting tool (the width of wear on the main 
and secondary clearance faces, the width of wear on the transversal 
edge, the depth of the wear crater on the rack face, and the width of 
the crater); and indirect methods that evaluate the wear characteristics 
of the tool by wear, Assuring, and breaking effects on the part or on the 
variation in time of some correlate measures with the cutting process. 

The offline evaluation of the tool state by direct methods is usually 
achieved by the tool's movement in numerical command to a point whose 
coordinates are known, followed by the cutting edge, rack, and clearance 
face’s palpation. The disadvantage of this method is that the tactile 
feelers can not perceive the integrity or the wear but only the palpation 
direction. 

Many more direct methods without contact used to evaluate tool 
status exist: 

Perception of the integrity of the cutting edge (drills, taps), using 
pneumatic transducers [142] 

Supervising the level of radioactivity of the tool cutting edge, pre- 
viously bombarded with hard particles [73] 

Wear and cracking evaluation using ultrasound sources 

Digitization of the wear surface image of the nonworn tool and its 
comparison, at time equal intervals, with the image of the real 
wear, calculating the distribution of the wear surface [97, 167] 



Present online methods for evaluation of the tool state cover over 
60% of total investigation techniques. They are less precise in determin- 
ing the amount of wear, but they can perceive the fissure and the crack 
of the cutting edge and do not interrupt the manufacturing process. The 
most frequently used methods are discussed in the following. 

2. 7. 1.1 Wear Control 

Used on the basis of processing time, it considers the estimated durability 
of the tool. 

2 . 7 . 1 . 2 Load Characteristics Evaluation of the 

Tool State 

The Skandvik Company produces tensometric jacks that must be used 
in the bearing of the lead screw, because of the special sensibility of 
the force on the lead direction, comparative to the tool wearing. The 
Prometec Company, Germany, produces special bearings for the main 
shaft, bearings containing tensometric marks mounted in special chan- 
nels of the external ring. The signal amplitude captured in this manner 
is proportional to the shaft load, implicit in the cutting tool. Figure 2.18 
schematically presents the formation of the signal for comparison and 
evaluation of wear and fracture of the tool in the cylinder bore manufac- 
turing process, using the Prometec technique (where Pq represents the 



Figure 2.18 Formation of the signal for comparison and evaluation of wear 
and fracture of the tool in the cylinder bore manufacturing process, using the 
Prometec technique. 



minimum load level during idle run, P m - m represents the minimum load 
level during work, P max represents the reference level for the limit wear, 
and Pf max represents the reference level for the tool fracture) [23]. 

2. 7. 1.3 Acoustical Emission (AE) Evaluation of the 

Tool State 

The acoustical emission signal can be correlated with specific cutting 
tool wear criteria. These criteria are multiple and depend on the pro- 
cessing type (roughing and finishing) and on the used tool. The most 
used criteria are wear on the clearance and/or rack face, the dimensional 
deviations of the part, and the roughness of the processed part. 

The characteristics of the AE signal modify at the same time as 
the wear; an important modification of the effective value of this signal 
indicates the reach of the wear limit. It is interesting to note that the 
cutting process can be usually continued, without rupture, but in very 
bad conditions of geometrical shape and quality of the surface. 

During the cutting process, the effective value of the AE signal 
has a specific evolution: progressive increase, followed by a substantial 
increase in the vicinity of the catastrophic wear phase (this is much 
more evident on tools equipped with metallic carbide plates), then an 
accentuated decrease until the moment of the cutting edge destruction. 
The same behavior is found both in the roughness and finishing cases. 

Figure 2.19 presents a recording of the effective value of the acous- 
tic emission signal on a milling machine with a mill having a single tooth, 
in three successive passes. The passes are separated by two zones of ab- 
sence of the AE signal; during the third pass the tooth has broken, a 
fact highlighted by a peak, and then an abrupt decrease of the effective 
value of the signal. 

Concerning the break of the cutting tool, the most frequent causes 
are the mechanical shock when meeting some hard inclusions in the 
material to be processed or successive mechanical shocks (on the turning 
of some polygonal and grooved parts), or the thermal shock as a result 
of interruption of the cooling liquid. Tool break is easily highlighted by 
supervising the effective value of the AE signal because of the transitory 
energy jump. A typical example of tool breaking, hard to highlight by 
other methods, is the breaking of the auger when boring (Fig. 2.20). The 
tool-part contact is well detected by supervising the effective value of the 
AE signal and is clearly separated from the breaking zone (between the 
dashed lines) and from the zone of ulterior rack of the tool. 



Figure 2.19 Recording of the effective value of the acoustic emission signal 
on a milling machine with a mill having a single tooth, in three successive 
passes. 


In conclusion, the AE method is applicable and already widely 
used for the detection of tool breakage, but monitoring tool wear in the 
same manner is a more complicated problem because of many factors. 
Finding the wear criterion that best correlates with the evolution and 
the value of some specific dimensions of the AE signal for roughing and 



Figure 2.20 A typical example of tool breaking. 



finishing operations in the case of different cutting processes remains 
another experimental study problem and requires the accumulation of 
an important database. 

2 . 7 . 1 . 4 Electrical Resistance Monitoring 

This monitors the tool-part contact or the thermoelectromotoric tension 
that arises during processing. The contact resistance in the tool-part 
zone is of 10~ 4 to 10 -2 compared to 10 2 to 10 3 fl, the resistance of 
the bearings and 10 to 10 2 fl, the electrical resistance due to the cooling- 
lubricating fluid. The fissure or fracture of the cutting edge produces an 
increase of up to 30% and then an abrupt decrease to zero of the voltage 
in the measuring circuit [122]. 

H. Wiele and P. Metz [165] highlighted the possibility of super- 
vising tool wear by estimating the thermal gradient. Placing two ther- 
mocouples under the cutting plate, they discovered that the difference 
between the contraelectromotor voltages of the two controlled points is 
proportional to the wear on the clearance face. The disadvantages of 
these two methods consist of the difficult emplacement of the transduc- 
ers, and also on the sensibility of the thermoelectromotor voltage on the 
variations of the tool geometry and of the parameters of the cutting 
regime. 

Recently, Romanian researchers have grappled with the domain of 
diagnosis of the cutting tool state in their doctoral theses. E. Carata 
[23] brought important contributions on the elaboration of a research 
methodology in order to determine online the tool wear on boring, by 
vibroacoustic methods. L. Bogdan [12] correlated the information on 
three parameters (noise, vibration, and temperature) in order to evaluate 
the state of the cutting tool. 

2.7.2 Diagnosis of the Cutting Process 

The cutting process generates elastic waves that transmit to the working 
piece and the tool and its support. An acoustical emission transducer, 
placed in the vicinity of the cutting zone, can detect them, and the 
perceived signal is a “signature” characteristic of the cutting conditions. 

At the beginning of the 1990s, T. Moriwaki [110, 111] proposed a 
supervising scheme of the acoustical signal in cutting, which today is 
widely accepted (Fig. 2.21). The inferior cutting frequency of the filter 
set used is 100 kHz, in order to avoid the vibroacoustic signals associ- 
ated with the system machine tool-clamping device-tool-part. The su- 
perior cutting frequency is 10 Mhz, in order to avoid electrical noise; the 
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Figure 2.21 Supervising scheme of the acoustical signal in cutting. 


natural nominal frequency of the transducer of acoustical emission can 
be approximately 500 kHz. This scheme has provided an identification 
method for acoustical emission sources specific to the cutting process; in 
addition, this scheme allows the determination of the characteristic fre- 
quencies and the realization of a correlation between the development 
of the cutting process and magnitudes characteristic of the acoustical 
emission signal. 

The main sources of acoustical emission in the cutting process are 
(according to Souquet, Weber, and Dornfeld; Fig. 2.22): 

Plastic deformation in the shearing zone 

Friction of the chip on the rack face of the cutting tool 

Cutting edge deposits 

Friction between tool and processing part 

Collision and friction between chip and part 

Fragmentation of the chip 

Friction of the chip on the chip breaker [70, 84, 107] 


5 



Figure 2.22 Main sources of acoustical emission in the cutting process. 


The analysis of the signals of these sources provided data for a math- 
ematical model that connected the effective value of the AE signal 
(AE#ms) with the parameters of the cutting regime (cutting speed, 
cutting depth). Experimental research on turning and boring has con- 
firmed the validity of the model for orthogonal cutting. In the milling 
case, difficulties occur because of the presence of shocks when each tooth 
touches the part. 

The proposed mathematical model has proved efficient both for su- 
pervising the evolution of the cutting tool wear and to detect its break- 
age, especially where the cutting conditions do not allow visualization 
of the tool (boring): 


AE rms — C sin a 


. T r , t\ cos a 1 ^ . sin 6 

T kbiV( — — — — — + -(l + 2/ z ) , 

sm<pcos(0a) 6 cos {(pa) 



where C is an experimental constant, a is the tool's attack angle, <p is 
the cutting angle, is the cutting stress, 6/ is the cutting length, li is 
the chip-tool contact length, ti is the thickness of the uncut chip, and 
V is the tangential cutting speed. 

The acoustical emission signal (AE rms), converted in the fre- 
quency domain, can look different depending upon the emplacement 
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Figure 2.23 Acoustical emission signal 



rms ), converted in the frequency domain. 


of the transducer, as shown in Figure 2.23. The most advantageous 
positions for emplacement are directly on the part (1) or on the part 
clamping device (3,2). 

Supervising the acoustical emission has imposed the grouping of 
the influence factors into three categories: parameters bound to the 
nature of the processed material (hardness, the capacity of being pro- 
cessed); parameters connected to the material and geometry of the cut- 
ting tool; and parameters of the working regime (the speed and depth 
of cut and feed). Concerning the factors from the first category, it has 
been observed that the level of acoustical emission varies proportionally 
with the material hardness, and the differences of the capacity of being 
processed modify only slightly the characteristics of the AE signal. For 
example, for an increase of hardness with 300 Brinell units, the level 
of acoustical emission can increase by approximately 15 dB. As to the 
factors in the second category, the estimation is that not the nature of 
the tool material but the way that materials are preponderantly degrad- 
ing (abrasion, fissure, chemical diffusion) affects the acoustical emission 
signal. The modification of the cutting tool geometry has a significant 
effect on the energy of the acoustical emission signal. From this point 
of view, the most important parameter is the rack angle; experimental 
research (Kannatey-Asibu and Dornfeld) has shown that the medium of 
the acoustical emission signal energy sensibly modifies on the increase of 
the rack angle (Fig. 2.24). One of the effects of the rack angle increase 



Figure 2.24 Medium of the acoustical emission signal energy depending on 
the rack angle. 


is the increase of the contact length chip-tool; consequently, a similar 
dependency will be observed also between this rack angle and the energy 
of the acoustical emission signal. Systematic tests (Moriwaki, Souquet, 
Deschamps, Weber, Herberger, Diei, and Dornfeld) concerning the in- 
fluence of the cutting regime parameters on the acoustical emission level 
have shown that the effective value of the acoustical emission signal is 
proportional to the square of the cutting speed (Fig. 2.25a) and also to 
the square of the depth of cut (Fig. 2.25b). The influence of the feed 
magnitude on the AE signal is negligible [32, 33]. During experimental 
research, it has been noticed that the use of cooling-lubricating liquids 
decreases the effective value (RMS) of the acoustical emission signal. 




Figure 2.25 Influence of cutting regime parameters on acoustical emission 
level; acoustical emission signal is proportional to: (a) the square of the cutting 
speed; (b) the square of the depth of cut. 



2.8 TECHNICAL DIAGNOSIS IN FLEXIBLE 
PROCESSING SYSTEMS 

Modern machine tools are integrated in flexible production systems that 
include parts and/or tool feeding systems, parts and/or tool store sys- 
tems, numerical or computer- assisted control systems, systems of adap- 
tive or optimal control, and so on, and each component of the system 
develops specific actions. From this point of view, a matrix representa- 
tion of the production system (Fig. 2.26) leads to a better understanding 
of its problems; one of these problems may be, for example, the predic- 
tion and diagnosis of the machine tool and/or cutting tools destined to 
fail [145, 35]. 

The normal functioning state of the flexible cutting systems sup- 
poses the optimal and simultaneous accomplishment of many fundamen- 
tal functions: 

Automatic processing 

Automation of the flows (materials, tools, parts) 

Control, supervisory, diagnostic, and automatic maintenance 

Informational integration 
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Figure 2.26 Matrix representation of a production system. 
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Figure 2.27 Expert system for the direction of flexible manufacturing systems. 


The laboratory of Machine Tools of the University of Aachen 
has devised an expert system for the direction of flexible manufactur- 
ing systems [163, 83]. In its frame, the diagnostic function appears on 
the most complex level, and the monitoring has a permanent character 
(Fig. 2.27). 

Accomplishing the diagnostic function involves a rapid and correct 
evaluation of the damage messages and deduction of the reasons for these 
damages, possible or existent. In this respect, an important amount of 
information is processed, considering that damage may also have many 
different causes. A correct diagnosis most often requires an oriented re- 
quest of supplementary information; statistical and heuristic data are 
requested and the system must have the capacity to adapt to the new 
experiments and accumulations. This is not possible using conventional 
programming methods because the combinations (message of damage, 
the cause) are practically infinite; that is why the expert system does 
not describe the problems as sequential algorithms, but as relations and 
functional dependencies. 

The presented expert system also has the capacity to recommend 
and even to execute a series of corrective actions as a result of estab- 
lishing the technical diagnosis. The reaction time is very short, the in- 
formation processing usually being done in real-time. 


2.9 CONCLUSIONS: THEORETICAL AND 
EXPERIMENTAL RESEARCH 
PERSPECTIVES 

Supervision of the functional state of machine tools and of the man- 
ufacturing processes developed by them have become subjects of wide 
interest in the last decade. The most important reasons for this phe- 
nomenon are: 

Unexpected stops of machines produce significant financial loss to 
companies, especially in those plants that use optimized plans 
for their production. 

Monitoring and diagnosis lead to the increase of the machines’ 
autonomy, allowing their intensive use in two or three shifts. 

Monitoring and diagnosis lead to increases in productivity and pro- 
duction quality. 


The cost of installation and maintenance of the monitoring and 
diagnostic systems is decreasing continuously, and at the same 
time the cost of the machine tools is increasing because of the 
increase of their complexity. 

By analysis of some technical parameters that correlate with the 
functional state of the machine tool (mechanical vibrations, noise and 
acoustical emission, temperature, force, cutting moment, etc.), a series of 
monitoring and diagnosis methods has been identified. Two classification 
systems of diagnosis exist: diagnosis methods of surface, which indicate 
only the functioning state and/or the existence of damage; and diagnosis 
methods of profoundness, which allow for the appreciation of the nature, 
location of damage, and time until the machine stops. 

Not all acquisitioned signals from the functioning machine tool 
contain useful diagnostic information. An example of this is the forced 
vibrations coming from a discontinuous cutting process, vibrations that 
provoke fluctuations of the cutting forces, and the mechanical resonance 
of the machine. Another example is the trepidation phenomenon (chat- 
ter) that induces vibrations in the entire machine level. Both examples 
present structural vibrations but they are regarded as major damage by 
an inadequate monitoring and diagnosis system. 

The literature has shown that relevant research exists concerning 
diagnosis of some elements of machine tool structure which are common 
to other types of machines and equipment (bearings, gearing, and belt 
transmissions). The diagnosis of cutting tools and the cutting process 
is an intensively investigated domain. The methods used in diagnosis 
are very diverse, and dependent upon the university or researcher who 
investigates them; often they do not have a systemic character or can 
not be reproduced. Considering these aspects, this research endorses the 
systemic approach of the diagnostic methods useful in the machine tool 
domain and their application on the kinematic chain level, especially 
for the feed kinematic chain, where the diagnostic information is miss- 
ing. This work proposes the theoretical and experimental verification of 
the physical model of the investigated feed kinematic chain, in order to 
remove the possibilities of false diagnostics provoked by an inadequate 
structural behavior. 

Theoretical and experimental research should be oriented to the 
determination by calculation of the characteristic frequencies of the func- 
tioning of the specific feed kinematic chains (bearings, ball screws, roll 
cam follower, lug) because of their role in identifying mechanical damage. 



Cognizant of the recent performance of computing systems, this 
work proposes the introduction of virtual instrumentation in the vibroa- 
coustic diagnosis of machine tools, by elaborating and providing users 
with virtual diagnosis apparatus that has increased performance capac- 
ity compared with traditional investigation equipment. 

Last but not least, the role of this work is to create a databank 
based on the findings of technical diagnosis by vibroacoustic methods 
so that further developments of machine tool diagnosis on the level of 
expert systems can be realized. 



3 


Systemization and Application of a 
Vibroacoustic Methods Package 
to Establish Machine Tool 
Technical Diagnostics 


3 . 1 THEORETICAL METHODS TO 
ESTABLISH TECHNICAL DIAGNOSTICS 

By clearly establishing the relations between the level of the supervised 
parameter or parameters and the functioning state of the technical sys- 
tem, pertinent methods of study have been obtained. 

Theoretical research uses statistical and / or probabilistic evaluation 
methods, mono- or multiparametric. The purpose of these methods is 
either to establish a series of decision algorithms on the basis of analysis 
of the possible state of the technical system, or to evaluate the influence 
of not mentioning the state of the technical system on the global cost of 
operation. 

3.1.1 Statistical Methods 

Statistical methods start from the evaluation of the minimum risk of 
spoilage. A limit level xo is adopted for the supervised parameter x and 
the repartitions of this parameter in both its states are supposed to be 
known: fi(x) for the normal state _D 1? and / 2 (:r) for the wrong state D 2 , 
as presented in Figure 3.1. It can be observed that adopting the limit 



level implies a certain level of risk in decision making, illustrated by 
the hatched zones from the D\ and D 2 domains, under the repartition 
functions. Two types of risk occurred: a, the risk of false alarm or the 
producer’s risk, and /3, the risk of not accomplishing the goal or the 
beneficiary’s risk. 

The Neyman-Pherson method is a statistical method that can be 
applied to technical systems with a single characteristic parameter. The 
method consists of establishing an a level, the maximum level allowed 
for the false alarm probability: 

pOO 

/ fi(x)dx = a 
J Xq 

From this a level, the risk value of the x 0 parameter is deduced. The 
same value may also result from the maximum b allowed level ( b < 0.05) 
in the case of an unaccomplished purpose: 




f 2 {x)dx = b 



When establishing the two levels (a and 5), the number of technical 
systems drawn out from operation must be higher than the number of 
those technical systems for which a spoilage is to be expected as a result 
of the inevitable errors when specifying the functioning state. 

If the probabilities that the system is in its normal state P(D\) 
or in its incorrect state P(D 2 ) are established and C 2 \ equals the 
cost implied by a false alarm, C\ 2 equals the cost implied by nonac- 
complishment of purpose, and Co equals the cost of nondetermination 


of the functioning state, the cost of the global risk is given by the 
relation: 


OO 


X 


a 


c = C 2 iP(D 1 ) I fi(x)dx + C 12 P(D 2 ) I f2(x)dx 

Xb 


— OO 


+ Co 


Xb rXb 

P(-Di) / fi(x)dx + P(P> 2 ) / f2(x)dx 

X a J X a 


( 3 . 3 ) 


Deriving the expression function of x a and x^ values of the x pa- 
rameter, the minimizing conditions of the global cost are determined as 


fl(Xa) _ ( C 12 -Co)P{D 2 ) 
f 2 (x a ) C 0 P(Pl) 

fi(x b ) _ C 0 P(D 2 ) 

f2(x b ) (C21 - Cq)P(Di) 



This method is used in cases where the costs of nonrealization of the goal 
or a false alarm are high. In this case, the existence of the nondetermined 
zone situated between the x a and x\, values of the supervised parameter 
x is admitted. 


3.1.2 Probabilistic Methods 

Probabilistic methods allow the establishment of each state character- 
ized by a number of parameters with discrete or continuum repartitions, 
using a diagram of a number of states of the technical system. On the 
basis of this correlation of state parameters, decision rules or parameters 
can be established. 

The Bayes method is a multiparametric probabilistic method that 
allows determination of the most significant diagnostic, by evaluating 
the state combinations of the significant parameters. If in order to de- 
termine the diagnostic, the functioning state is noted Th, Xj is a simple 
parameter, and if P(Di/xj) is the probability of the Di diagnostic in the 
known conditions of the influence of the Xj parameter on the system, 
and P(xj/Di) is the probability that the Xj parameter manifests in the 
Di state, then the probability that both events manifest (D^, Xj) is given 
by the following equation. 




P(D i ,x j ) = P(D i )P 


P(xj)P 


( 3 . 5 ) 



At a certain moment the diagnosed system can be in a unique state so 
the following relation is true. 

n 

E p (^) = 1 o-o 

i= 1 


If a concrete realization X is considered in the group of aq, aq, • • • > Xj 
parameters, then the Bayes relation becomes: 



It is considered that the Xj parameters are usually independent, which 
leads to the simplification: 


P(X/Di) = P(X 1 /D l )P(X 2 /D l ) . . . P(X n /Di ) (3.8) 


This simplification is used many times in practice, despite some in- 
evitable interdependencies. 

In order to simplify Bayesian analysis, diagnostic matrices can be 
made. In these matrices are inscribed the following probabilities deter- 
mined by statistical research: P(D/x), the probability of the x parameter 
manifesting in D state and P(D), the probability of finding D state. The 
most significant factor is recognized in matrix making P(D/X) = max. 

Considering the following equation, 



the basis of the estimation calculus of the technical diagnosis must be 
reevaluated if the probabilities P(D/x) are under (0.4 . . .0.5). 

Observations: 


1. The complex of X parameters, which were theoretically or ex- 
perimentally determined, may define a D{ state if P(Di/X) > 
Pi , where Pi is the limit probability for the Di state (Pi ~ 0.9 
is recommended). 

2. The analysis' volume increases exponentially with the number 
of functioning states and parameters considered significant for 
these states. 



3. The application of some rigorously controlled selection criteria 
leads to reduction of the number of significant values of X 
parameters. 


3.2 EXPERIMENTAL METHODS TO 

ESTABLISH TECHNICAL DIAGNOSTICS 

The functioning state of a machine tool may be correlated with divers pa- 
rameters such as: (a) vibrations, (b) noise and acoustic emission, (c) tem- 
perature, and (d) cutting force and moment. Diagnostic methods based 
on the use of the first two parameters are presently the most widespread. 

The acoustic or vibrator signal, captured by adequate transducers 
and used for diagnostic purposes, may be processed in many ways. The 
directly recorded vibration or acoustic signal is very rarely used in its 
primary form; usually the important elements used for analysis are the 
effective value or the power spectrum, thus the signal’s energy. The ad- 
vantage of methods based on the energy of the vibratory signal is that 
it allows the use of a common inexpensive instrumentation. The inex- 
pensive quality recommends these methods for research laboratories and 
production technical systems. More elaborate procedures use the phase 
of the captured signal. 

Use of the acoustic signal in diagnosis aims to create an ambient 
system that allows a correct acoustic measurement of the sound emitted 
by a machine tool, which is why measurements are done in an anechoic or 
reverberant room, well isolated from environmental sound. These con- 
ditions are hard to realize in a plant environment, but they may be 
produced in research laboratories. 

Another major disadvantage of using the acoustic signal is its inca- 
pacity to identify the diverse noise sources of the machine, and because 
of this disadvantage, if the machine does not pass the noise test, other 
diagnostic procedures will be used in order to determine the source and 
the time until repairing. This disadvantage may be used to correlate the 
vibration tests with the noise tests. It is possible to locate an accelerom- 
eter in the immediate vicinity of a noise source in a variety of locations 
on a machine, to determine if this source is a problem for the machine. 
Precise information necessary to diminish the noise level is also obtained. 

Use of the vibration signal for diagnostic purposes has proved to 
be very advantageous. The isolation of the machine from the vibra- 
tions of the plant environment is simpler and more efficient: a rub- 
ber or polyurethane foam carpet can serve as an excellent isolator as 



the environmental vibrations are of low frequency. From this point of 
view, it must be noted that machine tools need foundations that en- 
sure anti vibration isolation or at least use dampening pads for lean- 
ing. Environment-radiated noise in machine tools produces structural 
vibrations of incomparably smaller amplitude than those of the internal 
mechanisms; these vibrations can be ignored. The main advantages of 
the vibration measurements are twofold: the possibility of fault detec- 
tion and the possibility of locating the source of these faults, that make 
vibration the main signal for diagnostic study. 

There are major differences among diagnosis methods depending 
on their capacity to reach the imposed goal, identify the fault, and es- 
tablish the technical diagnostic. In this respect there are methods that 
indicate the functioning state and/or the existence of a fault, called basic 
diagnosis methods, and methods that estimate the fault kind, location, 
and time until spoilage, called profoundness diagnosis methods. A com- 
parative study of the two types of methods highlights the fact that the 
same physical effects and, many times, the same vibroacoustic signals, 
but treated with different mathematical algorithms, are differently val- 
orized depending on the desired precision degree of the diagnosis. 


3.2.1 Specific Parameters Common to 

Vibroacoustic Methods of Diagnosis 

In order to process the vibration signal in the time domain (Fig. 3.2), 
some parameters specific to the analysis of the vibroacoustic signals are 
introduced: 


The medium value (arithmetic) of the signal: 





(3.10) 


The effective value (effective of the root mean square, RMS): 



The peak value of the signal: 


(3.11) 


positive peak: x v + 
negative peak: x v - 


max .x(t) 
min .x(t) 



Figure 3.2 Vibration signal in time domain. 


The signal’s dispersion (cr-standard deviation of the same signal): 

i r T i N 

= I™ ™ / (x(t) - x ) 2 dt = — — - - ®) 2 
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(3.12) 

The function of amplitude distribution [the probability that the 
vibration amplitude is inferior to a given value x; see also 
Fig. 3.3 (left)]: 
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(3.13) 



Figure 3.3 Amplitude of vibration: (left) inferior to a given value; (right) 
instant amplitude is in a given interval. 





The function density of probability of the amplitude [the probabil- 
ity that the instant amplitude of the vibrator signal will be in a 
given interval; Fig. 3.3 (right)]: 



P(x) — P(x + Ax) 
A ; r — >-0 Ax 


(3-14) 


It can be deduced from the last definition that between the function 
of amplitude distribution and the function density of probability of the 
amplitude exists the relation: 



P(x) 

dx 


or P(x) 



p(x)dx 


(3.15) 


The autocorrelation function (estimation if the vibrator signal re- 
mains similar to itself): 

1 1 

C x{t)=— x(t)x(t + r)dt = — + X k ) (3.16) 

i.l 

Another parameter is added to those above in order to estimate the 
signal in the frequency domain: 

The spectral density of power (the amplitude density in the power 
spectrum) : 


2 

(3.17) 

This function is in fact the Fourier transform of the autocorrelation 
function from the time domain. With this function’s help the total power 
of the signal is obtained by integrating the partial powers of the spectral 
components: 
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(3.18) 


3.2.2 Diagnostic Surface Methods 
3.2.2. 1 Peak Factor Method 


This method evaluates the vibroacoustic signal in the time basis by 
recording the peak value and the effective value, followed by the calculus 
of the peak factor: 

A 

F v = — (3.19) 

X e f 


On the occurrence and development of a fault in a bearing, the 
shocks generated while rolling over the fault make the peak value increase 
substantially, but the shocks have a very small influence on the effective 
value of the vibrator signal. In consequence the value of the peak factor 
increases and this tendency has to be supervised. 

In the technical literature, on the basis of numerous tests, the fol- 
lowing characteristic values are given for the peak factor. 

F v < 10: good bearing; 

F v = 10, . . . , 20: fault conditions occur; 

F v = 20, . . . , 25: incipient fault; 

F v > 25: fault bearing. 

The effective value of the vibrator signal increases even as the shock's 
amplitude from individual faults remains constant, as the bearing is 
deteriorating and more and more faults occur. This phenomenon makes 
the peak factor “fall” to the initial value, toward the end of the bearing 
life cycle, making a deteriorated bearing appear in good condition. Such 
behavior will trick an uninformed user of this diagnostic technique. 

The peak factor method is simple and easy to use. It is especially 
useful in the case of monitoring a large number of measuring points when 
an early warning is not required, and the consequences of spoilage are 
not too great. In particular situations, the method can be completed 
by another diagnosis method in profoundness, for example, Cepstrum 
analysis. 

3. 2. 2. 2 Diagnostic Index Method 

This method has its basis in the use of normalized values of the effective 
and peak parameters of the vibrator signal. 



Normalization aims to report the value of parameters characteristic 
of the analyzed signal, measured at a certain moment £, to a reference 
value for that parameter. For diagnosis, it is useful for the reference value 
of the normalized parameter to be the value measured in the perfect 
operating state of the supervised system. 

Use of the normalization method allows the estimation of the func- 
tioning state of simple mechanical systems (e.g., bearings) and only by 
use of the normalized effective values or peak values of the vibration 
signal. For example, the fatigue of the radial ball or roller bearings 
can be correlated with the effective value of the signal acceleration as 
follows. 

x e f = 20 [mm/s 2 ]: normal function; 

x e f = 20, ... ,45 [mm/s 2 ]: slight traces of light settling; 

x e f = 45, . . . , 90 [mm/s 2 ]: strong settling until visible faults; 

x e f = 90, . . . , 150 [mm/s 2 ]: major faults; 

x e f > 150 [mm/s 2 ]: out of operation. 

The diagnostic error by this method is approximately 50% (see Ta- 
ble 3.1), and decreases to under 30% by normalization of the effective 
value indicator. 

The diagnostic index can be defined with the relation 

(3. 

so it has a value between 0 and 1. 



= Ze/(0)£(0) 

x ef (t)x(t) 


Table 3.1 Error Margin for the Diagnostic Index Method 



Without 

Error 

With 

Error 

Parameter 

normalization 

(%) 

normalization 

(%) 

Effective value, x e f 

X ef (t) 

<50 

Xef (t) 
X ef( 0) 

<30 

Peak value, x v 

x v (t) 

<50 

Xy ( t ) 

x v (0) 

<35 

Diagnostic index, K(t) 

x ef (t)xx v (t) 

<35 

x ef (0)xx v (0) 

X e f(t)xX v (t) 

<25 





The value of the diagnostic index correlates very well with the 
functioning state of the machine parts containing rolling elements: 

K(t) = 1, . . . , 0, 5: good functioning characteristics; 

K{t) = 0, 5, . . . , 0, 2: appearance of accelerating factors of damag- 
ing phenomena; 

K(t) = 0, 2, . . . , 0, 02: effects of damaging phenomena occur; 

K(t) < 0,02: risk of immediate work stoppage by spoilage. 

In the technical literature this method is quoted with an error range of 
25%, as presented in Table 3.1. 

3. 2. 2. 3 Kurtosis Method 

The Kurtosis method is a probabilistic diagnostic method that uses the 
acceleration of the vibrator signal. The amplitude of this signal may be 
represented by the function density of probability [relation (3.14)]. 

A bearing in good functioning condition is the source of some 
stochastic vibration that respects a normal (Gauss) distribution of the 
amplitude. The injurious processes are the cause of the occurrence of 
some supplemental components that change the initial character of the 
signal and implicitly modify the probability density function of the sig- 
nal. Thus, deviations from the normal (Gauss) distribution are born, and 
they can be quantified using the statistic moments of superior order. 

The Kurtosis factor is a fourth-order moment expressed as 

1 f°° 

( X ~ X ) A P{x)dx 

® J _ oo 

where x is the arithmetic mean of the signal, p(x) is the density of 
probability of the same signal, and a is the standard deviation. The value 
of the Kurtosis index for a Gaussian distribution of a signal is /? 2 = 3, in 
a large band of frequency (2.5 -r- 80 kHz) and with a maximum deviation 
of 8%. The increase of this index indicates the debut (/? 2 = 4 ... 6) and 
the existence (/3 2 > 6) of a mechanical fault, respectively. For higher 
values of the index (/? 2 = 9 ... 10), the machine has to be stopped and 
the faulty spare part changed. 

This method has been applied with good results in supervising 
the rolling bearings. Figure 3.4 presents the basic electronic circuit for 
the calculus of the Kurtosis index starting from the measured vibration 
signal on the box of a bearing. The role and type of each electronic 
block can be observed and also the successive transformations of the 
initial signal. The advantage of this method lies in renouncing the use of 


(3.21) 



Figure 3.4 Basic electronic circuit for calculus of the Kurtosis index. 





some comparative values that make the method useful without previous 
preparation. The Kurtosis method is quoted as one of the most precise 
surface diagnosis methods with an error range of 20 to 25%, according 
to the technical literature. 

3. 2. 2.4 Shock Impulse Methods 

Specific faults occur in the functioning of mechanisms that replace the 
slipping friction with the rolling friction such as some irregularities. The 
irregularities of rolling surfaces and bodies, which are due to mounting 
faults and/or wear damage, disturb the process of uniform rolling and 
lead to the occurrence of mechanical shocks. Thus, a wave of mechan- 
ical pressure occurs (shock impulse), which diffuses in the material as 
a spherical wave, with a speed specific to the bearing’s material. The 
residual arrow dampens very slowly compared to the process of acceler- 
ation by shock (Fig. 3.5). The maximum value of the amplitude of these 
large band sonic waves is dependent on the kinetic energy transferred 
and by the square of the shock speed. At constant rpm, the shock speed 
modifies proportionally with the body mass and the depth to which the 
rolling body penetrates. 

The magnitude of the damages can be established by measuring 
the shock impulses and by counting the sizes of the bearing (internal 
diameter D) and the rpm n. The level of maximal values cIBm and the 
white noise level dBc determine the logarithmic value dB sv (decibel 
shock value). The reference value is the initial value dBj , representing 
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Figure 3.5 Wave of mechanical pressure. 



the level of shock impulses of the new bearing perfectly mounted and 
greased. The reference value can be also established empirically using 
the relation: 

dB : = 20(log n + 0.6 log D - log 2150) (3.22) 

The normal value of shock impulses dB n is taken as the basis in order 
to evaluate the functioning state of the bearing: 

dB^ — dBsv — dBj (3.23) 

and the correlation is as follows. 

dBw = 0 . . . 20: good functioning state; 

dB^ = 20 ... 35: functioning state to the limit, because of the lack 
of grease, some mounting errors, or the appearance of the pitting 
phenomenon; 

dBw = 35 ... 60: critical state, risk of “out of operation.” 

The working principle of the method is illustrated in Fig. 3.6. The shock 
wave, which is captured by an acoustic emission transducer, is introduced 
in a narrow band-pass filter. The filtered signal is transformed by an 
impulse generator in an amplitude impulse proportional to the shock 
speed. 

The SPM diagnosis method has a statistic character and the diag- 
nostic error is under 10%. This method especially highlights faults due 
to pitting, material fatigue, and mounting errors. The use of this method 
has the following advantages. 
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Figure 3.6 Working principle of shock impulse method. 



The results of measurements are not influenced by machine tool or 
other equipment vibrations. 

The level of the shock impulses of a damaged bearing, compared 
to a new bearing, may increase by a factor of 1000 (60dB). 

The results of measurements of neighboring bearings can be sepa- 
rated as a result of dampening the high frequency vibrations. 

Figure 3.7 presents typical shock impulse diagrams drawn by the 
SPM method: (a) a bearing in good functioning shape, (b) impure vis- 




FlGURE 3.7 Typical shock impulse diagrams, drawn by the SPM method: 
(a) bearing in good functioning shape; (b) impure viscous grease; (c) fault 
or insufficient grease; (d) bearing not greased; (e) phenomenon of friction 
synchronized with rpm; (f) rhythmic functioning perturbations. 



cous grease, (c) a mounting fault or insufficient grease, (d) bearing 
not greased, (e) phenomenon of friction synchronized with rpm, and 
(f) rhythmic functioning perturbations, for example, by load shocks or 
pressure shocks. 

3.2.3 Profoundness Diagnosis Methods 
3 . 2 . 3 . 1 Spectrum Comparison Method 

Analysis in the frequency domain is based on obtaining the frequencies 
spectrum using the Fourier transform: 

/ CXD 

dt 

-oo 

(3.24) 

/ OO 

G(/y 2 ^ dt 

-oo 

Spectrum comparison is a frequency method that, on the basis of auto- 
matic comparison of the spectra coming from the same measuring point, 
but at different times, recognition of some faults and also, through ten- 
dency analysis, permits estimation of the time until final damage (out 
of use). 

The method is applied in the following steps. 

1. Obtaining the reference pattern. With the machine in perfect 
functioning state, more frequency spectra are recorded in the domain 
2 Hz ... 20 kHz, depending on the memory of the analyzer. The refer- 
ence spectrum is obtained from the mean of the recorded spectra (repre- 
sented as darker shading in Fig. 3.8a), and is kept as it is. The reference 
pattern results as an evolute of the reference spectrum, by widening the 
frequency peaks (represented by the thick line in Fig. 3.8a). This widen- 
ing of peaks is made to avoid the false signals that can occur on the 
translation of whole spectrum at small rpm fluctuations of the moni- 
tored element. 

2. Definition of the limiting profiles. Two other profiles are defined 
related to the reference pattern: the tolerance pattern (of the admissible 
levels) and the alarm pattern (of maximum levels where the interven- 
tion for readjustment is compulsory). These profiles are obtained either 
adding a tolerance of the amplitude, on short frequency domains to 
the reference pattern, or translating the whole reference pattern to a 
preestablished level. Usually, in order to obtain the tolerance pattern, 







Figure 3.8 Steps of the spectrum comparison method: (a) reference spec- 
trum and pattern spectrum, (b) tolerance pattern and alarm pattern, (c) cur- 
rent spectrum, (d) tendency chart. 


the translation is made vertical to the “2 x reference pattern” (repre- 
sented by the dash-dot line in Fig. 3.8b), and for alarm pattern to “10 x 
reference pattern” (represented by the thin line in Fig. 3.8b). 

3. Comparison of spectra. The current spectrum (black in 
Fig. 3.8c), captured in the respective measuring point after a time 
interval or with a rhythm on the user latitude, is compared by super- 
imposing on the tolerance pattern. After comparison, the frequencies 
that have surpassed this pattern and the value of the overages are high- 
lighted in a chart. It is necessary to stop the machine when the alarm 
pattern has been surpassed, in order to find and fix the damage. 

4. Drawing the tendency diagrams. Evolution of the frequencies 
that surpass the tolerance pattern can be analyzed on the basis of some 
successive comparisons for predetermined intervals. A diagram of the 
increase of the amplitude difference between the current spectrum and 


the tolerance spectrum can thus be drawn by points, a diagram which 
is usually linear (Fig. 3.8d). The point whose abscissa indicates the mo- 
ment of intervention (ti) for reestablishment is at the intersection of 
this line with the horizontal line representing the maximum admissible 
value for the analyzed frequency (value taken from the alarm pattern). 
The difference between this moment and the time when the last current 
spectrum was drawn represents the interval for which the machine can 
still be used in limited operation conditions. 

3. 2. 3. 2 Evolute Method 

The characteristics of part vibrations change while a fault begins to 
develop in a part that has rolling elements. Each time a rolling ele- 
ment meets a discontinuity on the rolling path, an impulse occurs. This 
impulse periodically repeats with a rate determined by the location of 
discontinuity, the geometry of the part, and the rpm of the driving el- 
ement. These repetition frequencies are characteristic frequencies; they 
are not always easy to find in a standard frequency spectrum because 
they are mixed with vibratory components of a much higher level. 

The impulses resulting from the concussion of the rolling elements 
with the discontinuities appear as an increase of the wide band in the 
domain of the superior frequencies of the vibratory signal. In conse- 
quence, separating this frequency domain using the Zoom function and 
the Fourier analyzers, or using a band-pass filter, a signal that contains 
these impulses is obtained. This signal is leveled in order to find the 
evolute of impulses in the time domain; then the signal is passed in the 
frequency domain and the evolute spectrum results. This spectrum has 
maximums at the impact frequency, thus the nature of the fault and 
also the element that has the fault can be recognized. This information 
is usually the key to the time approximation until stoppage due to dam- 
age. The absolute value of the peaks from the spectrum depends a great 
deal on filtration and on the excited resonance frequencies, and does not 
present interest for diagnosis. 

The maximums indicated by the occurrence and the evolution of 
the mentioned faults are not usually perceptible in the raw spectrum of 
the vibrator signal because their modular frequency manifests as lateral 
bands around the natural frequencies. 

If the periodic modulation function is also a harmonic oscillation 
then it can be written: 


f P (t) = A T + A m cos u M t, 


(3.25) 


where At is the amplitude of the carrier frequency, Am is the amplitude 
of the frequency of modulation, and ujm is the modulating pulsation. 
The following relation is deducted from the previous relation. 


x{t) = A t cose u T t 




COS {lot + c UM)t + COs(c dT ~ ^m)^\ 


(3.26) 


where lot is the pulsation of the carrying signal. 

The lateral bands that appear are visible in the real signal spec- 
tra only in very rare cases. Most frequently, they are superposed with 
neighboring components of the frequency. The classical models for sig- 
nal analysis can not draw sure conclusions from the analysis of such a 
spectrum, and that is why the extraction of the e volute curve is used. In 
order to extract this evolute curve, a modality to “posit ivate” the signal 
in the time domain must be found; this may be done using one of these 
methods: 


Forming the mean value of short time 

Forming the effective value of short time 

Directing and filtering of the signal 

Using the Hilbert transform 

Passing directly to Cepstrum analysis (Fig. 3.9) 

The mean value of short time can be calculated as 

t T 

= x*(t)dt 


where 


x*(t) = 0 if x(t) < 0 
x*(t) = x(t) if x(t) > 0 


or 
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where 


if X{ < 0 
if Xi > 0 


(3.28) 
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Figure 3.9 Diagram of amplitude difference increase between current spec- 
trum and tolerance spectrum. 






The effective value of short time can be calculated as 



(3.29) 
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(3.30) 
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The length of the time interval has a special importance; usually the 
following relation must be respected, 


ft > T > /m 


(3.31) 


where fc is the carrying frequency and Jm is the modulation frequency. 
It can be mentioned in conclusion that the time interval T must be at 
least as large as the period of the carrying vibration having the smallest 
frequency. 

The most exact method of extraction of the evolute is the applica- 
tion of the Hilbert transform of the time domain signal: 


H{x{t)} 



(3.32) 


When passing in the frequency domain applying the Fourier transform, 
the following relation is obtained, 


T{H{x{t)}} = sign (/)) 


(3.33) 


Thus, using Fourier analyzers, changing the phase in the frequency do- 
main is simple to realize; this change corresponds to a 90° rotation in 
the complex space. 

If we consider the original signal and the Hilbert transform of this 
signal as real and imaginary parts of a complex function xu\ then this 
is the evolute of the time signal: 

x H (t) = x(t) + jH{x(t)} (3.34) 


If the frequency analyzer does not possess the algorithm needed to obtain 
the Hilbert transform, the method of signal directioning and filtering 
can provide very good results. The signal captured by the accelerometer 
contains, as specified above, lateral bands of the fault frequency and 



of the rotation frequency around each characteristic frequency. In other 
words, the carrying signal containing the whole characteristic frequency 
is modulated by a signal containing the fault frequency fa, the rotation 
frequency / r , and their snbharmonics. 

The reduction of the frequency content of the signal is realized 
when passing through a band-pass filter centered on one of the char- 
acteristic frequencies. The bandwidth of this filter is taken in such a 
way that the lateral frequencies induced by the fault can also pass. To 
make this possible, the largest value of the fault frequency / m ^, and the 
bandwidth Z& > ±/ m ^ (usually ±2 fmd) must be approximated. Conse- 
quently, when leaving the filter, the signal will contain the characteristic 
frequency and its lateral frequencies fd(2fd) si f r (2f r ). The signal is 
then demodulated and passed through a low-pass filter to eliminate the 
eventual lateral frequencies that are around the first superior harmonic 
of the central characteristic frequency. 

Processed in the manner outlined above, the signal is then passed in 
the frequency domain by a Fourier analyzer, and the evolute spectrum 
is obtained, a spectrum that contains all the mentioned frequencies, 
including the frequency specific to the damage. 

Separation of fault damage is realized using a “normalization 
method” similar to that presented in the diagnostic index method. This 
separation will be possible if previously a spectrum of evolute in a perfect 
functioning state of the machine was obtained. Normalization consists 
of the distribution of amplitude values of the corresponding frequen- 
cies, frequency with frequency, between the spectrum with fault and the 
spectrum without fault, using a relation such as 

ma x{A d (i)} 

where A(i) is the amplitude of the frequency i, and the indices d and b 
have the meaning “fault” and “good.” 

3 . 2 . 3 . 3 Cepstrum Analysis Method 

Cepstrum was first proposed in 1963, and it was defined then as the 
“power spectrum of the power spectrum logarithm.” The reason for this 
definition is not entirely clear because even in the original work it is 
compared with the autocorrelation function that can not be obtained 
as an inverse of the Fourier transform of the power spectrum. Later, 
Cepstrum was defined as the “inverse of the Fourier transform of the 





logarithmic power spectrum,” clarifying the link with the autocorrela- 
tion function. At the same time, a function similar to Cepstrum has 
been defined as the “inverse of the Fourier transform of the complex 
logarithm of the complex spectrum,” and, in order to distinguish the 
Cepstrum defined above, it was named complex Cepstrum, and the first 
was renamed power Cepstrum. 

The cepstru denomination is paraphrased from the word spectrum; 
in the same mode a series of other terms is introduced: quefrency from 
frequency, ramonic from harmonic, gamnitude from magnitude (ampli- 
tude), lifter from filter, and others. 

The Theoretical Basis of Cepstrum Analysis. Using the T{ } sym- 
bol to indicate the Fourier transform of the quantity in parentheses, the 


definition used for the power Cepstrum is 

Cp( T ) = F~ 1 {LogF xx (f)} (3.36) 

and for the complex Cepstrum is 

Cc(r) =T~ y {log F x (f)} (3.37) 

where F xx (f) is the power spectrum of the time signal f(t ), 

Fxx(f) = \r{f x (t )}\ 2 (3.38) 

and F x (f) is the complex Cepstrum of f x (t) 

F x (f ) = = a x (f) + ib x (f ) = (3.39) 

expressed in real and imaginary components, respectively, in amplitude 
and phases, with f x (t) real. From the previous relation, the complex 
logarithm of F x (f) results in 

log F x (f) = log A x (f) + (3.40) 


The Cepstrum analysis applications imply the detection of periodic struc- 
tures from the spectrum (harmonics, lateral bands, echoes, reflections), 
followed by the separation of source from effects of the transmission path 
of the signal. The fundamental characteristic of this procedure is the rep- 
resentation, for any physical system, of output signal y(t) as a convolute 
of input signal x(t) with the system response in frequency function h(t): 


y(t) = x(t) * h{t) 


(3.41) 



Applying the convolution theorem, this relation transforms itself in a 
multiplication in the frequency domain, and applying the logarithm the 
relation becomes additive: 

Y(f) = X{f) xH(f) 

log Y(f) = log XU) + log H(f) (3.42) 

This additive relation maintains itself also in cepstru because of the 
linearity of the Fourier transform; so here, the effects of the source and 
transmission are additive: 

F- 1 {logY}=T~ 1 {logX}+T- 1 {\ogH} (3.43) 

The power Cepstrum and the complex Cepstrum offer the possi- 
bility of separating the components because of the excitation from the 
components due to the dynamic characteristics of the structure (the way 
covered by the signal), because of the logarithm properties. In addition, 
the power Cepstrum and the complex Cepstrum offer the advantages 
of the possibility of eliminating the last-mentioned components and re- 
constituting on this basis the excitation. At the same time, Cepstrum is 
less sensible to these dynamic characteristics of structure: for the same 
excitation, the spectrum of frequencies depends more on the point where 
the transducer is mounted on a structure, as Cepstrum is less influenced 
by the choice of the measuring point. 

Another important advantage of Cepstrum analysis is the fact that 
all the lateral bands of a certain frequency are mainly grouped in a single 
line (ramonica), this line containing significant information concerning 
the medium height of the lateral band. The use of Cepstrum analysis 
for machine tool diagnosis is based on the analysis ability to detect the 
periodicities from spectra (families of harmonics and lateral bands), and 
also on the insensibility of the analysis related to the path from the inter- 
nal source to external measuring point. These advantages have imposed 
the use of Cepstrum analysis for diagnosis of mechanical systems, es- 
pecially for gearing transmissions (reduction gears, gearboxes, etc.) and 
for bearings. 

3. 2. 3. 4 Acoustic Emission Method 

Acoustic emission is defined as the succession of elastic waves generated 
by the liberation of the internal energy stored in a structure. The acous- 
tic emission manifests in the high frequency domain (/ > 100 kHz) by 
elastic waves detectable as vibrations on the supervised structure. As 



a working method, acoustic emission represents a nondestructive tech- 
nique, capable of detecting when and where a fissure or crack occurs. The 
nature and the causes of these faults are then investigated by comple- 
mentary methods. Four main sources of acoustical emission exist: struc- 
tural dislocation movements, phase transformation, friction mechanisms 
(microfrictions, microconcussions), and fissure forming and extension. 

In the case of dislocations (i.e. , displacement of an imperfection 
of line in a crystalline network), that occur as an avalanche, the de- 
tected signal is of the continuum type, as for the phase transformations 
(i.e., forming of martensite in carbon steel), the signal is of the impulse 
type and can be detected for each transformed grit. 

The fissures occur in the material nodes and points where the local 
effort surpasses the fracture stress. Thus, new surfaces are formed and 
energy is released, energy that is partially converted in acoustic emission. 
The signal is of impulse type and high frequency. At the same time, the 
friction mechanisms entered in action also emit impulse type acoustic 
signals. 

In acoustic emission the signal amplitudes cover a wide domain; in 
relative units, these amplitudes are: 1 to 10 for structural movements, 
5 to 1000 for phase transformations, and 20 to 1000 for fissures. 

The propagation of acoustical emission is similar to that of radio 
waves. The source emits packages of spherical waves that are immedi- 
ately affected by the numerous intersected surfaces that create reflections 
and surface waves (Fig. 3.10). The nonhomogeneity of the propagation 
environment also distorts the wavefronts. In consequence, the mathe- 
matical relations that describe the real propagation phenomenon are 
very complicated, and they introduce serious difficulties in sources and 
effect localization and measurement. This fact limits the usable domain 
of methods based on acoustical emission for large structures made of 
steel, where a specific incertitude coefficient can be tolerated. 

Many types of treatment of the detected signal are available for 
the acoustic emission study; that is why at all times a general evaluation 
must be done first. 

Counting of impulses is necessary in impulse type signal treatment. 
This method supposes an evaluation of the number of impulses that pass 
over a previously established limit value. In order to do this evaluation, 
the measuring chain contains a discrimination device for the impulse 
amplitude level, followed by an impulse counter (per total or per time 
unit). The simple count of impulses (Fig. 3.11c) can be improved by 
an evaluation of the impulse area (Fig. 3.11a), which also takes into 



Figure 3.10 Propagation of acoustical emission. 
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Figure 3.11 Impulse count: (a) evaluation of impulse area, (b) combination 
of limits, (c) simple count of impulses. 


consideration the duration of the impulse, eventually by introducing a 
combination of limits (Fig. 3.11b). 

The meaning in amplitude of the impulses is realized in cases when 
the acoustical emission manifests itself by signals of the continuum type. 
Calculation of the effective value (RMS) is significant because it is pro- 
portional to the power/energy of the signal. This method is encountered 
sometimes under the name of energetic analysis of the impulse. 

It is useful to locate the source of the information provided con- 
cerning the source’s and the wave trajectory’s characteristics and mod- 
ifications, localization of the primary faults, and the reduction of the 
times of re-putting into function. In order to locate the source, some 
controlled elastic waves are created from ultrasonic sources. The ultra- 
sonic wavetrain is characterized by a factor of the simultaneous wave 
that depends on the frequency of the impulse succession, duration, and 
number of peaks that surpass a certain limit. In this way, breaking zones 
can be located with an anticipation in materials. 

In order to locate a source, three transducers are used. The time 
difference between the moments when the signal arrives at two trans- 
ducers determines a plain hyperbola, if the propagation speed in en- 
vironment is known. The intersection of obtained hyperbolas from the 
pair of transducers (1;2), (1;3), (2;3) defines the correct position of the 
source. 

3.3 RECOMMENDATIONS CONCERNING THE 
USE OF DIAGNOSIS METHODS FOR 
MACHINE TOOLS 

The multitude of noise and vibration sources produced by the complex 
dynamic structure of machine tools provides serious problems for diag- 
nosticians. However, it is expected that noise and vibrations, coming 
from functioning machines, are generated by the characteristics of the 
component mechanisms and subassemblies of the machine tool. At the 
same time, the working regime parameters and the characteristics of the 
processed materials have a great influence on these phenomena. 

In order to apply vibroacoustic diagnosis to machine tools, the 
natural sources of noise and vibrations must be identified, and also the 
parameters that best correlate with the functioning state of the machine 
tool or its subassemblies and elements. It is indicated that the moni- 
toring of the machine tool functioning state or its component elements 
may be realized using working techniques in the time domain, and the 


establishment of the technical diagnostic, using the method of analysis 
in the frequency domain. 

Very good results have been obtained using the Kurtosis method, 
which has the advantage of immediate application, without knowing the 
machine history. In the appendices to this book the use of this method 
of diagnosis of the functioning state for the majority of the components 
of feed kinematic chain physical models is illustrated. 

Concerning the profoundness diagnosis, the best results are ob- 
tained by the method of spectrum comparison that, besides indication 
of nature and place of fault, can also estimate the remaining time until 
a machine fails because of damage. 
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Theoretical and Experimental 
Research: Diagnostic Analysis 
of the Feed Kinematic Chain 
Mechanical System 


4.1 INTRODUCTION TO MECHANICAL 
SYSTEMS ANALYSIS 

The dynamic study of a mechanical system has as its goal the knowledge 
and, sometimes, the prediction of the dynamic behavior, in order to 
understand which way the modifications made to a system or to the 
disturbing sources affect the dynamic behavior of this system. Physical 
or mathematical models must be created in this respect, models that 
can approximate as well as possible the behavior of the real system. 

The creation of a system model usually supposes the use of an 
adequate combination of theoretical and experimental methods whose 
succession is determined by the goals of the research and the character- 
istics of the system. Theoretical analysis of the system is carried out to 
determine the dynamic properties of the system and is based on equa- 
tions that characterize the system, also known as “model construction 
equations.” The theoretical analysis implies certain steps that must be 
followed: 

1. Establishing some simplifying hypothesis for the system in 
order to reduce the analysis effort 




2. Establishing the sum equations for the masses, the energies, 
and/or the impulses that occur in the system 

3. Establishing the phenomenological equations in the case of ir- 
reversible processes (i.e., heat propagation) 

Generally, a system of ordinary differential and/or partial deriva- 
tive equations is obtained, and these equations represent the theoretical 
model of the system of given structure and parameters. This model is fre- 
quently very complex and can not be used as it is. Model simplification 
is necessary and can be accomplished by 

Linearization of partial derivative equations 

Approximation by ordinary differential equations of the partial 
derivative equations 

Reduction of the order of ordinary differential equations 

The resulting theoretical model contains the functional link between the 
physical data of the system and its parameters. The theoretical model 
is recommended only if sufficient elements are known (elements that are 
connected to the laws that characterize the system’s dynamic behav- 
ior), or if the theoretical model’s behavior must be simulated. It must 
be stressed that theoretical analysis enables the researcher to establish 
equations that describe the dynamics of the modeled system even when 
the system is in the state of design, not available for experiments. 

Experimental analysis of a system, a method also known as “identi- 
fication of the system,” proposes to determine the mathematical model 
on the basis of measurements of variables that characterize its evolu- 
tion in a certain regime. In this situation, we always start from knowl- 
edge about the system gathered through previous theoretical or other 
analysis. Then, the input and output of the system are measured and 
evaluated using an adequate method for identification, which is the 
link between the measured variables. It should be noted that the in- 
put magnitudes can be signaled by the normal function of the sys- 
tem or can be artificially introduced signals. The experimental model 
contains numerical values as parameters whose functional connection 
with the physical data remains unknown. This model, which gener- 
ally describes the momentary dynamic behavior of the system, is ob- 
tained with minimal effort, and can be used to lead or predict certain 
variables. 

Theoretical analysis can use the results of the experimental anal- 
ysis to verify the precision of the theoretical model or to determine 
the parameters of the model that cannot be determined otherwise. The 



experimental analysis can use the results of the theoretical analysis es- 
pecially for the model structure. 

The models that result from the two types of analysis can be com- 
pared. Reexecuting certain steps of the analysis can eliminate the even- 
tual nonconcordances. 

In conclusion, it is evident that the correct model of a system is 
necessary to perform an adequate combination of theoretical and exper- 
imental methods. One of the possible combinations of the two types of 
analysis is presented in Figure 4.1. The analysis steps and their succes- 
sion are influenced during the procedure. 

4.2 THE EXPERIMENTAL STAND 

The kinematic feed chain ensures the cyclical positioning on one of the 
generating trajectories (G, D) of the generator element (De or Ge)- The 
classical structure of the kinematic feed chains is presented in Figure 4.2, 
where the notations represent the following. 

OP represents the start /stop of machine movement. 

I is the reverse for the feed movement direction. 

Ks is the mechanism of periodic transmission of the movement 
(if the advance is intermittent). 

M r is the mechanism of feed magnitude adjustment. 

S is the safety mechanism (protection to overloads). 

Mt is the mechanism that transforms the rotation movement into 
translation movement. 

In the most recently proposed structures of the feed kinematic chains 
some of the above-presented mechanisms have disappeared, their roles 
being taken by the driving element which is usually a continuous current 
motor with a large domain for rpm adjustment. 

The problem of establishing the technical diagnostic in the kine- 
matic feed chains of modern machine tools necessitated building an ade- 
quate testing stand, one that contained all the necessary elements. Fig- 
ure 4.3 presents the structural kinematic schema of the studied kinematic 
feed chain. 

The driving element is a direct current motor, type SMUC-35, 
having the characteristics: 

The actual moment, M n = 35 Nm 

The nominal intensity of current, I n = 28 A 

The couple with the ampere coefficient, K R = 13 Nm/A 

The inertial moment, j = 0.065 kgm 2 
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Figure 4.1 Combination of theoretical and experimental analyses. 




The leading screw (which has a step of 10 mm and a double nut 
with balls (D = 60 mm, d = 30 mm, p = 10 mm, z = 103 bile, D w = 
6 mm, and kRb = 12 * 10 4 daN/mm) to compensate the axial play) is 
fitted in two identical bearing cases that contain a pair of axial bearings 
51307-P5 and one pair of radial bearings 6207-P6. 

The slide movement along the body’s gnidings is ensured by the 
cam followers with rolls type GRT 3 (L w = 14 mm, C Q = 10200 daN, 
/ = 1,25 (i m at 1000 daN, k cr = 80 * 10 3 daN/mm), mounted in “O”. 

The material used for slide and body is cast iron (Fmn having 
E = 1,6 * 10 4 daN/mm, G = 4,5 * 10 3 daN/mm, 7 = 7,3 kg/dm 3 , 
a = 18 * 10 -6 grd -1 ) and the leading screw is made from alloyed steel 
(40Crl0 having: E = 2,1 * 10 4 daN/mm, G = 8,1 * 10 3 daN/mm, 7 = 
8,2 kg/dm 3 , a = 0,16 * 10~ 6 grd" 1 ). 

The block schema from Figure 4.4 highlights the main component 
elements of the mechanic stand. 




Figure 4.3 Structural kinematic schema of the studied kinematic feed chain. 
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Figure 4.4 Main constitutive elements of the mechanical stand. 


4.3 THEORETICAL ANALYSIS OF THE 
KINEMATIC FEED CHAIN 
MECHANICAL SYSTEM 

The numerical analysis method proposes to model the physical phenom- 
ena in equations and describes the behavior of the physical systems by 
equations that express laws of mass, impulse, and shape conservation. 

Finite difference methods consist of replacing the differentials from 
the equations of the model with very small finite differences, with the 
restriction of the validity only in certain points of the analyzed domain 
(the points of the discretized network of the model). Because of this 
method the discretized network will have a rectangular shape, thus for 
models with curved surfaces the method will be very difficult to apply. 

Finite element methods (FEM) are based on local approximation, 
and on portions or subdomains that also consider the reports of relative 
dependence between the entire domain and the studied subdomain. The 
behavior of such a cut-up finite element is formulated by the equation: 

[k\ * {d} = {g} (4.1) 


where [k\ is the stiffness matrix of the finite element, {d} is the vector of 
displacements from nodal points, and {q} is the vector of the generalized 
forces that act on the elements by nodes. 



The boundary element method (BEM) supposes that displace- 
ments produced by an elastic concentrate force applied at the origin 
(node) are known. This method is based on the equality between the 
mechanical work which results from the motion of a force system on the 
movement coordinates of a second force system and the mechanical work 
produced by the second system on the movements of the first. By using 
the relations between forces and displacements the tensor components 
that actually represent the problem’s solution can be determined. 

Among these three methods, the finite element method is most pre- 
ferred because of the following qualities: large domain of applications, 
liberty to choose the type of discretization and to use simultaneously 
many more types of finite elements, discretization with variable geom- 
etry, and the possibility of performing calculus on substructures. The 
finite element method represents the “brick” of the discretized physical 
system, the shape, the type, and its equations being the discretization 
principles result. 

4.3.1 The Method and Images 3D Algorithm 

The FEM operates with three types of finite elements that are justified 
by geometry and number of independent coordinates: 

Unidimensional elements (bar type), that have as ends even the 
nodes of the discretization network. 

Bidimensional elements (plate type). The IMAGES 3D program 
includes two elements of this type: membrane (bending stressed 
only) and plate, which also has the thickness of the plate de- 
fined. 

Tridimensional elements (prism, tetrahedral, and hexahedral type). 

The increase in number of elements leads generally to the increase of 
the precision degree. However, over a certain value, the error does not 
decrease, no matter to what degree the number of the discretization 
elements would increase. 

Interpolation functions are used to write the relations concerning 
the deformation and stress state of the finite element. The most used 
interpolation function is the Hermite function. For each type of finite 
element the rigidity matrix is determined (on the basis of the principle 
of virtual mechanical work), in the local axis system. This rigidity ma- 
trix is used to express finite element behavior by the known equation: 
[k\ * {u} = {r}. 



The assembling operation of the rigidity matrices and of the nodal 
force vectors determines the system of equations that characterize the 
structure: [K\ * {U} = {/?}, where the relation’s terms represent the 
result of the aligning of the rigidity matrix and displacement and forces 
from the local system to the global axis system. 

In order to ensure the convergence of solutions to the most exact 
solution, certain conditions must be followed. 

The displacement models should be continuous inside the finite 
element. This condition is accomplished by the use of Gauss 
polynomial interpolation, which represents continuous functions. 

The displacements of the finite element should be compatible with 
the other ones. The displacement compatibility is satisfied only 
if the displacements of the points from every edge depend only 
on the displacements of the nodes that delimit the respective 
edge. 

The displacement models should consider the displacements of the 
entire discretized system and also the constant stress state of 
elements. This supposes avoidance of a very fine discretization 
because when constructing the rigidity matrix, the elements of 
the main diagonal will be of negligible value (zero), which leads 
to a null determinant and thus to a wrong value of the inverse 
of this matrix. 

4.3.2 Discretization of the Physical Model 

The theoretical working model was obtained applying the similitude 
theory for the physical model presented. The reduction scale for the 
length was chosen as A = 2.5 and the other scales have resulted: 

A / — — As — 2.5 

A^ = Xg = A a = A m = 1 (4.2) 

X F = Xf = 6.25 

In other words: using the same material, the linear dimensions of the 
system are multiplied by A, the forces by A 2 , and the elastic constants 
are multiplied by A for translation and by A 3 for rotation, then the am- 
plitudes and the frequencies are multiplied by A, and the efforts remain 
unchanged. 

The following discretization of structure was obtained (Fig. 4.5) 
based on the principles presented in the previous paragraph. The body 



Figure 4.5 Structure discretization. 


was made with a network of 613 points, the base finite element being 
the plate (membrane plates bending) under loads of compression and/or 
bending. The discretization is finer in the bearing zone to better estimate 
the maximum stress points of the structure. 

The body guidings are also plate elements type, but they do not 
have a median longitudinal point that represents the contact point with 
the roll cam follower. This point has been connected to the extreme 
points of the finite element by springs with an increased rigidity com- 
pared to the other elements (k = 10 7 daN/mm). In this manner, the 
guiding-rollers contact is better approximated. Finally, for the longitu- 
dinal body and its guidings, 463 plates with a thickness between 20 and 
40 mm have resulted. 

The longitudinal slide (the table) was discretized using a step of 
15 mm and 120 elements of plate type have resulted. The table’s rolling 
guidings had been “replaced” by spring type elements having the same 
rigidity as the elements that define the link between the body and the 



guiding plates. The leading screw was discretized in only 23 points be- 
cause of its homogeneity and its physicomechanical characteristics which 
are constant in all directions; a bar network that can be loaded to tensile 
and compression stresses results. 

A tubular shaped nut was made using tridimensional finite ele- 
ments (solids), which were connected to the table slide by means of 
plates made from the same material. 

A special problem was the connection between the body and ground, 
which was made by springs connected to the ground, having in count 
the damping and elastic characteristics of the base material. 

4.3.3 Statistical Analysis of Physical Model 

Static loads were applied in order to determine the deformation mode 
of the mechanical system, a system designed as an assembly of elements 
that introduce deformations in a certain proportion for each element, 
and also to determine the zones of the structure having maximum defor- 
mations. In order to find the deformation mode, the load of the model’s 
slide during the cutting process was simulated. The load was made with 
a static force that cumulated the weight of the semiproduct and the 
maximum value of the cutting force. The force was applied punctually 
on the slide decomposed on the following three directions: 1000 daN on 
OX; 2000 daN on OY; and 3000 daN on OZ. The total deformations of 
the constitutive elements were supervised on the above-mentioned three 
directions. 

Studying the displacements on the OX axis, it has been observed 
that the most considerable yield (Ax max = 1.14 fm i) is produced in the 
screw-nut mechanism, as expected. Figure 4.6 presents a cut-through 
model to visualize the deformations to this mechanism. Less impor- 
tant deformations have resulted on the OY direction until A ?/ max = 
—0.447 \i m value (compression) in the guidings zone. 

Hook’s law can be applied to the entire model because the loading 
does not surpass the elastic domain. This law introduces the propor- 
tionality relation between the unitary efforts and the deformations. As 
a result, the tensor of unitary efforts may be calculated and represented 
using the same algorithm. 

As previously mentioned, the deformation state was maximum at 
the screw-nut mechanism, so the nut stresses have been detailed as the el- 
ement with the most important yield. Considering the global stress state 
it can be noticed that a maximum value cr max = 4.51 10~ 2 daN/mm 2 is 



Figure 4.6 Section-through model. 


reached. The analysis continued by decomposing the tensor on the main 
directions (011,0 22,033 components). 011,022 have close values, but 033 
surpasses them by 150 %, being by far most important (the maximum 
tensile strength: 033, max = 3.86 10 ~ 2 daN/mm 2 ). 

Figure 4.7 presents the global stress state of the physical model. 
It is evident that the maximum stresses are in the zone of the move- 
ment transform mechanism and the guidings zone; the body appears 
nonstressed but the table slide is stressed in the rolling guidings zone 
and at corners. 

4.3.4 Modal Analysis of Physical Model 

The method of modal analysis follows to determine certain dynamic 
characteristics by the estimation of the weights matrix, rigidity matrix, 
and damping matrix. As a result of this analysis the natural pulsations 
of the structure, the existent damping in the system, and the type of 
natural vibration modes can be deduced. The first three natural vibra- 
tion modes and their corresponding natural frequencies were highlighted 
for the analyzed physical model using the IMAGES 3 D program. 

The first natural vibration mode (Fig. 4 . 8 ) affects the mechanism 
of movement transformation type screw-nut with balls and corresponds 
to a natural frequency of 259.6 Hz. High deformation of the free end 
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Figure 4.7 Global stress state of physical model 


of the nut is observed during the other subassemblages of the model 
(longitudinal slide, bearings, and body) to remain practically immobile. 

The second natural vibration mode appears at 380 Hz frequency 
(Fig. 4.9), where a sensible deformation of the screw-nut mechanism is 
recorded, a slightly longitudinal and transversal deformation of guidings 
occurs, and a little transversal deformation of the body also occurs. 

The third natural vibration mode (Fig. 4.10) occurs at 491 Hz 
and leaves the screw-nut mechanism almost undeformed. However, a 
powerful deformation of the guidings in a vertical direction and a slight 
deformation of the table slide are recorded. 

The following conclusions can be formulated as a result of theoret- 
ical analysis of the mechanical system of the kinematic feed chain. 

The most important yield from the kinematic feed chain belongs 
to the movement transforming mechanism. 

The nut represents the most stressed element from the point of 
view of tension intensity. 

The bearings rigidity influences the static and modal behavior of 
the leading screw. 




Figure 4.8 First natural vibration mode. 





Figure 4.9 Second natural vibration mode 
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Figure 4.10 Third natural vibration mode. 


The displacements of the slide-guidings-transforming mechanism 
ensemble are of major importance on transversal direction and 
less important on longitudinal direction. 

The theoretical model constructed approximates with very good 
results the real model of the feed kinematic chain of a milling 
and reaming machine AF 180. 

The FEM is susceptible to any factors that contribute to improve- 
ment of static and modal behavior of the analyzed structure. 

The analysis led to the identification of the elements with maxi- 
mum stresses in the structure, in order to monitor them during 
operation. 


4.4 EXPERIMENTAL ANALYSIS OF 
KINEMATIC FEED CHAIN 
MECHANICAL SYSTEM 


The vibration of a mechanical system is caused by external loadings or 
imposed movements that are variable in time, and are generally called 



excitations or perturbations. Movements of different points of the sys- 
tem and the dynamical efforts from its elements are manifestations of 
the mechanical nature of these perturbations, and are called responses. 
The response is conditioned both by the excitation parameters and the 
mechanical characteristics of the system. Solving the vibration prob- 
lems consists of establishing relations among excitation, response, and 
mechanical characteristics of the system. 

The experimental analysis (identification) of the mechanical sys- 
tem proposes a series of functions of excitation and corresponding func- 
tions of response, to facilitate a mathematical description or an analyti- 
cal model of the system. The relations between excitation and response 
are experimentally determined; the most frequently used data are the 
response curves in the frequency of the system obtained by excitation 
with test signals. On the basis of these curves, the identification of nat- 
ural frequencies and vibration modes is made, and also of the specific 
dynamic properties. It is important that during experimentation the in- 
fluence of other disturbing sources is reduced to a minimum. In addition 
the equipment used for structure excitation and response measurement 
should not modify the utilized mechanical system parameters. 

The properties determined under these conditions sometimes differ 
from those obtained during normal functioning, especially in nonlinear 
structure cases. Choosing correct types and levels of excitation leads to 
satisfying results. The dynamic characteristics of a linear system with a 
single input x{t) and a single output y(t), can be described in the time 
domain by the proportion function h(t), and in the frequency domain 
by the response in frequency function that constitute a pair of 

Fourier transforms: 

/ OO 

h(t)e~ iut dt (4.3) 

-OO 

In the case of determinist excitation, the response is given by the integral 
of convolution: 

y(t) = l h(t') x(t — t')d't = x(t) * h(t) (4.4) 

Jo 

In this case the proportion function h(t) represents the response to an 
excitation such as unity impulse (Dirac) x(t) = S(t). The following rela- 
tion gives the function response in frequency: 


H(iuS) = Y (ioS) / X (iuS) 


(4.5) 



where X and Y are the Fourier transforms of excitation and response, 
which can be written as 

/ OO 

x(t)e~ lujt dt 

-OO 

(4.6) 

/ OO 

y(t)e~ lUJt dt 

-OO 

the inferior limit being zero or the real systems. 

If the structure excitation has been made with a force of the har- 
monic type, the x(t) and y(t) functions will have the expression: 

x{t) = x v e iult , y(t) = y v e^ t+ ^ (4.7) 


which makes the response in the frequency function become 


H(icj) — 





From the above relation a very important conclusion results: the module 
of response in the frequency function (\H(ico)\) can be obtained from the 
amplitude-pulsation characteristic (y v /x v — uj); and p is obtained from 
the phase-pulsation characteristic (p — uj). 

These diagrams can be determined experimentally using the si- 
nusoidal excitation of constant amplitude and variable frequency. The 
two characteristics can be drawn either dot by dot, making measure- 
ments at discrete frequencies in the stationary regime, or continuum, 
using a frequency scavenging slow enough to permit the establishment 
of the regime response for each frequency. Using a live frequency ana- 
lyzer simplifies the gathering of responses, a single excitation signal being 
sufficient. The structure of frequency analyzers includes an analyzer for 
functions compartment at whose entrance the signals x(t) = x v sm cut, 
y(t) = y v sin (u;t + p) are applied and a reference signal z(t) = z v cos(c ut). 
This compartment performs the following multiplication, 


x(t)y(t) = — cos p — — cos(2 cut + p) 

Xy 


ry* ry* 


z(t)y(t) = — sin p + — sin(2 cut + p) 

X y Zy Xy Xy 



where the constant terms represent the real and imaginary parts of the 
response in frequency function. Elimination of variable terms is made 



with down-pass filters, making the mean of the products on an integer 
number of cycles of excitation. 

Processing the signals that correspond to the two components, 
the analyzer displays the polar (Nyquist) diagram of the response in 
frequency for the mechanical system (Fig. 4.11). Graphical analysis of 
these diagrams, as proposed by Kennedy and Pancu [82], remains the 
most exact method for determination of the dynamic parameters and of 
the type of natural vibration modes of a complex structure. To do this 
analysis, each buckle of the diagram is approximated with a circle and 
calculus relations established for systems with a single degree of freedom 
are used. Usually the work is done with the hypothesis of a proportional 
damping and, it is considered that, in the vicinity of no matter what 
natural frequency, the contribution of the nonresonator modes is either 
negligible, or constant (independent of the excitatory pulsation). 

Localization of the natural pulsation is made, in this case, using the 
criterion of extreme value of the imaginary component of the diagram. 
Thus, the M point is found (the extreme point on the imaginary axis), 
and, corresponding to it, the cu r frequency. 

To determine the damping, the BC diameter is drawn perpen- 
dicular on O' M\ this diameter intersects the approximation circle in 




Figure 4.11 Polar (Nyquist) diagram. 




Figure 4.12 Experimental installation. 


semipower points B and (7, at uj r ' and uj r n pulsations. The approxima- 
tion factor is calculated with the equations: 



+ LJ ‘r' 


(4.10) 


The modal masses are calculated with the equation: 



1 


(4.11) 


The experimental analysis of the physical model of the kinematic feed 
chain followed the stages highlighted previously. The experimental in- 
stallation used (Fig. 4.12) contained the impact hammer B & K 8202 
and the real-time frequency analyzer B & K 2034. The model excita- 




tion was made striking the longitudinal table in the vertical direction, 
and the structure response was captured from the body by a B & K 
4391 accelerometer. The selected frequency domain was between 0 and 
800 Hz, linear, considering the results of the theoretical analysis. Fig- 
ure 4.13 presents graphics of the magnitude and phase of the response 
in frequency; it must be noticed that the maximum values are located 
at 265, 374, and 474 Hz in the magnitude/frequency graphic. 
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Figure 4.13 Graphics of magnitude and phase of the frequency response. 
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Figure 4.14 Processed signal to obtain the real and imaginary components 
of the frequency response. 


At the next stage, the signal was processed in order to obtain the 
real and imaginary components of the frequency response (Fig. 4.14). 
Nyquist diagram was then drawn using these components (Fig. 4.15). 
The results indicate a stable structure from the dynamical point of view, 
as expected from the theoretical analysis made by the finite element 
method. The differences between the calculated natural pulsations using 
the FEM and those experimentally obtained are minimal, as Table 4.1 
shows. 

The experimental validation of the natural frequencies and the 
implicit natural vibration modes estimated by the finite element method 
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Figure 4.15 Nyquist diagram drawn using real and imaginary components 
of the frequency response. 



Table 4.1 Differences Between Calculated Natural Pulsations Using 
FEA and Those Experimentally Obtained 


Calculated natural frequency, f c [Hz] 
Measured natural frequency, / m [Hz] 


Deviation £ = 

/t \ /c *100 


[%] 


fm 




259.6 

380 

491 

265 

374 

474 

2.04 

1.60 

3.59 


lead to the conclusion that the analysis methods implied by this research 
have been correctly applied, and the results obtained are significant. 


4.5 CONCLUSIONS AND IMPLICATIONS 
OF DIAGNOSTIC ANALYSIS IN 
FUTURE RESEARCH 

Theoretical and experimental analysis of the physical model of the kine- 
matic feed chain have yielded the possibilities of identifying the points of 
maximum stress in the structure of the kinematic chain, and identifying 
the natural frequencies of the structure. In addition, reliable stability of 
the physical model has been confirmed, which makes this an apt model 
for the study of problems related to the technical diagnosis of elements 
from its structure. 

By analyzing the deformations of the physical model the following 
highlighted that the most important yield occurs where the zone of the 
mechanism is predisposed to failures, namely, the transforming mecha- 
nism screw-nut with balls and the rolling guidings of the longitudinal 
slide. These yields can influence the functioning of the mechanism and 
speed up the occurrence of the damaging processes that contribute to 
the mechanism's failure. The deformation values are closer to those of 
the preadjusted mechanical plays in these mechanisms that can lead to 
a fault functioning. 

Knowing the natural frequencies of the structure enables us to 
compare them with the characteristic frequencies (of operation) of all 
the mechanisms in the studied kinematic chain's structure. From the 
technical literature and also from the theoretical and experimental re- 
search for the determination of the characteristic frequencies of some 
mechanisms with bearings, ball screws, and roller cam followers, it has 


been noticed that these work frequencies are in a neighboring domain 
with that of the investigated natural frequencies, but usually inferior. 

Theoretically there is a possibility of resonance occurrence with 
the natural frequencies because of the important increase in the power 
spectrum (e.g., for an element with mechanical fault) that occurs not 
only at the characteristic frequencies, but also superior harmonics of 
these frequencies. This phenomenon leads to signals of false mechanical 
faults and unreasonable stoppages of machines, and happens especially 
in online diagnosis systems having low resolution. 

The results of this diagnostic analysis have oriented, from these 
points of view, future theoretical and applicative research concerning 
technical diagnostics for some elements of the general structure of kine- 
matic feed chains. 



Designing and Manufacturing a 
Virtual Instrument Vibroacoustic 
Method Diagnostic Package 


5 . 1 VIRTUAL INSTRUMENTATION: 

LabVIEW GRAPHICAL 
PROGRAMMING ENVIRONMENT 

The virtual instrumentation for vibroacoustic diagnosis, which is elabo- 
rated and presented in this chapter, uses the LabVIEW (Laboratory Vir- 
tual Instrument Engineering Workbench) program, version 2.5, made by 
National Instruments. The software uses a graphical programming lan- 
guage, G, to create programs such as block diagrams, without sacrificing 
any of the power or performance of a traditional programming language. 
LabVIEW uses terminology, symbols, and ideas that are familiar to en- 
gineers and researchers, and it is purposely based on graphical symbols, 
rather than a textual language, for describing programming activities. 
LabVIEW combines the technology of the most recent operating sys- 
tems with specialized programming techniques (OOPT, object-oriented 
programming techniques) to obtain a simple and flexible operating en- 
vironment. 

The LabVIEW software objective is structured on several main 
levels. 




1. The acquisition module provides the possibility of programming 
and controlling the hardware connected to the measurement/control sys- 
tem. This possibility is realized using a device driver software specific to 
each connected apparatus. At the present time over 450 of such kinds 
of drivers exist, made by drivers for instruments from over 45 different 
producers. It is important to note that all these drivers have a common 
programming interface (API, application programming interface) that 
allows the applied programs to be independent of the operating or cal- 
culus system. Lab VIEW is based on the NI-DAQ (National Instruments 
Data Acquisition) driver, ensuring the data management, buffer memo- 
ries, and other resources specific to each plate. The most used acquisition 
boards provided with this driver are: 

Multifunctional: AT-MIO-16D, AT-MIO-64F-5, PC-LPM-16, Lab- 
PC+ 

With rapid analogous input: EISA-A2000 

For signal dynamical acquisition: AT-A2150 

With numerical input/output: AT-DIO-32F, PC-DIO-42, PC- 
DIO-96 

With timing input/output (counting/time basis): PC-TIO-10 

2. The processing modulus converts the primary acquisition data 
into significant results. Lab VIEW has libraries of functions and routines 
for a multitude of programming requests and also a series of specific 
applications. Lab VIEW also includes the tools necessary for building 
and/or developing applications developed by the user, visualization of 
data flow, and fixing possible programming errors. Lab VIEW software 
offers a complete and powerful analysis ensemble for numerical data 
treatment, an ensemble that contains: 

Statistical process: histograms, calculus of mean and square mean, 
distribution calculus, error functions, rational and polynomial 
interpolations, and so on 

Linear and polynomial interpolations in ID or 2D 

Linear, exponential, and polynomial regressions 

Linear algebra: vector and matrix normalization, matrix mul- 
tiplication and inversions, determinant calculus, solving linear 
equations 

Signal generation: impulse, ramp, triangle, sinus, rectangular, alea- 
tory, Gauss distribution, white noise 

Processing in the time domain; integration and differentials, deci- 
mation, cut off, limit detection, impulse analysis 



Processing in the frequency domain: rapid Fourier, Hartley, and 
Hilbert transforms, spectrogram, magnitude and power spec- 
trum, inter- and autocorrelation, convolution and de-convolution, 
phase calculation 

Windows: Hanning, Hamming, the triangular, Blackmann, Kaiser- 
Bessel, cosinus, exponential, and others 
HR filters (infinite impulse response): Butterworth, Cebisev, Bessel 
for pass-up, pass-down, and cross-band 
FIR filters (finite impulse response): Parks-McClellan, window, 
pass-up, pass-down, cross-band, and user-defined filters 
Nonlinear filters: environment 

3. The presentation modulus is of an interactive type, the com- 
puter display having the function of a front panel, similar to traditional 
instruments. Using the mouse’s help, all the buttons, switches, inverters, 
and potentiometer from the front panel can be operated. In addition a 
series of options exists for stocking data on the hard disk, for transferring 
data to a network, or obtaining hard copies on plotter or printer. 

The Lab VIEW programming environment can be used both by in- 
experienced users (because of its graphical simplicity) and by specialists 
(because of its flexibility and the possibility of developing sophisticated 
applications). 

5.2 STRUCTURE OF Lab VIEW 
VIRTUAL INSTRUMENT 

The programs elaborated using the LabVIEW programming environ- 
ment are called Vis — virtual instruments — because of their similarity to 
traditional measuring and control instruments, not only in appearance, 
but also in the operating mode. In reality, they are analogous to the 
functions from conventional programming languages. 

Virtual instruments work with graphic blocks that are sequences of 
the program presented in a graphical shape adequate for their purpose. 
A graphical block accepts one or more types of input data. These data 
are processed and then the results are sent as output usable in other sets 
of graphical blocks. In this manner the continuity of signal transmission 
is ensured. It must be noted that the route of each type of data is 
presented in a certain color and line thickness. The virtual instrument 
can be in edit mode, when the instrument is modified with specific tools, 
or run mode, when the instrument is used for the purpose for which it 
was created. 



Any virtual instrument has three components: front panel, block 
diagram, and the graphical symbol with connectors. The front panel rep- 
resents the graphical user interface. This panel simulates the front view 
of a traditional instrument and can contain buttons, switches, indicators, 
and displays for graphical representations. In addition, it must specify 
the inputs and outputs of the virtual instrument using adequate indica- 
tors and control elements. The panel is built using the Controls menu, 
where diverse options of indicators and control elements are preconfig- 
ured. Using the operating tool (the hand with the straight forefinger), 
which is operated by mouse, we can intervene on each element from the 
front panel (on-off, adjustments, rescaling). Using the positioning tool 
(the arrow) we can select, reposition, and resize the elements. 

The block diagram represents the graphical solution of program- 
ming problems. A block diagram, that is, the equivalent of the work- 
ing algorithm (of the program) for the virtual instrument, accompanies 
each front panel, and therefore the block diagram can be considered 
as a code source in the graphical programming language G. The block 
diagram is made by symbols framed in buckles (e.g., FOR, WHILE) 
and/or case structure. The components are linked in between by wires 
or nodes in such a way that they follow the data flow in the diagram 
from input to output. The block diagram is built using the Functions 
menu, to place the components (arithmetical and logical functions, FOR 
cycles, WHILE buckles, TRUE/FALSE cases), and also using the link 
tool to make the tracts between components. It must be noted that the 
Lab VIEW application automatically places on the block diagram the 
terminals associated with the indicators and control elements from the 
front panel. 

The Icon/Connector is the simplified graphical representation of 
a virtual instrument; this representation indicates only the instrument 
destination and the type and placement of its inputs and outputs. Con- 
centrating a virtual instrument in such an Icon/ Connector, it can be used 
as a subinstrument in any other block diagram connected to the other 
components of the operating schema. When they are defined, the icons 
and/or connectors are displayed in an alternative manner in the upper- 
right corner of the virtual instrument from the working window, in order 
to find the link mode between inputs and outputs. These connections 
are analogues to the parameters of a subroutine or of a function, and 
correspond to the indicators and control elements from the front panel. 

Lab VIEW is a modular and hierarchic programming environment. 
Modular because any block diagram can be concentrated in an icon and 



thus used as a subvirtual instrument (subVI), and hierarchic because 
any subinstrument may be used in the block diagram of an instrument 
of superior level. The number of virtual subinstruments from a hierar- 
chy is practically unlimited, both on horizontal and vertical. These two 
characteristics of Lab VIEW software facilitate creation, understanding, 
and maintenance of virtual instruments. A complex application may 
be divided into a series of simple tasks and a virtual instrument can 
be built to solve each of these tasks. Finally these instruments are in- 
cluded in a virtual instrument of superior level, which solves the com- 
plex application. Each instrument can be run individually, separated 
from the rest of the application, so the eventual malfunctioning can be 
rapidly located. More than that, subinstruments of inferior level can ac- 
complish tasks that are common to other applications, so they can be 
reused. 

5.3 VIRTUAL INSTRUMENTATION FOR 
SURFACE DIAGNOSIS 

As presented above, surface diagnosis methods are useful for supervis- 
ing the operating state of machine tools; they signal the overpassing of 
some signal levels which are considered normal, without providing other 
information concerning the cause or the nature of the abnormality and 
without evaluating the time remaining until the final stop. 

Three virtual instruments have been elaborated and tested for sur- 
face diagnosis: 

FACVARF, based on evaluation of the peak factor (Fv) 
INDDIAG, based on evaluation of the diagnostic index [K(t)\ 
INDKURT, based on evaluation of the captured vibrational signal 

The theoretical bases of diagnosis methods utilized were presented in 
Chapter 3 of this work, and are not redescribed here. 

Being designated for “offline” type surface diagnosis, these virtual 
instruments operate in a working space different from the space where 
captioning of signals is done. Thus grouping of the acquisition files has 
been in view, depending on measuring points; that is why temporary 
positioning of each signal acquisition has also been in view, so that the 
file format saves both the signal and the moment of acquisition. 

The three instruments can be used separately or simultaneously for 
the analysis of the same signal; they do not exclude each other — rather 
they are complementary. Following the idea of complementarity, the 


vibration signal reading zone works as a function of time and as a func- 
tion of finding out the characteristic values common to the three virtual 
instruments (Fig. 5.1). Thus, in a FOR cycle, on the basis of the path 
built by using the graphical block 2, the hard disk is read, and block 1 
(List Directory) returns two matrices that specify the names of acquisi- 
tion directories and folders met. Graphical block 3 (File/Directory Info) 
provides information concerning the directory or the folder (such as size 
and the date of creation or last modification); the date in absolute time is 
converted into a calendar date through block 4 (Get Date/Time String); 
this is used to order the function of time for different acquisitions. 

The graphical block 10 (Citvib.vi) is a specially built virtual in- 
strument, here as a sub- VI, necessary for acquisitions reading as func- 
tions of time; it also has the possibility of amplifying the signal by a 


DIRECTOR 



Figure 5.1 Vibration signal reading zone is common for FACVARF, 
INDDIAG, and INDKURT virtual instruments. 





factor introduced by the operator. Figure 5.2a presents the block dia- 
gram of this virtual instrument; the icons and connectors are shown 
in Fig. 5.2b. 

In the inferior zone of cycle FOR’s windows the arrangement of the 
work matrices is accomplished with the help of graphical blocks 6 (Array 
Max & Min) and 7 (Max & Min). The matrices are then assembled 
in a single cluster, near the information referring to date in graphical 
block 8 (Bundle). The continuation of processing the acquisitioned signal 
is different for each of the three aforementioned virtual instruments, 
depending on the diagnostic method used. 

The FACVARF.VI instrument (block diagram presented in Fig. 5.3) 
takes the data from the sorting block 11 (Sort ID Arry), and displays on 
the front panel (Fig. 5.4) the graphic of two recordings of the processed 
signal, in the time domain, as follows: the recording considered as refer- 
ence (having the machine in a perfect state of operation) and the current 
recording. For each of the two recordings the maximum value (RMS) is 
indicated; the value of the peak factor was calculated previously and 
thus, a first evaluation of the functioning state of the supervised ele- 
ment became possible. Data processing is possible by disassembling the 
information in graphical block 10 (Unbundled) and then connecting it to 
the blocks that make the link with the front panel (DBL). The ensemble 
of statistical processing of the signal is integrated in a WHILE buckle. 
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Figure 5.2 (a) Block diagram of the virtual instrument; (b) icons and 
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Figure 5.3 FACVARF.VI i 
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Figure 5.4 Front panel of FACVARF.VI instrument. 







A third graphic from the front panel of the instrument follows the evo- 
lution in time of the peak factor. A point is marked on the tendency 
graph for each new acquisition, illustrating a short history of the oper- 
ating state, as this state is captured in the respective measuring point. 
The graphic is scalable depending on the measuring unit. 

The INDIAG.VI instrument (block diagram in Fig. 5.5) performs a 
processing series similar to that presented previously, but this time using 
the algorithm specific to the diagnosis index. The front panel (Fig. 5.6) 
also contains the reference recording and the current recording, near the 
graphic of the evolution of diagnosis index. Both instruments present 
a short legend on the front panel next to the graphic of evolution of 
diagnosis index to make the reference and the current recording easier. 

The INDKURT.VI instrument is based on the diagnosis algorithm 
specific to the Kurtosis method. It is more complex compared to the 
other virtual instruments. In the block diagram (Fig. 5.7) the instrument 
contains a sub- VI (Kurt /32.VI), integrated in a WHILE buckle, for the 
calculus of the fourth-order statistic moment (/ 32 ) and the graphical 
representation of the probability distribution density of the captured 
signal. 

The symbol and connectors of the graphical block Kurt /32.VI are 
defined according to Fig. 5.8a. The calculus relation of the statistical 
moment is shown in Section 3.2. 2.3, and for its instrumentation a few 
virtual instruments from the Lab VIEW software library have been used 
(Fig. 5.8b). Those virtual instruments are the following: Standard De- 
viation. VI, to calculate the arithmetic mean of the signal and standard 
deviation; Gaussian White Noise. VI to test the instrument; Numeric 
Integration. VI, to integrate the entire domain; and Histogram. VI, to 
calculate the probability density of the same signal. Figure 5.9 presents 
the block diagram of the Kurt /32.VI virtual instrument. 

The INDIKURT.VI virtual instrument displays (Fig. 5.10) a table 
of files on the front panel; those files have the available acquisitions on 
the stocking directory, listed in terms of the acquisition data. When a 
file is selected from this table certain data processing is performed: 

Representation of the captured signal, as a time function (left-up 
window) 

Calculus of the maximum and effective values for the peak factor 
and Kurtosis index 

Graphical representation of the probability density of the signal 
(left-down window) 
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Figure 5.6 Front panel of INDIAG.VI instrument 
















Figure 5.7 INDKURT.VI instrument — block diagram 














Figure 5.8 (a) Symbols and connectors of graphic block Kurt /32; (b) vir- 

tual instruments from the LabVIEW library used to calculate the statistical 
moment. 

Knowing the fact that a mechanism in good functioning state is 
the source of some stochastic vibrations that respect a normal (Gauss) 
distribution of amplitude (in this case the fourth-order statistic moment 
is placed in the range of value 3), an evident correlation can be made be- 
tween the increase of the fa index value and the deviation of probability 
from the Gaussian shape, when a defect occurs and develops itself. 

Taking into consideration the previous determinations, it has been 
noticed that the first instrument (which notes the abnormal state pro- 
duced by the occurrence and development of a damage) is INDKURT; 
confirmation is done gradually by INDIAG and FACVARF. 

5.4 VIRTUAL INSTRUMENTATION FOR 
PROFOUNDNESS DIAGNOSIS 

A precision vibroacoustic diagnosis cannot be accomplished using only 
statistical parameters associated with signal processing in the time do- 
main. The most rapid solution to avoid this deadlock is to pass the signal 
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Figure 5.9 Block diagram of Kurt (32 virtual instrument. 
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Figure 5.10 Front panel of INDIKURT.VI virtual instrument. 




in a frequency range using a Fourier transform. Interpretation of the fre- 
quency spectrum produces, as seen in previous chapters, important data 
about the causes and the nature of the damage. 

A first virtual instrument created for profoundness diagnosis, named 
MECOSPEC, has at its base the spectrum comparison method, and it 
has certain advantages: 

The diagnostic method is included within the category of those 
which ensure from the beginning an increased precision of defect 
identification. 

Identification of frequencies characteristic at defect is made auto- 
matically, excluding the subjective factor. 

It is possible to estimate the reserve time interval until the machine 
comes to a stop because of a defect. 

Acquisition intervals can be programmed with a fixed or variable 
step. 

A diagnosis specific database can be created. 

Although sophisticated, the method is simple to use. 

Because of the complexity of virtual apparatus created to apply this 
method, the data processing has been organized on working blocks with 
the vibration signal. 

The “acquisition and calculus of power” block, whose block dia- 
gram (partial) is presented in Figure 5.11, makes the link between the 
AT-MIO-16F acquisition plate and computer. 

Captured signals using the Ach3p.VI virtual instrument are called 
and processed by another virtual instrument (Files. VI) especially cre- 
ated to manage the data files, whose front panel (Fig. 5.12) is made like 
a record of the measuring points. The virtual instrument associates a 
series of characteristics with each measuring point of the machine; those 
characteristics are useful for storing captured data such as point in- 
dex, number of averages and averaging control, sensibility of acquisition 
transducer, and spectrum type (acceleration, speed, or displacement), 
as can be seen in the block diagram (Fig. 5.13). This supervising record 
must be completed by the user (operator). The signals are then passed 
through a Hanning window (Hann 2 A 17.VI) to avoid a cutting-up effect 
when the signal passes from the time domain to the frequency domain, 
and the power spectrum is obtained using the rapid Fourier transform 
(FFT). 

The front panel of this first signal processing block (Fig. 5.14) con- 
tains a graphical representation in the time domain of the captured vi- 
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Figure 5.11 Partial block diagram of “Acquisition and calculus of power spectrum” block. 
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Figure 5.14 Front panel of first block for signal processing. 




bration signal, and also contains the power spectrum graphic of the same 
signal, on an adjusting domain between 10 Hz and 100 kHz. In order to 
make possible the reading of values from the two graphics, CO pointers 
can be used to scavenge the graphics (using the four buttons disposed 
as a rhombus). The abscissa and the ordinate of the point indicated by 
the pointer are numerically displayed in the foot of graph. The right 
column groups the elements for adjusting the represented parameters 
influencing the block diagram. 

The block “Comparison of power spectra” is presented as a block 
diagram in Figure 5.15. Comparison of spectra has as its starting point 
the assignment of the reference spectrum obtained with a machine in 
a perfect operating state; all of the spectra obtained from the signals 
captured previously are compared with this power spectrum. In order 
to avoid the false alarms caused by possible spectrum slipping (e.g., at 
small variations of working rpm), the reference spectrum is boarded 
in a virtual subinstrument (Bord 3p.VI), whose block diagram is pre- 
sented in Figure 5.16. The boarding width can be manually adjusted 
from the block’s front panel. Transformation of the reference spectrum, 
after boarding, in the warning and alarming patterns, respectively, is 
made by translating the reference spectrum to adjustable levels using 
the buttons of the front panel. The 2x levels are recommended for the 
warning pattern and the 10 x for the alarm pattern. 

The front panel (Fig. 5.17) “Comparison of power spectra” block 
contains two graphical representations as well: upwards, the currently 
captured spectrum superposed on the boarded reference spectrum; down- 
wards, the difference derived as a result of comparison between the two 
spectra. Using the adjusting elements from the right column, the board- 
ing width, warning level, and alarm level can be varied. 

For the profoundness diagnosis, as presented in Chapter 3, a method 
with excellent results in identification of defects of mechanisms with 
rolling elements is the evolution method. A virtual apparatus that sim- 
ulates the application of this method is presented in the following para- 
graphs. 

Figure 5.18 presents a FOR buckle as a block diagram where two 
sinusoidal signals are composed and modulated. After the modulation 
an exponential damping is applied to the resulting signal. A pulse sig- 
nal is obtained whose amplitude and frequency are adjustable from the 
front panel's buttons (Fig. 5.19) and which simulates the excitation sig- 
nal introduced by the occurrence of a defect (running of ball or roller 
over irregularities/pinches occurring on one of the running tracks). The 



Figure 5.15 Block diagram of “Comparison power spectra” block. 
















Figure 5.16 Block diagram of Bord.3p.VI virtual instrument. 

impulses are identical and evenly distanced at period T. The spectrum 
of this series of impacts would be a spectral line that includes all of the 
harmonics of the repetition frequency 1/T. These harmonics have the 
largest amplitudes in the vicinity of the resonance frequency and their 
diagnosing by simple evaluation of spectra is difficult. 

In the next step, the critical zone, which contains the structural 
resonance that has been excited by the impact due to defect, is ex- 
tracted from the frequency spectrum by zoom-in. Applying the Hilbert 
transform, evolution of the signal in the time domain is generated by 
the calculation of the time function amplitude. This evolution can be 
passed and analyzed in the frequency domain to establish the frequency 
or frequencies of impact. On this basis the type of defect and the place 
of occurrence can be identified, as well as its dimensions. 

The instrument’s front panel offers the representation of the sim- 
ulated pulse signal (upward graphic) and the corresponding evolution 
from the time domain (downward graphic). The processing of this last 
signal in the frequency domain is similar to what was done in the “ Acqui- 
sition and calculus of power spectrum” of the MECOSPEC apparatus, 
and it has not been included in the simulation. 

Results from the work presented so far show that the evolution 
method has large perspectives in identification and separation of the ex- 
istent modular sources in a machine tool; it allows a precise diagnosis of 
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Figure 5.17 Front panel of “Comparison power spectra” block. 




Figure 5.18 Block diagram of METINFA 



/irtual apparatus. 
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Figure 5.19 Front panel of METINFA virtual apparatus. 




defects in bearings, rolling guides, and gears, even though many defects 
are simultaneously present. 

These virtual instruments have been tested in the laboratory and 
then compared with traditional measurement schemes, built with tra- 
ditional instrumentation. For example, the virtual instrument for the 
acquisition and calculus of the power spectrum has been mounted in par- 
allel with the real-time frequency analyzer Briiel & Kjaer 2034, belonging 
to the Eurotest laboratory, in order to process the same vibration sig- 
nal. A perfect similitude has been noticed between the displayed power 
spectra, with a better resolution of the virtual instrument for lower fre- 
quencies. All of the virtual instruments described have been found to be 
adequate from the functional and operating accuracy point of view. 



Theoretical and Experimental 
Research: Diagnosis of Some 
Feed Kinematic Chain 
Structure Elements 


6.1 VIBRATION AND NOISE SOURCES 
DURING FEED KINEMATIC 
CHAIN OPERATION 

Modern mechanical systems have necessitated the development of re- 
search and creation of performance mechanisms. These mechanisms are 
capable of reducing stresses and energy loss due to friction, and produc- 
ing a substantial increase of output. Integration of these mechanisms in 
the kinematic chains of machine tools leads to performance increase (pre- 
cision, dynamical behavior, and reliability) of these chains and further- 
more to the possibility of their integration in modern flexible production 
systems. 

The helical joint with rolling elements has an outstanding place 
among the above-mentioned mechanisms. Introducing the rolling bodies 
between screw and nut, the sliding friction is replaced by rolling fric- 
tion; as a result, positioning precision and output increase. Noise and 
vibration accompany the functioning of feed kinematic chains with the 
movement-transforming element-type screw-nut with balls. The noise 
and vibration, together with the other sources, can negatively influence 




the manufacturing process, even though they are not very powerful. The 
noise and vibration can also negatively influence the dynamic and reli- 
ability of machine tools. 

Figure 6.1 presents the main noise and vibration sources that can 
appear during the functioning of modern feed kinematic chains, grouped 
in internal and external sources depending on the origin point. In the 
case of external sources it can be stated with certitude that the most 
important weight belongs to the cutting process. In the case of internal 
sources, only a careful analysis can disclose the weight of each source in 
the global level of noise and vibrations. 

A first step of this analysis is the determination by calculus of 
the characteristic frequencies for each mechanism that is considered as 
a possible source. It follows the recording of vibroacoustic signals com- 
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Figure 6.1 Main noise and vibration sources that can appear during func 
tioning of modern feed kinematic chains. 







ing from the previously mentioned sources, and the treatment of these 
sources using an FFT algorithm. The obtained spectrograms can indi- 
cate the weight of each characteristic frequency, and, by interpolation, 
the importance of each source in the assembly of the feed kinematic 
chain. 


6.2 CALCULUS OF THE CHARACTERISTIC 
FREQUENCIES OF BEARINGS 

A few simplifying hypotheses are necessary to calculate the characteristic 
frequencies from bearings: 

The contact between the elements and rolling paths is done without 
slipping. 

The influence of inertial and weight forces is negligible. 

The effects of centrifugal force and gyroscopic couple are ignored. 
The pressure angle olb is considered the same for the inner and the 
outer rings. 

If the O point from the bearing axis (Fig. 6.2) is taken as the 
reference point, each point of contact between the rolling body and the 
two rolling paths has an angular speed equal to that of the median point 
of the body; that is: 



Figure 6.2 Bearing schematics. 


The rotational speeds of the contact points can be written as 


Va = -^WMo{Dm + Db cos a B ) 

x (6-2) 

Vi = -LOMo(Dm — Db COS OLb) 

If we write the angular speed of median point M function of O point as 


umo = Vm/ ( 


then the above relations become: 


Ve 
V i 




Db_ 

Dm 


cos a e 




Dm 


cos as 




The E point is the center of instantaneous rotation of rolling bodies and 
each point of these bodies rotates around E with the angular speed oje- 
The M point results from the following formula for the rotational speed, 


Vm = ~^ADb 

(6.5) 

and for the I point 


Vi = ^ADb 

(6.6) 


This speed should be equal, according to the first hypothesis with the 
peripheral speed of the inner rolling path in the I point: 

Vj = Trn(Dm ~ Db cos a#) (6.7) 


From the last three relations, by substitution, the rotation speed of the 
median point M results. This speed has the expression: 



The rotational peripheral speed of the bearing cage is equal to the speed 
of the median point M from the rolling bodies, V c = Vm- 

Coming back to the expressions of the rotational speeds of the 
points from the rolling paths, inner and outer (6.4), and replacing the 



median point speed with the previously obtained relations (6.8), the new 
expressions result: 


Ve = -tttiD 
2 


m 1 I t-) COS Ot B 

\ i-/m 


D b 


Vi = - 7 rnD 


1 - -pp COS a B 


(6.9) 


In this moment, the rotational frequencies characteristic to rolling bear- 
ings can be deduced: 


Rotational frequency of the inner ring (rotational) 

fi = fn = ™/60 

Rotational frequency of the cage: 


( 6 . 10 ) 


fc — 2 f n 


i Db 

1 — — — cos a b 


( 6 . 11 ) 


Rotational frequency of the rolling bodies in the cage 


f _ 1 .D m 

Jm ~ 0 Jn j-, 

Z Lf h 


1 ( Db 

1 I jz— cos a B 
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( 6 . 12 ) 


Rotational frequency of the rolling bodies on the outer ring: 


JE = ^JnZ ( 1 - £j— COS a B 


Rotational frequency of the rolling bodies on the inner ring: 


Ji = 2 f^ z { 1 + C ° S aB 


(6.13) 


(6.14) 


where z is the number of the rolling bodies. 

In the case of using balls as rolling elements, an experimental 
relation exists that also permits the determination of their resonance 
frequency: 


fb = 


0.848 E 

D b 2 p 


(6.15) 


where E is the elasticity modulus and p is the density of ball material. 
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Figure 6.3 Vibration signal of bearing located in front of the lead screw, in 
perfect functioning state, recorded in time domain. 


Figure 6.3a shows the vibration signal of the bearing located in 
front of the lead screw, in perfect functioning state, recorded in the 
time domain. Passing this vibration signal through a FFT of order 10, 
a spectrum is obtained in the frequency domain. Figure 6.3b shows this 
spectrum. 

6.3 CALCULUS OF FREQUENCIES 
THAT CHARACTERIZE THE BALL 
SCREW-NUT MECHANISM 

In the context of theoretical research regarding the vibroacoustic diag- 
nosis of machine tools integrated in flexible manufacturing systems the 
following problem was posed: knowledge of characteristic frequencies of 
the main movement transforming mechanism in feed kinematic chains — 
the screw-nut with balls mechanism. 

Some profoundness diagnosis methods (i.e., the evolute method) 
stem from the idea that a series of mechanical effects is manifesting in 
the frequency spectrum by large band increases near the characteristic 
frequencies and their harmonics from the vibrational signal. It is neces- 
sary to previously determine by theoretical procedure the characteristic 
frequencies, as they are not always easy to find in the standard fre- 
quency spectrum. In order to theoretically determine the characteristic 
frequencies an original calculus method has been set. 

The study of ball movement on the rolling path is done starting 
from some simplifying hypothesis similar to the hypothesis presented in 
the previous paragraph: the contact between the balls and the rolling 
paths is done without sliding; and the gyroscopic couple influence, and 
those of inertial, weight, and centrifugal forces are ignored. To study 
the complex character of ball movement in the rolling channel cre- 
ated between screw and nut, a fixed triorthogonal axis system OXYZ 
(OZ direction superposes the screw axis) and a mobile system oxyz are 
considered. Regarding the mobile system, the oy direction is tangent to 
the helix described by the ball center during movement (Fig. 6.4). In the 
general case of ball movement on the rolling canal, it is considered that 
this canal is made from two trays of ogival profile, driven by axial forces 
of opposite directions. Two contact points will exist between the ball 
and the interior tray (the leading screw's canal) and one contact point 
between the ball and the external tray (the nut’s canal). These contact 
points determine a plane that also passes through the center of the ball. 
The rolling movement of the ball will be done also around the rela- 



Figure 6.4 Helix described by the ball center during movement in the screw- 
and-nut mechanism. 


tive instantaneous rotation axis (A) that passes through the previously 
mentioned contact points of the ball with the leading screw (Fig. 6.5), 

P cos#, - 7 p sin#, 6 ^^tga\ (6.16) 


where D m is the medium diameter, measured to the ball center. 

The tangent to the helix in the corresponding point is expressed as: 


- dP / d6 ( • /> * • x 

r = — = ( — sin # cos a. cos # cos ot, sin ct ) 

\dP/dO\ 


(6.17) 


and the interior normal is: 

dr/de ■ a a\ 

n = WM = { ~ ' 


(6.18) 


For the mobile axis system Oxyz, the vector of which the module 
is 1 for the Oy axis is r, for the Ox axis is — n, and for the Oz axis is 
V = —nxr. 

Let us consider now S\ and S 2 the contact points between the ball 
and the leading screw, S 3 the contact point with the nut (Fig. 6.6), and 



Qm*Dbtos£ 






D b the ball's diameter. The location of these points depending on the 
mobile reference point is given by the relation: 


rsi = — y cos A” + -y sin A* 7 


(6.19) 


Depending on the fixed reference point, the position is Rsi = R + fsi, 
and the coordinates can be written as 


X Si = 


s\ . -^b q s\ . -^b • n • 

— — cos 6 + — cos pi cos 0 + — sm Pi sin n sin 6 

Z Z z 

Ysi = sin 0 + Db cos Pi sin 0 — ^ sin Pi sin a cos 9 


( 6 . 20 ) 


Zsi — 0 — c sin Pi cos a 


For Si and S 2 points from the screw, the speed is 


USi = Xsi 


= 0 


-Dtti • n -Db o ' r\ 1 ~D b 'O' r\ 

— — sm 6 — cos Pi sm 9 + sm Pi sm a cos 0 


VSi = Ysi 


( 6 . 21 ) 


= 0 


■Dtti q -D b 'O’ • /) 

— — cos 0 H cos Pi cos 0 + — — sm ft sm a sm 6 
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u>Si = Z S i = 0—^tga 


Simplifying the calculus for 0 = 0, and knowing that 0 = the follow- 
ing speeds result. 

D b . 

usi = u s — sm Pi sm a 


. Db a 

VSi = ^5 I -W- + “W COS ft 


( 6 . 22 ) 
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VJSi = ^5 


m 


tga i = 1,2 


The Si and S 2 points are usually symmetric related to the Ox axis, 
so it can be written: 


Pi = 7T - ft ft = 7T + ft ft = 27T - P' 


(6.23) 



and the previous relations become: 
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(6.24) 


U 5 2 = ids 
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cos /3 


W S 2 = ids-^-tga 

The S3 point belonging to the nut has no speed (the screw is considered 
in a rototranslation movement). If C is the center of a ball, then its 
vector-radius is: 


R c — Aiifcl + ^2-^5 2 + \zRs3 


(6.25) 


where the coefficients have the values 


. /? + /?' 

sin — - — 


Ai — 


2 sin /3 cos 




. /?-/?' 
sin 


(6.26) 


Ao — 


2 sin P cos 


/? + /?' 


A 3 — 


sin ft 


cos (3 — cos P' 


The speed of the C center can be written as a vector 
v c = A1T5I + A 2 Ts 2 


(6.27) 



or detailed on components: 


D b cos f3' sin /3 sin a 
U ° UJs 2 cos /3 + cos P' 



D b \ cos P' 

— cos p 

2 J cos P + cos P' 



Dm COS P' 

2 g cos P + cos P' 


(6.28) 


Considering the simplifications tg 2 a = 1 and sin 2 a = 1 (a < 1), the 
value of the ball center speed becomes: 




D b cos P) 


cos P' 

2 (cos P + cos P') 


(6.29) 


The relative speed of the C point in a reference point with the origin in 

— > — * — * 

Si, which is united with the screw, is the difference AV = Vc — Vs 1, 
and the value of this speed is 

|AV| = u,(D m - Ptcos/?) 2(cos °“' 3 cosg , ) (6.30) 

Knowing that the balls are consecutive, the characteristic frequencies 
can be determined at the moment of: 


The rolling frequency of the balls on the leading screw tray: 



w s (£>m - D b COS (3) 
2D b (l + cos /3') 


(6.31) 


The rolling frequency of the balls on the nut’s tray: 



U s (Dm - Db cos/3) 
2Dfr(l + cos P) 


The rotation frequency of the balls themselves: 


fb = u s 


COS P' 

COS P + COS P' 


(6.32) 


(6.33) 


Figure 6.7a presents the vibration signal of the screw-nut with balls 
mechanism, recorded in the time domain. The signal was then passed 
in the frequency domain through a FFT of the tenth order. Figure 6.7b 
presents the signal’s frequency spectrum. 
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Figure 6.7 Vibration signal of the screw-nut with balls mechanism, recorded 
in time domain. 



6.4 CALCULUS OF CHARACTERISTIC 
FREQUENCIES FOR THE TANKETTE 
WITH ROLLS 

In the case of tankettes with rolls, the single simplifying hypothesis refers 
to the existence of the nonsliding contact between the rolls and the 
rolling paths. 

Figure 6.8 shows that the tangential speeds of the rolls on the 
rolling paths are equal and of opposite direction, and in direct relation 
with the displacement speed of the mobile element (the saddle): 

V s = —Vi = V ( 6 . 34 ) 


In this case, the angular speed of the balls is 


2V 


lj = 


D 


( 6 . 35 ) 
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which leads to a frequency 


fb = 


V 


7T D r 


( 6 . 36 ) 


The number of the active rolls is approximated as 
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( 6 . 37 ) 
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Figure 6.8 Tangential speeds of the rolls on the rolling paths are equal and 
of opposite direction. 




In the above-mentioned conditions, it can be approximated as the rolling 
frequency of the rolls on the rolling paths (inferior and superior paths) 


_ VL 

Js — Ji — Jb^a — 

71 U, 


or the specific rotation frequency of the rolls: 


fr = 


V 


71 D 


Observations: 


(6.38) 


(6.39) 


1. The identification on the basis of previous relations of a dam- 
age on the rolling element is burdened by the fact that the 
roll has a passive stroke bigger than that of the active stroke: 
Lp = L + 7t(Dq + Db) (relative to the roll center), whose vi- 
bration signal is absent. 

2. The identification of the damage from the lower rolling path 
(guiding) is possible only in the time interval when the guiding 
is covered by the tankette. 

Figure 6.9a shows the vibration signal of one of the tankettes with 
rolls of the stand, signal recorded in the time domain. The signal was 
then passed in the frequency domain by a tenth-order FFT. Figure 6.9b 
presents the frequency spectrum of the signal. 


6.5 PRELIMINARY EXPERIMENTAL 
RESEARCH 

As shown in previous chapters, the main sources of noise and vibration 
in feed kinematic chains of modern construction are the bearings of the 
leading screw, the mechanism of movement transformation-type screw- 
nut with balls, and the rolling guidings of the saddle. At the same time, 
in the diagnostic analysis made in Chapter 4 of this book, the main 
sources of noise and vibration were identified as being the weak points 
from the feed kinematic chain structure. The diagnosis of the mechanical 
system of the feed kinematic chain must designate these sensible points 
exactly. 

The common feature of these mechanisms is that they have the 
same principle of operation: transformation of sliding friction into rolling 
friction inserting intermediary elements — the rolling bodies. More than 
that, the damages of these mechanisms during functioning have, most 
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Figure 6.9 (a) The vibration signal of one of the tankettes with rolls of 

the stand, signal recorded in the time domain; (b) frequency spectrum of the 
signal. 




Figure 6.10 Causes of damage to bearings and leading screw. 


often, the same causes as those generically presented in Figure 6.10. 
These common points allow a unitary treatment from the viewpoint of 
the diagnostic methods used. Thus, this chapter presents applications 
of the monitoring and diagnosis methods to one or another analyzed 
mechanism that does not exclude the use of the same methods to other 
mechanisms of the feed kinematic chain. 

6.5.1 Experimental Research with General 
Use Instrumentation 

The first determinations accomplished in the frame of this work were 
done with a general use instrumentation, of portable type, having 
medium performances, manufactured by RFT (Germany) and B&K 
(Denmark). Functioning state supervision of the 6207-P6 bearings on 
the test stand was done on a bearing testing machine whose kinematic 
schema is presented in Figure 6.11. An accelerometer B&K 4344 was 
connected on the cover of the tested bearing. The B&K 4344 accelerom- 
eter was connected to a general use B&K 2511 vibrometer. Peak and 
effective values were measured at regular intervals, each time calculating 
the peak factor (peak factor method). Figure 6.12 presents the results. 
When the bearing was dismounted, the occurrence of pitting phenom- 
ena on the rolling path of the inner ring could be noticed after the peak 
factor became larger than 23. 



Figure 6.11 Kinematic schema of the bearing testing machine. 
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Figure 6.12 Peak factor method. 










In parallel with measurements of vibrations, the noise emitted by 
bearings was also studied. For this purpose, a reverberant chamber was 
created setting panels around the testing machine for bearings. An RFT 
0023 sonometer with a microphone with capacitor was also used. The 
preliminary calibration of the microphone was done with the pistophone 
RFT 0008. A good correlation was highlighted this way (between the 
limits ±2 dB) between the level of acceleration signal of vibration and 
the emitted noise, as can be observed in Figure 6.13. 

The technical literature claims that the peak factor method can 
not precisely locate the damage because it is strongly influenced by the 
elements of the mechanical transmission of movement. This statement 
was tested by remounting the fault bearing on the test stand, at the 
entering end of the leading screw. Vibration measurements were done 
using the same equipment as used for the previous experiments (peak 
value and effective value of acceleration signal had been done), to both 
bearing cases of the leading screw. The peak factor was calculated, with 
the result that the difference between calculated value for the case with 
the fault bearing and the case with the good bearing was only of one 
unit, which confirmed the imprecision of source localization using the 
peak factor method. 
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Figure 6.13 Good correlation between level of vibration acceleration signal 
and emitted noise. 
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Figure 6.14 Method of diagnosis index has indicated the appearance of the 
injurious process and of the damage before the peak factor method. 

After the experimental determinations, it can be noticed that the 
peak factor method is simple, easy to handle, and a method that uti- 
lizes general equipment. This method is useful especially in the case of 
monitoring a large number of measuring points, when an early warning 
is not required, and the consequences of the damage are not consider- 
able. In particular situations, the peak factor method can be very well 
completed by a profoundness diagnosis method such as, for example, 
Cepstrum analysis. 

The diagnosis index method was applied in parallel with the peak 
factor method on the 6207-P6 bearing tested on a bearings 5 testing ma- 
chine. A first determination of statistical magnitudes that characterize 
the vibration signal for the bearing in perfect functioning state was useful 
to identify the reference and effective values of the signal. It is necessary 
to know these values because the diagnosis method is of “with normal- 
ization 55 type, which means the connection to the reference values. The 
diagnostic index method has indicated the appearance of the injurious 
process and of the damage before the peak factor method (Fig. 6.14). 

6.5.2 Experimental Research Using 
Frequency Analysis 

A new series of experiments was done after buying an assisted system 
of data acquisition, which is a process manufactured by the firm S.C. 
Emco S.A. This system captures vibration signals and process and has as 




components (Fig. 6.15) a commanded load amplifier (A), for transducers 
of piezoelectric type, a block of simultaneous acquisition with 16 channels 
(B), and an interface of 8 bits for a computer (C) having a minimum 
configuration of 286 XT, 1 Mb RAM, and a VGA monitor. The software 
that governs this system (VIB-01, then VIB-02 version) was adapted to 
the requests of monitoring and diagnosis together with the specialists of 
Emco. 

Hypersignal application of the firm Hyperception was adapted later 
to the above-mentioned system only for signal processing. 

This previously mentioned software is powerful enough to allow 
processing the signal in the time domain, passing the signal in the fre- 
quency domain using a fast Fourier transform and multiple transforma- 
tions in the frequency domain. 

Several options exist in the time domain: display the captured sig- 
nal (waveform display or digital oscilloscope); define and use FIR (fi- 
nite impulse response) and HR (infinite impulse response) filters; use of 
convolution, correlation, and autocorrelation functions; and fast Fourier 
transformation; and in the frequency domain: amplitude representations 
and phase-frequency representations, respectively; power spectrum; spec- 
trograms in two and three dimensions, respectively; pole-zero diagram; 



Figure 6.15 System of data acquisition. 


transfer function for two signals; and inverse fast Fourier transformation 
(IFFT). 

The Emco system of data acquisition and process, assisted by com- 
puter through VIB-02 application has been used for detection of a dam- 
age on the guidings with roll tankettes of the stand, by the Kurtosis 
method (Fig. 6.16). The simulated damage was the occurrence of pit- 
ting phenomena on the inferior rolling path of the tankettes. A B&K 
4344 accelerometer was fixed on the longitudinal saddle in the near- 
est neighborhood of a tankette. The electronic block that processes the 
captured signal to calculate the Kurtosis index value was presented in 
Chapter 3 of this book (Fig. 3.4). A block of graphical process was at- 
tached to the electronic block to visualize the vibration signal and the 
probability density of the signal’s amplitude. 

Figures 6.17a and b present the vibration’s acceleration and the 
probability density of the amplitude of this signal for guiding in a good 
functioning state. Figures 6.18a and b present the same thing as Fig- 
ure 6.17 but for the damaged guiding. The modifications of the accelera- 
tion signal of vibration, and, more important, of the probability density 



Figure 6.16 Emco system of data acquisition and process, assisted by com 
puter through VIB-02 application. 
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Figure 6.17 Vibration’s acceleration and amplitude probability density of 
this signal for guiding in a good functioning state. 


function can be monitored. In relation to this function, the calculated 
Kurtosis factor is, in the first case, f3 2 = 3.48, and in the second case 
(3 2 = 9.34. 

Determination by calculus of characteristic frequencies of some el- 
ements from the kinematic chain structure of machine tools allows the 
identification of certain types of faults by analysis of the frequency spec- 


A* io> 

4.10 
2,28 
0.46 
- 1.36 
-3.18 
-5.00 



* to * 
4,55 - 



1 J 

- 3,02 0.96 * 10 * 


Figure 6.18 Vibration’s acceleration and the amplitude probability density 
of this signal for the damaged guiding. 



trum. As presented before, the most useful representation in this case 
is the power spectrum, which is correlated directly with the vibration 
signal energy. 

The bearing testing machine presented in Figure 6.11 was used 
again to identify the faults that occur in the functioning of radial ball 
bearings, through power spectrum analysis. This time, the Emco system 
with the Hypersignal software was used to capture and process the vi- 
brator signal from the case of the radial single-row ball bearings 6207 
and 6209. 

Table 6.1 presents the constructive parameters considered for the 
calculus of characteristic frequencies of 6207 and 6209 bearings, and also 
the values of these frequencies, depending on rotational frequency f n of 
the inner ring. The characteristic frequencies were calculated using the 
relations (6.11) to (6.13) previously presented. 

Figure 6.19 presents a power spectrum that was drawn on a 6207 
bearing in a state of perfect functioning. This power spectrum can be 
considered as the reference spectrum. After rolling a few days in the 
absence of any lubricant, the first signs of damage occurred and the 
interpretation of the power spectrum of the captured vibrational signal 
could be started (Fig. 6.20). It can be noticed that the peak of greatest 
amplitude corresponds to 78.1 Hz frequency. 

Considering that the rotational speed of the inner ring is equal 
to the rotational speed of synchronism of the driving electrical motor 
iP'ME — 1420 rpm), the frequency of the inner ring is f n = 23.66 Hz. 
In these conditions, the frequency characteristic to balls rolling on the 
external tray is equal to f e = 77.42 Hz, calculated with the relation to 
Table 6.1. The calculus of the other frequencies characteristic to bearing 
functioning did not offer any similitude with the significant increases of 
amplitude in the power spectrum. It can be stated that the frequency of 


Table 6.1 Constructive Parameters Considered for the Calculus of 
Characteristic Frequencies of 6207 and 6209 Bearings 



Constructive parameters 

Characteristic frequencies 
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fe 
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Dm 

D b 
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(fn) 

(fn) 

(fn) 

(fn) 

6207 

54.5 

9.922 

8 

0 

3.272 

4.728 

2.655 

0.409 

6209 

65 

12.000 

8 

0 

3.262 

4.738 

2.616 

0.408 


POWER SPECTRA ANALYSIS 

ffa&vjjmf FRAME 0 C0»08fl ) 



FREQUENCY 

Figure 6.19 Power spectrum drawn on a 6207 bearing in a state of perfect 
functioning. 
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Figure 6.20 First signs of damage on the power spectrum drawn for the 
6207 bearing. 



balls rolling on the external tray and its superior harmonics that occur 
in the last captured power spectrum signalize the beginning of a fault on 
the external rolling path — a fact confirmed by the subsequent evolution 
of the fault. 

The same identification method of mechanical defects was applied 
also in the case of the movement-transforming mechanism screw-nut 
with balls, for a damage simulated on the rolling path from the leading 
screw. 

Table 6.2 presents the values of the constructive parameters that 
enter in the calculus of characteristic frequencies of the screw-nut with 
balls mechanism [relations (6.31) through (6.33)], and also the values 
of these frequencies, depending on the rotation frequency of the leading 
screw. 

Figure 6.21 shows a power spectrum resulting from a vibration sig- 
nal captured on the “clean” zone of the leading screw using an Emco ac- 
quisition system and Hypersignal software. Figure 6.22 presents a power 
spectrum of a new signal captured rolling the nut over the area with a 
simulated damage (pitting) on the leading screw. The maximum values 
from this last spectrum correspond to a frequency of 234.4 and 254 Hz, 
respectively. 

For the measured rotational speed of the leading screw of 960 
rot/min, the rotation frequency is f n = 16 Hz, which leads to f s = 
237.95 Hz and ft = 253.63 Hz, respectively, values of the characteris- 
tic frequencies. It is noticeable that even the power spectrum signals 
the occurrence of a damage by the significant increases in amplitude at 
characteristic frequencies; this time, it is not possible to locate the dam- 
age on one of the rolling paths. This behavior can be due to the play 
between the mechanism’s elements. 

In conclusion, the method of interpretation of the power spectrum 
is an efficient and rapid method of identification of mechanical defects — 
although their localization is, in some cases, more difficult. These reasons 


Table 6.2 Values of the Constructive Parameters that Enter in the 
Calculus of Characteristic Frequencies of the Screw-Nut with 
Balls Mechanism 

Constructive parameters Characteristic frequencies (Hz) 

Dm D b [3 ft f s f p f b 

30 3.00 35° 20° 14.872 f n 15.825 f n 3.355 f n 
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Figure 6.21 Power spectrum resulting from vibration signal captured on 
the “clean” zone of the leading screw using an Emco acquisition system and 
Hypersignal software. 
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Figure 6.22 Power spectrum of a new signal captured rolling the nut over 
the area with simulated damage (pitting) on the leading screw. 


led to the elaboration of the spectrum comparison method, which has 
shown excellent results in profoundness diagnosis. The theoretical pre- 
sentation of this method can be found in Chapter 3 of this book. 

6.6 EXPERIMENTAL RESEARCH IN VIRTUAL 
INSTRUMENTATION 

6.6.1 Structure of Data Acquisition and 
Process System in Virtual 
Instrumentation 

A data acquisition system has to be able to accomplish three funda- 
mental functions: conversion of physical signal in a signal that can be 
measured; measurement of the signals generated by transducers/sensors 
to extract information; and construction of data analysis and data pre- 
sentation in a useful form. A data acquisition and process system uses 
a computer (PC) as a process controller. Figure 6.23 presents the gen- 
eral structure of such a system, and this structure can be described as: 
transducers for conversion of the measured signal into a signal of elec- 
tric nature; circuits to adapt the signal for isolation, conversion, and/or 
amplification of the signal coming from the transducer; a system of mul- 
tiplexors and analogue-digital converters; and a calculus system with an 
adequate software. 

The analogue-digital board (A/D) has the role of transforming the 
analogue signal (a continuous function in time), passed by the adapting 
circuit, in a numerical (discrete) format that can be processed by the 
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Figure 6.23 Data acquisition and process system using a computer (PC) as 
process controller. 



computer. The analogue-digital conversion is a comparison operation; 
the signal is compared to a reference value, which is then represented as 
a codified digital number. There are a minimum and maximum number 
of samples that must be captured in order to optimize the accuracy of 
measurements. In the data acquisition and process system architecture, 
the role of the analogue-digital interface is very important because this 
interface must have a few essential functions for the user: 

High speed data transfer to PC on the DMA (direct memory ac- 
cess) channel 

Electronic part for starting the hardware and software 
Amplification with programmable gain and noise filtration 
Memory buffer 

The operational system through which the virtual instrumentation pre- 
sented in Chapter 5 was elaborated, tested, and then applied, contains: 
an IBM PC 486 computer (100 KHz, 16 Mb RAM), data acquisition 
board AT-MIO-16L-9 National Instruments, conditioning amplifier 2625 
B&K, piezoelectric vibrations transducer 4344 B&K, LabVIEW 2.5 Na- 
tional Instruments software, and a real-time frequency analyzer 2034 
B&K for testing virtual instruments. This system was provided by the 
S.C. Eurotest S.A. Laboratories, Seism-Vibrations department, together 
with their support and collaboration all throughout the experimental re- 
search. 


6.6.2 Experimental Research in 
Surface Diagnostics 

The virtual instruments for surface diagnosis FACVARF and INDDIAG, 
whose structure was presented in Chapter 5, have been used to monitor 
the evolution of mechanical damage on the radial ball bearing 6209. The 
experimental setup is similar to that presented in Figure 6.11, only the 
virtual instrumentation took the place of traditional equipment. 

Figure 6.24, showing diagnosis in virtual instrumentation by the 
peak factor method, presents the front panel of the FACVARF virtual 
apparatus, which exhibits the reference acquisition graphic in the time 
domain, for the new bearing, and under this graphic, another graphic 
of a subsequent recording when the mechanical damage already started. 
In the graphic of peak factor evolution the spectacular decrease of the 
peak factor value from F v = 10.1 to F v =8.2 is due to bearing lubrica- 
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Figure 6.24 Diagnosis by peak factor method on the front panel of the FACVARF virtual 
apparatus. 


tion, which had been mounted without lubricant to hasten the damage 
occurrence. 

The INDDIAG virtual apparatus for the same monitoring and di- 
agnosis operation was used simultaneously (Fig. 6.25). Data acquisitions 
were also made with the bearing in perfect functioning state after the 
damages started, and a diagnostic index evolution graphic was auto- 
matically drawn. Both front panels display a legend to facilitate the 
interpretation of results. 

The Kurtosis method in virtual instrumentation (virtual apparatus 
INDKURT) was used in many diagnostic situations both in the case of 
ball bearings and the screw-nut with balls mechanism, because it is pos- 
sible to estimate without knowing the machine’s “history,” the function- 
ing state of the analyzed element. Figure 6.26 presents this method for 
diagnosis of the same bearing. The value of the Kurtosis index @2 = 3.59 
and the probability density distribution of the amplitude of the captured 
signal confirm that the bearing is in perfect functioning state. 

Other front panels of the INDKURT virtual apparatus are pre- 
sented in the appendices of this book for other diagnostic situations on 
the basis of the vibrator signal. 

6.6.3 Experimental Research in 
Profoundness Diagnostics 

The virtual apparatus MECOSPEC was presented in Chapter 5. This 
apparatus is destined to establish the technical diagnostic using the spec- 
trum comparison method. Given the precision of this apparatus, noted 
during the laboratory tests, this virtual instrument was used during 
experimental research to monitor the functioning state of some needle 
bearings. The purpose of this action was to see if the calculus relations 
presented for determination of critical frequencies for the ball bearings 
were verified also in case of needle bearings. 

Table 6.3 presents the constructive parameters that compete with 
the calculus of characteristic frequencies of bearings, and also the values 
of these frequencies, depending on the rotation frequency f n of the inner 
ring. The characteristic frequencies were calculated using relations (6.12) 
to (6.14). 

Figure 6.27 presents the captured signal in the time domain and 
passed in the frequency domain by a fast Fourier transform for a radial 
needle bearing NA 4909 in perfect functioning state. The resulting power 
spectrum is considered a reference spectrum and is used to obtain the 
warning template and alarm template, respectively. 
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Figure 6.25 Diagnosis on the front panel of the INDDIAG virtual apparatus. 
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Figure 6.26 Kurtosis method in virtual instrumentation (virtual apparatus 
INDKURT). 









Table 6.3 Constructive Parameters that Compete with the Calculus of 
Characteristic Frequencies of Bearings and the Frequency Values 



Constructive parameters 


Characteristic frequencies 

Bearing 
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Ie 
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fru 

type 

Dm 

d 

z 

n 

C fn ) 

(fn) 

(fn) 

NA 4907 

45 

3.5 

27 

0 

12.45 

14.55 

6.39 

NA 4908 

51 

3.5 

29 

0 

13.00 

15.50 

7.25 

NA 4909 

56.5 

4 

25 

0 

11.62 

13.38 

7.03 


The current spectrum is compared to these templates in the next 
section of the virtual diagnosis instrument. This way the warning tem- 
plate obtained by raising by the power two and side plating to the oc- 
tave of the reference spectrum serves as a comparison element for a wear 
spectrum taken after the debut of mechanical damage (Fig. 6.28). The 
difference between the two spectra is graphically displayed to allow a 
direct read. 
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Figure 6.27 Captured signal in time domain and passed in frequency do- 
main by FFT for radial needle bearing NA 4909 in perfect functioning state. 
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Figure 6.28 Warning template obtained by raising by the power two and 
side plating to octave of the reference spectrum serves as comparison element 
for wear spectrum taken after debut of mechanical damage. 


During determinations only a few and the same frequencies have 
increased overpassing the alarm template, a phenomenon which is pre- 
sented in Table 6.4. The characteristic frequencies calculated as a func- 
tion of the input rotational speed n = 960 rpm (/ n = 16 Hz) applied to 
the inner ring are: 

Rolling frequency of needles on the outer ring, f E = 185.92 Hz 
Rolling frequency of needles on the inner ring, fj = 214.80 Hz 
Rolling frequency of needles in the cage, f m = 112.48 Hz 

It can be stated that the relations for the calculus of the characteristic 
frequencies for one-row ball bearings can be applied also to radial needle 
bearings comparing the above calculated frequencies with those exper- 
imentally obtained (Table 6.4) by the intensive wear of the bearing on 
the test stand. 

Figure 6.29 graphically presents the evolution of overpasses to the 
three frequencies mentioned in Table 6.4. It is considered that the evo- 
lution of the frequencies that signal the mechanical damage is approx- 
imately linear in the regime wear domain, so the points resulting from 



Table 6.4 Comparing the Calculated Frequencies with Those 
Experimentally Obtained 


Overpasses (dB) 
at read no. 

7 

9 

11 

13 

15 

17 

Frequency (Hz) 
183 

0.5 

0.9 

2.8 

5.3 

7.7 

9.4 

208 


0.3 

1.2 

2.9 

3.6 

4.8 

230 





1.6 

2.0 


measurements can be interpolated to a line. At any moment, the time 
left until the stop by damage can be estimated by intersecting the inter- 
polated line with the alarm level for the researched frequency. 

In the last version of the virtual apparatus diagnosis MECOSPEC 
(Fig. 6.30) the spectrum comparison and supplementary processing of 
the warning template overpass is replaced with a representation type 
WATERFALL, which is a tridimensional representation made for the 
power spectra coming from successive captured signals at the same mea- 
suring point. At the spectrum axis, frequency and amplitude, the time 
axis is added. The spectra are successively represented on the time axis 
depending on the acquisition moment (Fig. 6.31). On this representation, 
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Figure 6.29 Evolution of overpasses to frequencies mentioned in Table 6.4. 
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Figure 6.30 WATERFALL representation on the last virtual apparatus di- 
agnosis MECOSPEC. 
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Figure 6.31 In the WATERFALL representation, the spectra are succes- 
sively represented on the time axis depending on the acquisition moment. 


TUT 



with the aid of a time cursor any of the represented frequency spectra 
can be selected, and with the aid of a frequency cursor the evolution in 
time can be monitored for any frequency of the vibration signal from 
the power spectra. 


6.6.4 Advantages of the Use of Virtual 

Instrumentation in Technical Diagnosis 

The use of virtual instrumentation in monitoring and establishing of 
technical diagnosis offers special facilities to users as follows. 

Visual Programming. This type of programming is primarily a 
rapid, sure, and easy path to the elaboration of special applications, 
adapted to imposed requests, easy to test and to maintain. It is the type 
of programming recommended to those who have not had the opportu- 
nity to familiarize themselves with a “classical” programming language. 
Instead of writing hundreds of rows of program, the user can build a 
program as a functional diagram choosing from the menu a series of 
visual elements that he or she logically interconnects. The rapidity and 
simplicity of the work stimulate conversion of traditional techniques of 
diagnosis and help find new virtual instrument techniques. 

High Rate of Sampling. The rate of sampling (expressed in sam- 
ples per second or sometimes in Hz) represents the measure of speed 
by which the A/D board can scan the input channel, and can identify 
the discrete value of the signal depending on the reference value. The- 
oretically, a data acquisition system must sample with a speed at least 
twice that of the largest frequency that can exist in the input signal; if 
this rule is not realized, a completely different wave shape will be ob- 
tained, having a smaller frequency, a phenomenon called aliasing. The 
high sample speeds that are reached at present rapidly load the com- 
puter memory but this phenomenon is no longer a problem due to the 
technological development of computers. The diagnostic problems need 
a medium level of sample rate that can be solved with 8 or 16 Mb RAM. 

Good Resolution of Measurements. The resolution (expressed in 
percentages or bits) defines the smallest variation of the input signal that 
can be detected by the acquisition system. In virtual instrumentation 
a resolution comparable to that of dedicated traditional instruments 
produced by prestigious companies can be obtained. 



Possibility of Choice of Sampling Mode. In applications of fre- 
quency analysis using FFT, any deviation in the period of time between 
sampling produces considerable errors. The acquisition system with vir- 
tual instrumentation realizes the start of the A/D converter (triggering) 
directly by the clock from the hardware or an external clock on the acqui- 
sition board. More than that, the possibility of choosing the pretrigger 
or postt rigger modes exists. 

Data Transfer Speed. Using the direct memory access transfer 
mode the system takes over the data from the acquisition interface and 
puts these data directly in the computer's memory. DMA transfers are 
directly controlled by hardware and are extremely rapid. Ultrafast ac- 
quisition systems that use memory located directly on the acquisition 
plate also exist, so these systems will no longer be limited by the com- 
puter thoroughfare speed. In these conditions a “real-time” data process 
becomes possible, an accomplishment hard to achieve and expensive to 
realize with traditional instrumentation. 

Possibility of Using the Multiplexor. The multiplexor is an elec- 
tronic device that disposes more input channels, an output channel, and 
digital control inputs. Using these control inputs the needed input chan- 
nel can be selected and connected to the output channel. The multiplexor 
is useful for supervising machines with a large number of measuring 
points, using the same measurement instrument. 

Use of Work Programs with Annexed Apparatus (Driver-E Soft- 
ware). These work programs with annexed apparatus constitute a soft 
level that directly programs the acquisition hardware, administers the 
functioning of the acquisition hardware, and ensures integration with 
computer resources. They hide the complicated details of hardware pro- 
gramming, ensuring an easy to understand interface for the user. 

Data acquisition systems tend to be more and more mass con- 
sumption goods because of the continuous technological improvement 
of this peak domain. These conditions determine the transformation 
software, an important factor related to data acquisition and process 
differentiation. 
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A Neural Approach to Establishing 
Technical Diagnosis for 
Machine Tools 


7 . 1 INTRODUCTION TO NEURAL 
NETWORK THEORY 

The end of the 1980s signaled a period of research regarding the de- 
velopment of a new approach to data processing within computational 
systems. This approach earned several names: connectionism, neural net- 
work, neural processing, and parallel processing; all of these terms are 
synonyms for that kind of information processing that tries to simulate 
the thought processes of the human brain based on experience, rather 
than the classical method of that based on algorithms. 

7.1.1 Information Neural Processing 

Data neural processing means to elaborate and study networks with 
adaptable nodes; these networks store experiential knowledge, and are 
available to be utilized for the purpose for which the network was cre- 
ated. The word “adaptable” used for the network nodes represents that 
property of a node which is able to react correctly even if the stimulus 
received is not precisely the same as the ones already learned. 
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Figure 7.1 Adaptable node, main element of a neural network. 


The adaptable node, the main element of a neural network, can 
be represented as an electronic device with a manifold input channel 
(Fig. 7.1). It consists of one output channel and two special connections: 
one for introducing data into the learning process, and the other as a 
switch of learning/utilizing nodes. The F node function is the way in 
which an output parameter is associated with each input dataset. For 
example, the function of an adaptive three-input node can look like a 
logic table in which the response is 1 only when the input dataset has 
a single 1 value (Table 7.1). This representation of an adaptive node 
makes obvious the association with the neuron, the adaptive node of the 
brain. The human brain has an average of 10 11 neurons, organized in 
complex structures. 

The neuron (Fig. 7.2) consists of a cellular body called the peri- 
karion, which contains the nucleus, and two sorts of protoplasmic pro- 
longation: the axon (cylindrical shaped, long, usually single) and the 
dendrites (many and short). The dendrites represent the input channel 


Table 7.1 The Function of an Adaptive Three-Input 
Node Can Look Like a Logic Table 

xl 0 0 0 0 1 1 1 1 

x2 0 0 1 1 0 0 1 1 

x3 0 1 0 1 0 1 0 1 

Y 0 110 10 0 0 



Figure 7.2 Neuron schema. 

into the neuron, while the axon is the output channel. An electric activity 
characterized by short and fast impulses (about 100 impulses/second) 
has been noticed in the axon when the neuron “emits.” Neurons are 
interconnected through the ends of the dendrites, called synapses. One 
neuron can receive 5000 to 15,000 input signals from other neurons’ 
axons. Synapses can be excitators, if they help the neuron to emit, or 
inhibitors, if they discourage the neuron from emitting. 

The informational model of a neuron was first proposed in 1943 by 
the neurophysiologist Warren McCulloch and the logician Walter Pitts. 
This model is still the basis of information neural processing. According 
to this model, the synaptic modifications are continuous, and the neuron 
takes into consideration all of the synaptic signals, both excitator and 
inhibitor; then, it sums up their effects and determines whether to emit 
them through the axon. Number 1 is associated with the emission state, 
and 0 with the repose state. 




If X denotes the neuron state and the synaptic connection effect 
is represented by a weight W, then the effect of a synapse upon the 
neuron is given by the product X • W. The weight W can have values 
within the interval — the negative values characterize the in- 
hibitor synapses. Since there is a multitude of synapses, an index j has 
to be attached, which designates Xj and Wj as the input and the weight 
of the synapse j, respectively. The neural model constantly adds those 
effects and compares them with a threshold value T ; if the sum exceeds 
the threshold, the neuron emits. For the McCullogh and Pitts neuron, 
this emission rule can be described mathematically by the relation: 

X ± W ± + X 2 W 2 + • • • + XjWj + • • • + X n W n > T. (7.1) 

The electronic interpretation of this relational model is given in 
Figure 7.3. The basic components are: the summator amplifier, which 
provides an output voltage proportional to the sum of all products X-W, 
and the voltage comparator, which generates a voltage equal to 1 if 
the output voltage of the summator exceeds the threshold voltage T. 
The value of the adjustable weights W is set automatically during the 
learning process. 

On the basis of the model presented so far, two sorts of neural net- 
works have been theorized: feedforward and feedback networks. Within 
these networks neurons are distributed in one or more layers that cannot 
be accessed directly; they are called hidden layers and have free access 
only at the input channel and, respectively, at the output channel of 
the network (Figure 7.4). The feedforward networks operate between 
the input channel and output terminals, learning to associate output 
data with the input data. In feedback networks, the information can 
run within a loop from the input to the output channel and vice versa, 
creating a so-called internal input channel. Both types of network oper- 



FlGURE 7.3 Electronic interpretation of Eq. (7.1) relational model. 




Figure 7.4 Neurons distributed in one or more layers that cannot be accessed 
directly. 


ate under emission rules that determine whether the neuron will emit 
for each dataset. These rules are associated with each connection in the 
network, and they determine the weight of each of those connections. 
As a consequence of the application of the emission rules, the network 
becomes able to generalize, that is, to provide an adequate response for 
input datasets that were not introduced during the learning process. 

The learning algorithms, also called learning laws, are based on the 
following condition: when the weights are recalculated for a new input 
dataset, no discontinuity with what was done for previous sets should 
occur. The biological basis of these laws derives from the hypothesis 
that, if a neuron is part of a network in which a synaptic input emits 
continuously at the same time as the neuron, then the weight of that 
synaptic connection will increase. The disadvantage of this theory is 
that it leaves the responsibility of the neuron’s emission to the action 
(eventually learned) of some synapses. Although the technique proposed 
by Hebb was useful for the investigation of the neural network, another 
rule proved more suitable: the delta rule, which can be described through 
the following steps: 

1. Select an input dataset. 

2. When an error is detected in the network’s response, calculate 
the deviation from the desired output parameter. 

3. Adjust the active weights (i.e., those which are emitting) and 
the threshold value, in order to partially correct the error. 

4. Return to step 1 until no input dataset causes errors. 
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Figure 7.5 Perceptron electronic model: Aj are the preprocessing units, 
called associative units. 


The McCulloch and Pitts model was improved by Frank Rosenblatt 
by adding fixed preprocessing units; their duty is to extract specific 
features from the input signals. Thus the perceptron was born, which is 
defined as a model recognition device. Figure 7.5 shows the perceptron 
electronic model; Aj are the preprocessing units, called associative units. 
In terms of these units, the perceptron order is defined; the preceptron 
order is equal to the number of inputs of the associative unit with the 
largest number of inputs. 

7.1.2 The Learning Process in 
Neural Networks 

The information neural processing paradigm principle states that neu- 
ral networks can be trained merely through examples from the exterior. 
In other words, any learning algorithm utilized in multilayer networks 
is based on output error estimation. This is the so-called hard learn- 
ing problem, and it is one of the most important problems in neural 
networks. 

A major step in overcoming this problem was made in 1982 by 
John Hopfield, an American biologist and chemist; he presented neu- 
ral networks as content addressable memories (associative memories). 
Through his work, he made two major contributions. He developed a 
type of network analysis that uses the concept of energy, concluding 
that a network reaches, while operating, its energetic minimum, after 





which the output signal set does not change (i.e., stability occurs). In 
addition, he showed that learning rules, such as the delta rule, can be 
utilized to adjust the network parameters purposely to create the ener- 
getic minimum. The Hopfield neuron model has the characteristic pa- 
rameters Vj, the output signal (V = 0 if the neuron does not emit; 
V = 1 if the neuron emits); the weight of the connection of the 
neuron i with the neuron j (T = 0 if the neuron i is not connected to 
the neuron j ); and £//, the threshold value for which the neuron emits. 
In short: 

Vi becomes 1 if Tjj Vj > Ui 

Vi becomes 0 if Tjj Vj < Ui 

The energy of a neuron can be calculated with 

Ei = -Vi - Ui) (7-3) 

The amount between brackets is called the activation energy and is de- 
noted by Ai. In order to have a stable state of a neural network, none of 
the nodes should be activated in such a way that the emission conditions 
change. So, if a neuron is emitting (V = 1), its activation will be posi- 
tive, so that the emission is not further stopped; the same is true if the 
node is not emitting ( V ^ 1); its activation will be negative, so that it 
does not further emit. Consequently, when the neural network changes 
its state, either it keeps the same energetic level, or it goes down to a 
lower energetic level. When lower energy levels are no longer accessible, 
the network remains stable at its most recent state. 

From an energetic analysis perspective, it becomes clear that a 
learning procedure application, such as the delta rule, is simply a way 
to decrease the state energy to a minimum. The Hopfield model has 
its disadvantages: the most important is that the neural network “gets 
stuck” in false energetic minimums. The way that a neural system runs 
out of false minimums was discovered by Geoffrey Hinton; it consists 
of “noise” utilization (i.e., applying an uncertainty degree to the state 
energy). Intuitively, this method can be illustrated by representing the 
network state as a ball on a waved surface (Fig. 7.6). If the ball has 
an internal property which makes it “jump,” it is probable that it will 
spend the longest time in the deepest valley it can reach. 




Figure 7.6 Representing the network state as a ball on a waved surface. 


At the crux of any neural network activation remains the phe- 
nomenon of increasing the aleatory motion of gaseous molecules while 
temperature increases, discovered by Ludwig Boltzmann at the end of 
the nineteenth century. Analogously with Boltzmann’s research, the un- 
certainty degree, introduced to assess a neural network state, was named 
temperature. Thus, at zero temperature the network will behave accord- 
ing to the Hopfield model; at higher temperatures, an uncertainty degree 
proportional to the temperature has to be introduced in the neuron acti- 
vation function. This procedure has the advantage of helping the network 
to run out of its false minimums but there is another side of the coin: the 
network no longer remains stable. Hinton proposed a “thermal regime” 
be utilized to enable this trouble to be outrun: starting the network at 
a high temperature and then cooling it gradually until it reaches the 
“thermal equilibrium.” In this way, the network has the greatest chance 
to end in a state associated with the lowest minimum to be reached for 
specific input data. This manner of approach to neural networks has 
been called “simulated tempering.” 

The introduction of temperature in the neuron activation function 
is done through a probabilistic function, called the emission probability 
Boltzmann function: 

1 

1 + exp (-A/T) 

The plot of this function is represented in Figure 7.7 for “temperatures” 
of 0.5 and 0.25 (these are arbitrary units, not related either to Celsius 
or Kelvin degrees or to other measure units of temperature as a physical 
parameter). The result of the introduction of temperature as a mod- 
erator in the neuron activation function can be pursued in changes of 
the way in which the electron’s emission is interpreted. According to 
the Hopfield model, if the A activation were negative, the neuron would 
not emit \p( 1) = 0]; if the activation were positive, the neuron would 




I 



have emitted \p( 1) = 1]; this is shown in Figure 7.8a. The new model, 
named the Boltzmann machine, suggests that, for negative values of the 
A activation, there is a probability p( 1) = 1 that the neuron emits; this 
probability goes to zero as A becomes negatively smaller and smaller. 
Similarly, there is a probability _p(0) that the neuron does not emit, 
even if the activation is positive; this probability gets lower when the 
activation gets higher (Fig. 7.8b). 




Figure 7.8 (a) Hopfield model; (b) new model, named the Boltzmann 

machine. 


The Boltzmann machine showed that the probability that the net- 
work comes to an end in a certain state depends nearly exclusively on 
that state energy. That means that, through the learning process, the 
energies corresponding to the minimums of the system have to be well 
controlled. 

Starting from the assumption that not all units of a network are 
defined by the learning dataset, Hinton showed the necessity of hidden 
units, essential for solving the hard learning problem. The learning pro- 
cedure, based merely on information as to whether the visible units be- 
have correctly, has to ensure that hidden units develop their own correct 
weights and thresholds, too. If we assume a neural network has v visible 
units, it can have 2 /y possible states. If Si, S 2 , . . . , S r are input datasets 
for training the network through visible units, state probabilities can be 
calculated: P + (Si), P + (S 2 ),..., P + (S r ), where the “+” sign is used 
to indicate that these are desired probabilities. On the other hand, the 
“— ” sign is used for the same probabilities, but results from the free 
(untrained) run of the network, that is, P~(S 1 ), P~(S 2 ), ■ . • , P~(S r ). 
The function introduced to measure “the distance” between these two 
probability sets is: 

G = E p+ (^) ln ^|f y ( 7 - 5 ) 


The first term of the previous relation, P + (5 a ), makes the states with 
the largest occurrence probability have the largest effect upon the sum; 
the natural logarithm becomes zero when the two probability sets are 
identical [P + (5 a ) = P~(S a )]- 

It can be shown that the “distance” changing rate between the two 
probability sets depends on the trained, respectively, untrained, visible 
units 5 temperature and average state probabilities, after the relation: 




The meaning of this relation consists of the following: in order to lower 
the rate G by changing the weight uyp all that needs to be known is the 
local information (p^ — p ^, • ). If this term is positive, the weight Wij has 
to be increased; in the opposite case, it has to be decreased. Obviously, 
the process is finished when G — 0, so the network learned to reproduce 
the state probabilities, and it is considered completely trained. 



As shown before, there are two types of neural networks: feedfor- 
ward and feedback. For feedback networks it is simple to introduce into 
the loop the signal associated with the output error, and this is done 
on purpose to diminish this error while training. As far as feedforward 
networks are concerned, they are the most utilized ones, and they can- 
not work freely because of the lack of reaction to inform; therefore, a 
training method for the hidden units is necessary, through a process of 
propagation of the measured error on a backward direction, from the 
output channel. 

The hidden layers can be considered the location in the network in 
which input parameters are partially processed and labeled before the 
final result is reached in the output layer. In these layers, representa- 
tions are formed, which are not provided during the training process. 
A generalization of the delta rule forms the basis of the hidden units’ 
behavior through the converse error propagation process. The goal of 
this method is to minimize the overall output error £_p, defined as the 
half-sum of the squares of all neurons’ output errors: 





where t p j is the target output of the neuron j for the input dataset p ; 
o P j is the actual output of the neuron j for that input dataset. 

The training process through converse error propagation has two 
steps. The first is the forward step, during which the input data are 
applied and let run to the output channel. The output parameters are 
calculated and compared to the target parameters (which have to be 
known). During the backward step, the errors (those resulting as a con- 
sequence of the comparison in step one) are propagated backwards to the 
input layer. The purpose is to recalculate the neural connection weights. 
Another forward step follows, then another backward step, until the error 
minimizes within a preset limit. For this manner of feedforward neural 
network training, it has been noted that the most suitable neuron activa- 
tion function is the same type as the function presented in relation (7.4) 
and Figure 7.7, also called the sigmoid function. A feedforward type 
neural network with three layers (input, hidden, and output), having 
a sigmoid activation function, and being trained through the converse 
error propagation method, can be represented schematically as in Fig- 
ure 7.9. It should be noted that a larger number of neurons in the hidden 
layer can guarantee a better result in network training and usage. 




Figure 7.10 Structure of Windows Neural Network. 








7.2 UTILIZATION OF NEURAL NETWORKS IN 
MACHINE TOOLS DIAGNOSTICS 

7.2.1 Multicriteria Application for 

“Good/Defective” Classification 

The problem of identifying defective elements in the kinematic structure 
of a machine tool can be solved by utilizing feedforward type neural net- 
works. The following deals with elaboration and training of a neural 
network for bearing diagnosis. The data used for network training come 
from the experimental research presented in Chapter 5. Windows Neu- 
ral Network is a user interface, written in Visual Basic, for building 
and utilizing feedforward neural networks, fully connected and trained 
through the converse error propagation algorithm. The structure of the 
neural networks with which this application works (Fig. 7.10) contains 
the input layer, one or more hidden layers, and the output layer. The 
output channel of each neuron is connected to all neural input channels 
of the next layer, which can introduce BIAS units (tendency or predic- 
tion units), which facilitate the network's training. There is no activation 
function in the input layer; this layer can only distribute the input data 
to the first hidden layer. The hidden layers and the output layer have 
an activation function, which can be: 

A linear function: 

f(x,T)=xT (7.8) 

A hyperbolic tangent function: 

f(x,T) = th(xT) (7.9) 

A sigmoid function: 

f(x,T) = l[l + exp(— sT)] (7.10) 

From a mathematical point of view, the network runs the following 
algorithm. 

For the hidden layer output: 
h (j) = ^w{i,j)xi{i) 

i 

s{i) = f(h{j)) 


(7.11) 


For the output layer output: 
h'(k) = ^ w(j, k)xs(j) 

i (7.12) 

o(fc) = f(ti(h)) 

where i(i) represents the network input channels; o(k), the network out- 
put channels; w(i,j), the weight of the connection of the neuron I with 
the neuron j in the next layer; and /, the activation function. 

The objective function has to minimize the final error: 

Pmax ^'max 

SRMS =^2^2 fttP’ k )] 2 ( 7 - 13 ) 

p— I k—I 

where t(p, k) is the output target value for the output dataset p, and 
o(p, k) is the actual output value for the same dataset. 

The DIAGNO neural network for identifying defective bearings has 
the structure shown in Figure 7.11. The network's dimension is given by 
the number of layers and the number of neurons in each layer; conse- 
quently, the built network is a 3 x 4 x 1 type, with BIAS units in 
layers 2 and 3. The neuron activation function is the sigmoid function 

(Fig. 7.10). 
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Figure 7.11 DIAGNO neural network for identifying defective bearings. 


The selection criteria for data needed to train the network are the 
peak factor criterion, the diagnosis index, and the Kurtosis index. Data 
corresponding to the evaluation through these criteria come from the 
database of vibration signals stored and processed in the time field (see 
Chapter 5); they are provided to the network in three ways. The train- 
ing datasets, containing input data and the corresponding output target 
parameter, are shown in Table 7.2. These data cannot be introduced 
in the network as they are, since they could not be processed by the 
activation functions. For instance, the sigmoid function can only take 
values between —2 and 2 at the input channel; other values would sat- 
urate the neuron and lead to an output with value 0 or 1 all the time. 
To avoid this inconvenience, the input data should be logarithmically 
or linearly normalized, as the variation field of those data is larger or 
narrower. This operation is done automatically, on demand, from the 
main menu of the application; the normalization can be global (to all of 
the network’s nodes) or individual (node by node). 

The next step is represented by the network’s converse propagation 
algorithm adjustment, through the two parameters that connect the new 
weights with the derivatives of the old weights, after the relations: 

1 ) = rjdW ( i,j , t) + adW (i, j, t - 1 ) 

W(i, j, t + 1) = W{i, j, t) + dW (■ i,j , t) 

where r/ is the learning parameter and a the moment parameter. The 
best results were obtained for 77 = 0.2 and a = 0.5, as the front panel of 
this application shows (Fig. 7.12). It should be mentioned that during 
the network run there was no need to introduce “noise” and “tempera- 
ture” to help the network avoid getting stuck in local minimums. Since 
temperature is a multiplicator of the activation function argument, not 
to take it into consideration means T = 1. 

The structural network training process runs in epochs. An epoch 
represents the number of input datasets in terms of which weights are 
calculated. Thus, if the length of an epoch is 1, a continuous training 
is done (weights are recalculated after each dataset); if the length of an 
epoch is equal to the number of datasets, a simultaneous training is done 
(weights are recalculated after each pass over all datasets) . Simultaneous 
training is useful when there are few input datasets, as in the present 
case. At the same time, this sort of training is faster than continuous 
training. Training this neural network needed 555 iterations, that is, 
555 passes over the input datasets. Training ended when the target was 


(7.14) 



Table 7.2 Training Datasets, Containing Input Data and Corresponding Output Target Parameter 

Set number 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Peak factor 

10.0 

7.5 

20.0 

27.0 

6.0 

25.5 

26.0 

29.0 

14.0 

24.5 

Diagnosis index 

.500 

.650 

.170 

.013 

.800 

.018 

.015 

.010 

.400 

.019 

Kurtosis index 

4.1 

3.7 

7.2 

9.2 

3.4 

9 

9.3 

12.0 

5.8 

8.8 

Target output 

1 

1 

0 

0 

1 

0 

0 

0 

1 

0 
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Figure 7.12 Front panel caption. 




Table 7.3 Weights Wij Calculated for Each Neural Connection 


Total calculated weights: 21 

Between input layer and second hidden layer: 



2.567564 

-0.552004 

-0.623598 

0.690139 


1.273501 

0.300688 

1.931408 

0.233710 


-3.382554 

1.183681 

0.871934 

1.333970 


2.154515 

-0.137544 

-0.532012 

1.897794 


Between second hidden layer and third output layer: 



-2.092981 

-2.371565 

4.373098 

1.461245 

-0.007646 


reached with an error smaller than the preset one, in this case. It can 
be noticed in the window with the training results (Fig. 7.12) that the 
actual root mean square error is £rms = 0.000379604. 

As a consequence of training the network with experimental data, 
weights have been calculated for each neural connection (Table 7.3), 
through the converse error propagation algorithm. Weight distribution 
in the network is shown in Figure 7.13 in the shape of a 10-interval his- 
togram. The number of weights in each interval is represented in terms 


Weight Distribution 


Count 




Weight Value 


Figure 7.13 Weight distribution in the network. 




of the weights' value, between —4 and +4. This representation provides 
information about the quality of the result obtained in the calculus of 
weights: normally, for a network with the same activation function in all 
layers, the weight distribution should be as close as possible to a Gauss 
distribution. Figure 7.13 shows the weight distribution for a DIAGNO 
neural network, correlated with data in Table 7.3. The shape of this his- 
togram gets closer to a Gauss distribution, so the weights are correct. 
In these conditions, a graphical representation of the actual network 
outputs can be plotted at the same time as the target outputs. In Fig- 
ure 7.14 the network’s target outputs, as they were set (Table 7.2), are 
drawn with a solid line, and the trained network’s actual outputs with 
a dashed line. 

The error between the target output and the actual output can 
also be plotted, as in Figure 7.15. It can be clearly seen that the neu- 
ral network observes some trouble in interpreting datasets 2 and 8 (see 
Table 7.2). Indeed, during the training process, the two input datasets 
had been introduced artificially by the author within categories “defec- 
tive” and “good.” 

At this time, the DIAGNO neural network training process can be 
considered finished and the network can be utilized in the classification. 
Tests of the built neural network have been carried out on a group of four 
datasets known from experimental research; these datasets are presented 
in Table 7.4. It should be noted that this time the input target value 
for each input dataset was not indicated; the neural network will have 
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Figure 7.14 Network’s target outputs, as they were set (Table 7.2) (solid 
line) and trained network’s actual outputs (dashed line). 



Figure 7.15 Error between target and actual output. 


/ Error 


to perform the set classification, in the manner it learned during the 
training process. The weights of neural connections calculated during 
the training process were used as they were during the test. After the 
normalizing operation, the test data were plugged into the network. The 
training parameters were kept the same (77 = 0.2 si a = 0.5); it was not 
necessary to use noise or temperature to avoid the network getting stuck 
in false minimums. Data for running the classification test for the four 
bearing types from which input data were taken are shown in Table 7.5. 
The root mean square error in the test was smaller than the one that 
resulted in training the network, as was expected. 

Figure 7.16 illustrates the result of the neural network test: indeed, 
the first two test datasets belonged to two radial ball bearings 6209 which 
were working perfectly; the last two datasets were taken from the same 
bearings, having pitting on the outer race track due to an intensive wear 
on the stand described in Chapter 5. 


Table 7.4 Test of the Built Neural Network Carried Out on a 
Group of Four Datasets Known from Experimental Research 


Set number 

0 

1 

2 

3 

Peak factor 

6.50 

11.00 

26.70 

25.40 

Diagnosis index 

0.780 

0.480 

0.013 

0.017 

Kurtosis factor 

3.50 

4.50 

9.70 

9.10 


Table 7.5 Data for Running the Classification Test for Four Bearing 
Types from Which Input Data Were Taken 



RMS error: 0.000103 
Good pats: 100.0% 


Set 

number 

Output actual value 

Output adopted value 

Error 

0001 

0.994650 

1.000000 

0.000029 + 

0002 

0.993171 

1.000000 

0.000047 + 

0003 

0.012820 

0.000000 

0.000164 + 

0004 

0.013081 

0.000000 

0.000171 + 


In conclusion, the DIAGNO neural network is able, after a pre- 
liminary learning process, to distinguish perfectly working bearings and 
defective bearings with remarkable accuracy. Other tests carried out 
with the same neural network, over other experimental datasets, gave 
the same results every time. 

7.2.2 Application for Neural Diagnosis of the 
Working State 

The goal of a diagnostic operation is to detect “ante factum” a defect, 
in order to prevent a machine’s being out of operation due to damage. 
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Figure 7.16 Result of neural network test. 







Starting from this consideration, a neural network utilized for monitoring 
and diagnosing should be able to evaluate the real working state of the 
supervised element. Therefore, a neural network named DIAGNOZA has 
been built, having the structure shown in Figure 7.17; this network has 
been trained to classify supervised bearings into four categories: 

Bearings in a perfect working state (target y = 1) 

Bearings that have conditions for defects to occur (target y = 2) 
Bearings in a limited working state (target y = 3) 

Defective bearings (target y = 4) 

The network is of the multilayer feedforward type, size 2x5x1, with 
BIAS type prediction units, able to be trained through the converse error 
propagation algorithm. The neuron activation function is a sigmoid one 
[relation (7.10)], the same for neurons in all layers. The network is fully 
connected, meaning each neural output channel from a layer is connected 
to all neural input channels in the next layer. 

The two network input channels were given data from the experi- 
mental research and interpolations of those data (Table 7.6); the classi- 
fication criteria considered this time were the peak factor criterion (Fv) 






Figure 7.17 Structure of DIAGNOZA neural network. 



Table 7.6 Data from Experimental Research and Interpolations of Those Data 


Set no. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Peak factor 

25.25 

14.50 

7.36 

26.64 

24.21 

6.64 

26.35 

22.89 

17.12 

Kurtosis index 

9.79 

5.12 

3.61 

11.23 

8.73 

3.47 

10.68 

8.13 

6.25 

Target 

4 

2 

1 

4 

3 

1 

4 

3 

2 

Set no. 

10 

11 

12 

13 

14 

15 

16 

17 

18 

Peak factor 

19.07 

23.64 

27.01 

8.71 

5.72 

8.25 

20.85 

25.98 

12.14 

Kurtosis index 

6.70 

8.46 

11.88 

4.09 

3.06 

3.88 

7.45 

10.19 

4.88 

Target 

2 

3 

4 

1 

1 

1 

3 

4 

2 

Set no. 

19 

20 

21 

22 

23 

24 

25 

26 


Peak factor 

15.63 

26.83 

6.12 

6.00 

14.07 

20.03 

25.55 

21.79 


Kurtosis index 

5.73 

11.57 

3.24 

3.14 

4.81 

7.22 

9.59 

7.89 


Target 

2 

4 

1 

1 

2 

3 

4 

3 






and Kurtosis index criterion (^ 2 ). Data were linearly normalized in or- 
der to be accepted by the neuron activation function. The best results 
were obtained for the values of the learning parameters 7 / = 0.2 and 
a = 0.5, without using noise or temperature as aiding factors in the 
network. On the front panel of this application (see Fig. 7.18), the re- 
sults of the network training can be read, after using the converse error 
propagation method. Thus, after a relatively large number of iterations 
(6140), the working root mean square error decreased under the preset 
value e = 0.001 for all of the input datasets. 

As a consequence of the training process, 21 weights of the neural 
connections were evaluated: 15 connections between layers 1 and 2, and 
6 connections between layers 2 and 3, also taking into consideration the 
connections of the BIAS prediction units in layers 2 and 3. These weights 
are presented synthetically in Table 7.7. The histogram of this weight 
distribution is shown in Figure 7.19. The envelope of these histograms 
is close to a Gauss distribution, which confirms the weight calculus was 
correct. 

The neural network running on the input data was performed using 
the simultaneous training method (the length of an epoch is equal to the 
number of datasets), so the weights were recalculated after each pass 
over all 26 datasets, and this allowed the training time to be shortened. 
Table 7.8 presents the results of the run: the output actual values, the 
target output values, and the interpreting error for each training dataset. 
A graphical plot of the actual and target outputs is shown in Figure 7.20; 
Figure 7.21 shows the output error histogram of the trained network. It 
should be noted that the neural network had classification trouble with 
sets 5 and 25, where the concordance between the evaluation criteria 
was not too good. 

Training of the DIAGNOZA neural network can be considered suc- 
cessfully finished; DIAGNOZA may now be utilized in the problems of 
diagnostic classification for which it was built. The test for the DIAG- 
NOZA neural network was performed on a group of six experimental 
datasets (Table 7.9), coming up after monitoring the forced wear of the 
same radial ball bearings type 6209 (see Chapter 5). 

After data were linearly normalized, they ran in a trained neural 
network. The output root mean square error was 0.000307, therefore 
under the limit of 0.001 imposed for all sets. The network identified 
accurately the working state of the tested bearings; furthermore, during 
the test a false target was indicated for one of the input datasets, but the 
classification performed by the network did not change. This validated 
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Figure 7.18 Results of network training can be read after using converse error prop 
agation method. 










Table 7.7 Twenty- One Weights of the Neural Connections Evaluated 


Total evaluated weights: 21 


Between layer 1 and layer 2 


1.754083 

4.748973 

3.361797 

1.547162 

-3.090095 


-1.154256 

5.315055 

-1.410271 

-6.397412 

0.540312 


-0.962436 

-2.021535 

3.961644 

4.818911 

-5.167025 


Between layer 2 and layer 3 

-2.278813 1.963403 5.296241 -4.803792 -3.292156 0.073231 


that a good training of the network was achieved. Figure 7.22 represents 
the actual values (the white bar) and the adopted values (the black bar) 
of the network output, in the test case. 

In conclusion, training of DIAGNO and DIAGNOZA neural net- 
works with data obtained from experimental research led to excellent 
performance of recognition and classification of elements in the machine 
tool structure. The tests performed — the most representative ones have 
been presented here — confirm the capabilities of this type of information 
processing. 
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Figure 7.19 Histogram of weight distribution presented in Table 7.7. 


Table 7.8 Output Actual Values, Target Output Values, and Interpreting Error for Each Training Dataset 


RMS error: 0.000150 


Good pats: 100.0% 

Set 

Actual 

value 

Target 

value 

Error 

Set 

Actual 

value 

Target 

value 

Error 

0001) 

3.966598 

4.000000 

0.000124+ 

0002) 

2.027430 

2.000000 

0.000084+ 

0003) 

1.008047 

1.000000 

0.000007+ 

0004) 

3.988856 

4.000000 

0.000014+ 

0005) 

3.087048 

3.000000 

0.000842+ 

0006) 

1.003739 

1.000000 

0.000002+ 

0007) 

3.987121 

4.000000 

0.000018+ 

0008) 

2.956078 

3.000000 

0.000214+ 

0009) 

1.933849 

2.000000 

0.000486+ 

0010) 

2.038465 

2.000000 

0.000164+ 

0011) 

2.988344 

3.000000 

0.000015+ 

0012) 

3.990406 

4.000000 

0.000010+ 

0013) 

1.042064 

1.000000 

0.000197+ 

0014) 

1.002214 

1.000000 

0.000001+ 

0015) 

1.024447 

1.000000 

0.000066+ 

0016) 

3.014108 

3.000000 

0.000022+ 

0017) 

3.982717 

4.000000 

0.000033+ 

0018) 

1.977742 

2.000000 

0.000055+ 

0019) 

2.026888 

2.000000 

0.000080+ 

0020) 

3.989686 

4.000000 

0.000012+ 

0021) 

1.002707 

1.000000 

0.000001+ 

0022) 

1.002615 

1.000000 

0.000001+ 

0023) 

2.003102 

2.000000 

0.000001+ 

0024) 

2.939797 

3.000000 

0.000403+ 

0025) 

3.905148 

4.000000 

0.001000+ 

0026) 

3.020487 

3.000000 

0.000047+ 
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Figure 7.20 Graphical plot of actual and target outputs 
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Figure 7.21 Output error histogram of trained network 


Table 7.9 Test for DIAGNOZA Neural Network Performed on 
6 Experimental Datasets 

Set No. 1 2 3 4 5 6 


Peak factor 8.53 23.14 16.05 25.73 18.56 23.9 

Kurtosis index 3.97 8.29 5.48 9.81 6.50 8.6 





| Target-1 

II Net-1 


Figure 7.22 Actual values (white bar) and adopted values (black bar) of 
network output in the bearings test. 


7.3 FINAL REMARKS AND PROSPECTS FOR 
UTILIZING NEURAL NETWORKS FOR 
MACHINE TOOL DIAGNOSIS 

As noted, utilization of neural networks allowed the exploitation at a 
higher level of the data library coming from experimental work and 
running. The success rate of the neural diagnosis depends mainly on the 
following factors. 

1. Neural network architecture: The number of the network’s hid- 
den layers can be neither too high nor too low. A network without hidden 
layers becomes a simple linear separator and does not reach its goal; on 
the other hand, a too large number of hidden layers makes the learning 
process inefficient. Usually, the number of hidden layers is 1 or 2. Dur- 
ing the DIAGNOZA neural network elaboration, the same input datasets 
and the same initial aleatory weights were used for two networks: type 
2x5x1 and type 2 x 4 x 4 x 1. It was noticed that — although both 
networks had a 100% success rate — the learning process took longer in 
the case of the network with two hidden layers (the rapidity of conver- 
gence cannot be a disadvantage, since the learning process is offline). 
The classification errors were nearly one size order larger for more than 
half of the training datasets. 


2. The training dataset size: It is well known and experimentally 
proven that neural network performance increases when the training 
dataset is larger. It is recommended that, when evaluating the necessary 
number of datasets for training, the complexity of the problem to be 
learned by the neural network be taken into consideration. 

Information processing by means of neural networks proved to be 
a viable alternative for classical monitoring and diagnostic techniques. 
A comparison between the success rates of these techniques and those 
based on neural processing (Table 7.10) proves that the latter are better. 

As noted, the neural networks elaborated were trained with data 
that had already undergone a previous process within criterial estima- 
tions (the peak factor criterion, the diagnosis index criterion, or Kurtosis 
criterion). But these networks are also capable of learning from unpro- 
cessed data, provided the target is indicated correctly. Such an approach 
could be done by building a neural network that can be trained by means 
of the captured signal; this network can be represented by a power spec- 
trum. Figure 7.23 illustrates this spectacular change from neural net- 
works fed with structured data (Fig. 7.23a) to a neural network able 
to recognize defects through a global representation of vibration signals 
(Fig. 7.23b). 

In previous chapters, an analysis of characteristic frequencies was 
described for some elements of a kinematic chain structure. It was shown 
that, within the power spectrum obtained by applying the fast Fourier 
transform to the signal captured in the time field, the increase of ampli- 
tude for these frequencies is the result of certain typical defects’ nucle- 


Table 7.10 Comparison Between the Success Rates of 
Classical Monitoring and Diagnosing Techniques and Those 
Based on Neural Processing 


No. 

Monitoring/ diagnosing technique 

Success rate (%) 

1 . 

Peak factor method 

50-70 

2. 

Diagnosis index method 

62-70 

3. 

Kurtosis method 

67-75 

4. 

Spectrum comparison method 

70-85 

5. 

Envelope method 

75-85 

6. 

Diagnosis neural networks 

95-100 

7. 

Classification neural networks 

98-100 



Figure 7.23 (a) Neural networks fed with structured data; (b) neural net- 

work able to recognize defects through global representation of vibration 
signals. 


ation. The nature and proportion of these defects can be identified by 
supervising the characteristic frequencies. 

It is possible to build a diagnostic and classification system for 
defects on the basis of power spectrum recognition, by means of a prin- 
ciple schema similar to the one in Figure 7.24. In the signal captured 
and then processed in the frequency field, characteristic frequencies and 
eventually other characteristic parameters should be monitored; these 
parameters should take the shape of a state vector. A neural network 
trained with this sort of dataset provides at the output the interpreta- 
tion of those data from a diagnostic viewpoint. Using a simple algorithm, 
these output parameters are related to different types of defects for the 
monitored/diagnosed element; the probability of the occurrence of these 
defects is indicated. 


Defect 

types 

Identifying 

algorithm 

Output data 

Neural 

network 


Input data 
(characteristic 
parameters) 

Captured 

signal 



Figure 7.24 Principle schema of diagnostic and classification system for 
defects on basis of power spectrum recognition. 


Elaboration of these neural networks represents a major step in 
developing diagnostic expert systems. Since the algorithmic method of 
learning is generally accepted for expert systems over the method of 
learning from experience, two components are necessary to build neural 
networks: a memorized list of rules and a set of procedures that allows 
conclusion interpolation by utilizing those rules and experimental data. 
Under these circumstances, the information neural processing would be 
the most advantageous way to assimilate the experience. 







Final Conclusions and 
Original Contributions 


8 . 1 GENERAL ASPECTS AND 

CONCLUSIONS— RECOMMENDATIONS 

The performance of modern machine tools and the utilization of elec- 
tronics in machining processes headed the growth of the field of technical 
diagnosis, a very important field of research. 

Machine tools are complex dynamic structures that contain vari- 
ous sources of noise and vibration. Consequently, it is supposed that the 
noise and vibration are determined by the characteristics of the com- 
ponent mechanisms and subsystems. Also, both the parameters of the 
working regime and the characteristics of material workpieces have a 
great influence on these phenomena. 

Research for the evaluation of machine tools is now channeled in 
some meaningful directions: 

Diagnosis of wear and wreckage of cutting tools 

Diagnosis of bearings, especially those for main shafts 

Diagnosis of adjusting mechanisms such as the gear box or feed box 




In order to use vibroacoustic evaluation for machine tools, this book 
describes the specific sources of noise and vibration for machine tools, 
the requirements undertaken for the diagnostic systems, and the steps 
needed to build those systems. The book emphasizes the advantages 
of diagnostic systems utilization related to the cost of investment for 
their implementation, and it also draws some conclusions regarding the 
efficiency of the utilization of diagnostic systems for machine tools that 
are integrated into manufacturing technological systems. 

This research starts with a documented study of the fundamen- 
tals of establishing a technical diagnostic for machine tools. It illus- 
trates the diagnostic methods and techniques applied to the elements 
of the machine tools’ kinematics chain: rolling and sliding bearings, 
gear and belt transmissions, and other rolling parts. The researchers 
are interested in the identification and interpretation of these vibra- 
tion sources while, from the acoustic point of view, the modern ten- 
dency is towards the supervision of acoustic emission (/ > 100 kHz), 
a method which gives information from the microphenomena field of 
physics. This book presents current results obtained in the field. It has 
studied the diagnosis of electrical drive engines, the main sources of 
noise and vibration, from the mechanical, magnetic, and aerodynamic 
phenomena points of view. The principal aspects for diagnosis of the 
cutting tools and the machining process are displayed in this research 
because the machine tool can not be separated from the context in which 
it does work. 

In order to determine the relationship among source-signal-work- 
ing condition, one should choose adequate parameters, diagnosis crite- 
ria, proper algorithms, and admissible levels for the chosen parameters. 
These principles have been followed during the elaboration of differ- 
ent diagnostic methods. The authors grouped these in two categories: 
(a) surface diagnosis methods (the peak element method, diagnosis index 
method, kurtosis method, and the shock impulse method); and (b) pro- 
foundness diagnosis methods (the spectral comparison method, evolute 
method, Cepstrum analysis method, and acoustic emission method). 
While the first group shows the working condition and/or the existence 
of a defect, the methods from the second category can be used in order 
to appreciate the type, the place of the defect, and the time until the 
working condition reaches the end. The book illustrates the theoretical 
basis of the above-mentioned methods and refers to their utilization do- 
main and the way they can be used for the diagnosis of machine tools, 
cutting tools, and the machining process. 



The physical model of the feed kinematic chain has been analyzed 
using the diagnosis method in order to see if it is adequate for study- 
ing diagnostic problems of the elements included in its structure. The 
working process on the machine is stable if the structure of the ma- 
chine tool has an adequate dampening capacity in order to minimize 
the forced and self-excitation vibrations; otherwise the working process 
becomes unstable and the machine chatters. The dynamic stiffness of 
the machine tool is very important for determining system stability. 

The results and small differences between the estimates of proper 
structure frequencies using the analytical technique (finite element 
method) and experimental technique (Nyquist diagram) confirm the ad- 
equacy of the model. Moreover, analyzing the deformations of the phys- 
ical model demonstrates that the most important deformation shows up 
in the area of the mechanisms that are inclined to defective wear: ball 
screw-nut transformation mechanisms and rolling slides of the longitu- 
dinal carriage. These failures can influence the working of the mecha- 
nisms, hastening the appearance of the failure processes that can drive 
the mechanism wreckage. However, if the structure’s proper frequen- 
cies are known, one can compare these frequencies with those of all the 
mechanisms from the studied kinematic chain. Based on theoretical and 
experimental research and studying the literature in the field, it can 
be shown that these working frequencies are situated in a smaller field 
compared to the proper frequencies. 

For an element that has mechanical defects, significant increases 
in the power range show up not only for characteristic frequencies but 
also for their higher values; it is possible for resonance phenomena with 
proper frequencies to appear. This can lead to the discovery of mechan- 
ical defects and to stoppage of the machine tool without any reason. 

The development of virtual instrumentation is necessary because 
of the rapid transformations in the field of measurement devices, and 
because of the extraordinary performances of computers; this should 
be taken into consideration when establishing the technical diagnosis. 
Introducing computerized machining techniques in the field of signal 
acquisition and evaluation means an important gain regarding the ra- 
pidity, reliability, complexity, and flexibility of data measurement and 
processing. 

This book presents original virtual devices that have been devel- 
oped by the authors for surface and profoundness diagnosis. The devices 
were built using the structures presented in Chapter 3. The front panel 
with instructions was shown and the detailed block diagram presented 


for every built virtual device. The authors have identified the main noise 
and vibration sources from feed kinematic chains. The theoretical ap- 
proach allows the calculation of characteristic frequencies of the motion 
elements from these kinematic chains. 

Preliminary experimental research, started in 1992, used tradi- 
tional acquisition and machining drafts based on the physical model; 
they operated on a testing stand that was conceived by making small 
adjustments to an existing machine. Experiments have been done using 
a multichannel acquisition system (produced by S.C. Emco S.A. with 
an eight- byte board and software created by the company specialists). 
These tests have shown the correlation between the appearance of failure 
factors, then of mechanical defects, and finally the indications of diag- 
nostic methods. The FFT transformation of acquisition signals allowed 
the comparison of the frequencies induced by defects and characteris- 
tic frequencies of the studied mechanisms. It was demonstrated that 
knowing the characteristic frequencies during the analysis of acquired 
vibroacoustic signals is very necessary. 

The utilization of virtual instrumentation for experimental research 
using virtual devices has settled a comparison with the same diagnostic 
methods but accomplished by traditional instruments. The conclusion 
is that the virtual instrumentation is the best because it can be easily 
utilized, and it is very precise. Also, the results indicate that virtual 
instrumentation is satisfactory for use in vibroacoustic diagnosis. The 
virtual instruments achieved are operational and can be found in S.C. 
Eurotest S.A. laboratories; this company offered the necessary logistic 
support throughout the experiments. The methodology and measure- 
ments plan indicated that experimental research imposed by the diag- 
nostic techniques should be done in normal conditions and accelerated 
tests. The utilized diagnostic methods proved that they were efficiently 
chosen, leading to a correct and precise evaluation of the working con- 
ditions of the controlled elements. 

During the trials, the unmistakable correlation between the vibrat- 
ing signal and the working condition for some constitutive mechanisms 
of feed kinematic chains was emphasized. As to the use of noise as a 
diagnosis basis, the main disadvantages of signal catching and process- 
ing were presented. At the same time, the correlation between noise 
and the vibratory signal global level was pointed out. The simultaneous 
utilization of different diagnostic techniques and the accumulated expe- 
rience led to a firm conclusion: monitoring the working condition of the 



machine tool or of the constitutive elements should be realized using 
time domain techniques and establishing the technical diagnosis should 
be accomplished using frequency domain techniques. 

The book utilizes parallel processing theory in order to capitalize 
the research database. With this type of information processing expe- 
riential knowledge can be stored, which is why the neural network was 
created. Two diagnostic neural networks have been elaborated: one for 
choosing the good bearings from bad ones, and the other, more ad- 
vanced, which is capable of classifying the bearings' working condition 
into four categories, ranging from good to bad. Information processing 
using neural networks is a viable competitor to classic techniques for 
monitoring and diagnosis. A comparison between the success ratio of 
these techniques and those based on neural processing will definitely 
favor the latter. 


8.2 ORIGINAL CONTRIBUTIONS 

From the theoretical point of view this book displays the following orig- 
inal aspects: 

The identification of noise and vibration sources from feed kine- 
matic chains 

The elaboration of a reckoning method for characteristic frequen- 
cies of the ball screw-nut mechanism 

The method and calculation of characteristic frequencies for ball 
tankettes 

The introduction of virtual instrumentation for machine tool diag- 
nosis 

The development of three virtual instruments for surface diag- 
nosis, based on verified diagnostic methods, using time signal 
processing 

The elaboration of two virtual instruments for profoundness diag- 
nosis, based on competitive methods from the signal frequency 
processing domain 

The introduction of information neural processing for machine tool 
diagnosis 

The generation of a multicriteria neural network having data from 
experimental research, which creates a network capable of rec- 
ognizing perfectly working bearings from deficient bearings 



The elaboration of a neuronal network able to perceive the working 
condition of the bearings 

The realization of a diagnostic analysis, using the finite element 
method, for the physical model of the feed kinematic chain 
The systematization of theoretical and experimental diagnostic 
methods in order to establish the technical diagnosis of machine 
tools 

The elaboration of some conclusions and recommendations regard- 
ing the diagnosis of vibroacoustic methods and their application 
for machine tools 

From the point of view of experimental research, the original con- 
tributions refer to: 

The investigation of the correlation between some statistical pa- 
rameters of the vibration signal and the working condition of 
some ball bearings and needle bearings 
The experimental identification of characteristic frequencies of some 
elements from feed kinematic chain structures and their compar- 
ison with the theoretic frequencies 
Research aspects of determining the technical diagnosis using fre- 
quency analysis of the vibration signal 
Diagnostic analysis of the physical model stability using a Nyquist 
diagram and the identification of the weak points in the system 
The testing and utilization of original diagnostic virtual devices, 
both for surface diagnosis and profoundness diagnosis 
The testing of the neural networks with experimental unsystem- 
atized data and the checking of their estimation capacity from 
the point of view of the working condition 

8.2.1 Perspectives for Continuing Research 

Based on the virtual instrumentation obtained, the concept of vibro- 
acoustic diagnostic equipment may begin. The virtual instruments can 
be transformed into implemented programs that, together with an ade- 
quate hardware interface, can be easily used by everyone. 

The theoretical and experimental research on working frequencies 
for the mechanisms of machine tool kinematic chains should continue 
based on the research methodology illustrated in the book. Research for 
correlating many different parameters in order to establish the technical 



diagnosis in machine tools and the machining domain should also be 
continued, as should the development of a database regarding vibro- 
acoustic diagnosis and its use in processing information through neural 
network techniques for learning by practice. 

Conceiving and developing diagnostic expert systems specialized 
for machine tools, including the diagnosis of cutting tools and the ma- 
chining process is another area for future work. 



Appendix 1: Surface Diagnosis 
Through Virtual Instrumentation 


Al.l SURFACE DIAGNOSIS 

A1 . 1 . 1 INDKURT Virtual Instrument Applied 

to Some Elements of the Feed 
Kinematic Chain Structure 
(Tested Physical Model) 




Figure Al.l Signal captured from the ways (guidings) of the saddle — ways type tankettes with rolls. 
Diagnosis by Kurtosis method. 






















































Appendix 2: Profoundness 
Diagnosis Through Virtual 
Instrumentation 


A2.1 PROFOUNDNESS DIAGNOSIS 

A2 .1.1 MECOSPEC Virtual Instrument Applied 

to a Radial Ball Bearing Coming from a 
Sheet Metal Rolling Mill 
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Figure A2.1 Signal acquisition and spectrum comparison. 
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Figure A2.2 Representations of the acquisitioned spectrums. 





Appendix 3: General Notions 
Regarding the Diagnosis of the 
Functioning State of Machine Tools 


A3.1 COMPELLED VIBRATIONS 

Compelled vibrations appear because of kinematics and/or dynamic fac- 
tors. They have a permanent occurrence and major consequences on the 
working of the technological equipment, and on the equipment, techno- 
logical process, and workpiece quality. 

The usual classification of this type of vibration incorporates (a) 
compelled vibrations that depend on the working process, including such 
factors as stock variation, the periodic variation of the chip cross-section 
(milling, stitching), workpiece material hardness variation, and work- 
speed variation; being dependent on the working process characteristics 
these vibrations are very difficult to avoid; and (b) compelled vibrations 
that are not dependent on the working process, which appear because 
of a deficiency mounting of the equipment, technological and assembling 
lack of precision of the parts, and because of some specific features. This 
type of vibration can be observed during idle running of the equipment. 

The following are some influence areas of factors that are indepen- 
dent of the working process. 

The floor vibrations have a complex spectrum of periodic and ran- 
dom components and also shocks with a frequency range of 1 , . . . , 40 Hz 



with amplitudes between 0.5, . . . , 15 mm. The functioning of electrical 
engines (especially three-phased asynchronous engines and continuous 
current engines) induces a wider spectrum, 16, ... , 530 Hz, with ampli- 
tudes up to 0.4 mm. 

A special influence on gearing working is made by pitch errors, pro- 
file errors, the wheels' eccentricities, and the axle deformations on which 
they are mounted. Two critical zones have been observed in experiments: 
one corresponding to the specific frequency of the elastic system of the 
teeth and the other corresponding to the specific frequency of the elastic 
system gear-shaft. 

The driving belt transmissions introduce vibrations corresponding 
to specific throbs. They can be computed using the formula: 

t in 

un= — \ — 
t V m 

where l is the length of the belt, m is the weight of the belt, S is the 
belt tension, and n is the pulsation number (n = 1,2,3,...). 

The vibrations caused by the movement of the bearings depends 
mainly on shaft rotation frequency and on the number of rolling ele- 
ments. The functioning of hydraulic systems, of cam mechanisms, Malta 
cross, and the like also produce compelled vibrations of frequencies and 
amplitudes which can be estimated. 

Compelled vibrations, within or outside the process, appear and 
take place simultaneously, which means that vibrational phenomena 
have a complex nature and their consequences should be evaluated tak- 
ing into account all the stages of the working process. 

A3. 2 MEASUREMENT OF VIBRATIONS, 

SPECIFIC QUANTITIES 

The vibrator signal picked up in the measurement points using transla- 
tors is, usually, an aleatory sum of periodic and nonperiodic vibrations. 
Ways of processing and evaluating this signal are described below. 

A3 .2.1 Periodic Determinist Vibrations 

The pure harmonic movement (Fig. A3.1) is characterized by the math- 
ematical equation: 


x{t) = X v * sin (cut + 6) 



Figure A3.1 Pure harmonic movement. 


where X v is movement amplitude (the peak value), 6 is phase difference, 
and lj is movement throb. 

By derivation one can obtain the movement speed and acceleration, 
respectively: 

dx / 7T \ 

v(t) = — = ujX v cosM + 0) = V v sin [cut + 0 + —) 

cp x 

a(t) = r- = —cd 2 X v smtcdt + 6) = A v sin (cdt + 0 + 7r) 
at z 

The speed and the acceleration of the movement are also harmonic, 
having the same throb as the displacement, with a phase difference of 
7 t/ 2 and 7r, respectively. 

The characterization of periodic determinist vibrations is possible 
not only by throb and amplitude but also by defining some features 
related to the process progress during one period: 

Absolute average value (mathematical), Xa‘- 


X A = 


x(t)\dt 


Effective value (square average), X e f (RMS): 




The configuration element, Ff : 



The peak element, F v : 



Note: in the case of pure harmonic movement between Xy, Xa, 
and X e f for Ff and TV, respectively, the mathematical relationships are: 

X e f = ——=Xa = — 

J 2 y /2 a /2 

F - 71 
f 2V2 

Fy = V 2 = 1.414 (« 3 dB) 

The deformation energy (elastic) of the elastic system, W p is 



Among these parameters the most important is X e f because it is 
proportional to the vibration power; this can be seen from the previous 
mathematical formula. 


= 1.11 (« 1 dB) 


A3 .2.2 Aleatory Vibrations 

In this case movement is irregular and not repeated in time; the vibra- 
tion should be monitored constantly (theoretically, infinitely); however, 
this is impossible. In practical terms it can work at specific intervals of 
time, which are called “achievements” of the periodic process. All these 
“achievements,” picked up in similar conditions, form the aleatory pro- 
cess itself. This method requires a very long time; as a result, one can 
introduce some probabilistic features that are capable of characterization 
of the aleatory vibrations. 



The performances can be interpreted in terms of time or overall 
(Fig. A3. 2). The amplitude in a specific time t\ is an aleatory feature, 
and it is characterized by the statistical values during the considered 
“achievements” : 


%i(h) = [xi(t 1 ),x 2 (t 1 ),x 3 (t 1 ) ...} 

If the same probability is considered for every “achievement” xk(t) the 
following features can be described: 


The average value of the amplitude, on time t : 


m(ti) = lim 


1 

n 


^Xkih) 


k—1 


The autocorrelation function, which can appreciate in what measure the 
aleatory process remains identical with itself: 




h h b 


leveling in time 





Figure A3. 2 Aleatory vibrations. 


overall leveling 




The intercorrelation function: 


ip'(t\,T) = lim -Viifti) • y k (h + r) 

k=l 

which estimates the way two aleatory signals x(t), y(t) are similar. 

A3. 2. 3 Ergodic and Stationary 

Aleatory Vibrations 

Stationary processes are those processes for which the average value and 
autocorrelation function do not depend on the specific time tl: 

rn(tl) = m(t2) = • • • = m 
ip(tl, t) = ip(t2, T ) = ■■■ = 

For every “achievement” of a stationary aleatory process one can calcu- 
late the average value and the autocorrelation function, respectively: 

i r T/2 

m x (k) = lim — / Xk(t)dt 

1 JT/2 

1 r T ' 2 

ip{k,r) = lim—/ x k (t) ■ x k (t + t) 

1 JT/2 

This means that an ergodic aleatory process can be characterized by a 
“single achievement” which is a great advantage in practice. 

Other features for these processes inclue the repartition function 
of the amplitudes; F(x) is defined by the probability that the movement 
amplitude is inferior to a given value x : 

F(x) = P (— oo < x(t) < x) = P(x{t) < x) 

This function is a monotone increasing function (Fig. A3. 3) and it helps 
determine the probability that the movement amplitude can be found 
in a given range (a, b): 

P(a < x(t) < b) = F(b) — F(a) 

In addition, the amplitude probability density is the limit ratio between 
the probability that the momentary amplitude of the aleatory movement 


im 



can be found stuck between a given interval and the magnitude of this 
interval when it tends to zero: 


p(x) = lim 


P(x < x{t) < x + Ax) 

Ax 


In conformity with the previous relation it can be inferred that 


p{x) = lim 


F(x + Ax) — F(x) 

Ax 


dF(x) 

dx 


The most well-known analytical form for the amplitude probability 
density function is the Gaussian (Fig. A3. 4). 


A3. 2. 4 Noise and Acoustic Emission 
A3. 2. 4.1 Sound and Noise 

Sound is the sensation that is perceptible by the human ear as a result of 
rapid fluctuations of air pressure; it represents the mechanical vibration 
of an elastic medium in which the energy can be propagated from the 
source by progressive sound waves. The noise is usually described as a 
sound or a sum of undesirable sounds. It is considered to be a byproduct 
of daily activity. 

The characteristic of sound waves is that substantial particles os- 
cillate with respect to an equilibrium position and the wave propagation 
speed (the sound speed) is significantly higher than the oscillation speed 



Figure A3. 4 Gauss form for the amplitude probability density function. 
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Figure A3. 5 Dependency of wavelength by frequency. 


of the particles. The distance traveled by the front wave during an oscil- 
lation period is called the wavelength: A25 = c/ f = c * T, where c is the 
propagation of the sound speed, / is the frequency, and T is the period 
or the duration of a complete oscillation. 

The dependency of the wavelength on the frequency is shown in 
Figure A3. 5, where the medium of sound propagation is air; for this 
medium the speed propagation is given by: 



where 7 O is the ratio between the specific heat on constant pressure and 
the specific heat on constant volume, p is static pressure, and p is the 
mass of the volume unit of the medium. Depending on the distance be- 
tween source and receiver, the sound waves are considered to propagate 
in the shape of progressive plane waves or progressive spherical waves 
while the source can be point like or linear. 

Decibel. The introduction of a decibel measurement scale (dB) 
for a subjective evaluation of sound power is based on Weber-Fechner 
physiological law. According to this law, the subjective sensation is pro- 
portional to the decimal logarithm of the excitation when the reference 
level of the acoustic intensity is zero. The lowest acoustic pressure that 
the human ear can sense is 20 mPa, which is 5 * 109 times weaker than 
normal atmospheric pressure. 

The sonorous pressure, expressed in mPa, varies in a very large 
range from 20 to 108 mPa. Therefore, it is difficult to do mathematical 
calculus using such a scale. The decibel scale (Fig. A3. 6 ) avoids this dif- 
ficulty because it uses an audibility threshold of the value 20 mPa. This 
value is considered to be 0 dB. Consequently, every time the acoustic 
pressure (in Pa) is multiplied by 10, 20 dB will be added to the decibel 
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Figure A3. 6 Decibel scale. 


level; hence, 200 mPa corresponds to 20 dB, 2000 mPa corresponds to 
40 dB, and so on. So, the scale in decibels compresses the values from 
20 to 20 million mPa in a range of 0 to 120 dB. 

Another useful aspect of the decibel scale is that it offers a better 
approximation of the noise threshold by human perception; the human 
ear reacts to relative changes of threshold. On this scale, 1 dB variation is 
the same relative variation no matter where it is placed on the scale, 1 dB 
being also the lowest variation that humans can sense. An increase of 
6 dB represents twice the threshold of acoustic pressure and an increase 
of 10 dB is necessary in order to obtain a sound twice as strong. 

Physiological Perception of the Noise. The human auditory ap- 
paratus, the ear, allows the perception of sounds produced by different 
sound waves with frequencies between 16 Hz and 20 kHz (the audibility 
domain). The maximum sensitivity of the ear is in the range of 2000 and 
6000 Hz. 

If two sounds have frequencies / 1 and / 2 it is said that they are 
separated by the interval /1//2. If /2 > / 1 those two frequencies have 
the bandwidth A / = fl — /2. In the acoustic industry, bandwidths 
having an octave and a tierce (a third of an octave) which have the 
corresponding intervals of 2 and \/2 = 1.26, respectively, are very im- 
portant. The central frequency fc of an octave is the frequency that has 
the limited frequencies fl and / 2 such as fc = \/2 * fl. When band- 
widths have an octave, the normalized central frequencies are 31, 63, 
125, 250, 500, 1000, 2000, 4000, 8000, 16,000 Hz and upward. 

In order for a sound to be perceptible, it is necessary that its 
sonorous intensity have a specific minimal level that depends on sound 
frequency and on the sensibility of the ear. The lower limit of the acous- 
tic pressure, for a given frequency that can be heard by a human being 
is called the audibility threshold. It is considered to be a sound with the 
frequency of 1000 Hz and sonorous pressure of 2* 10 — 5 N/m2. For lower 
frequency sounds (<1000 Hz) the audibility threshold increases. 

On the other hand, very strong sounds create pressure on the 
eardrum that can induce pain. The threshold where the pain appears 
is 2 * 10 N/m2 for the frequency of 1000 Hz. The power of the sounds 
represents a feature according to which the sounds can be arranged from 
weak to strong. Subjective perception of a sound or noise strength de- 
pends on the acoustic pressure level and on its spectral characteristic. 

Figure A3. 7 illustrates Fletcher-Munson izosonic curves. Standard 
sound has the frequency 1000 Hz for an acoustic pressure of 1 dB. 
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Figure A3. 7 The Fletcher-Munson izosonic curves. 


In order to appreciate the acoustic pressure level the term moderate 
acoustic level has been used. Hence, the measurement apparatus is 
equipped with balanced filters given by the A, B, C, or D curves in 
Figure A3. 8. The most used moderate acoustic level is represented by 
the A curve, especially in industry and transportation. In the aerospace 
industry the moderate acoustic level is illustrated by the D curve. 
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Figure A3. 8 Balanced filters for measurement apparatus. 



A3. 2. 4. 2 Acoustic Emission 


Acoustic emission is the sequence of elastic waves generated by the re- 
lease of the internal energy stored in a structure. It becomes manifest in 
the higher frequency domain (/ > 100 kHz) by elastic waves detected 
as vibrations on the structure surface. Acoustic emission represents a 
nondestructive method used in order to perceive when and where a flaw 
or a crack appears. The nature and the causes of these shortcomings are 
determined using complementary methods. 

There are four main sources of acoustic emission: 

Movements of structural dislocation 

Phase transformations 

Friction mechanisms (microfrictions, microcollisions) 

Formation and development of cracks 

In the case of dislocations (movement of a line imperfection within a 
crystalline structure) which develop like an avalanche, the signal is con- 
tinuous while in the case of phase transformations (the formation of 
martensite in carbon steel) the signal is impulse type and can be de- 
tected for every transformed grain. 

Flaws come into view in the material where the stress outruns 
strength tension. New surfaces appear and there is a release of energy 
that is partially transformed in acoustic emission. The signal is impulse 
type having high frequency. At the same time, friction mechanisms also 
emanate acoustic signals. The signal amplitudes in acoustic emission 
cover a large area; in relative units these amplitudes are 1 to 10 for 
structural movements, 5 to 1000 for phase transformations, and 20 to 
1000 for flaws. For machining with machine tools these sources are: con- 
tinuous and discontinuous chip forming; deformation of the workpiece 
material; cracking of the workpiece or of the tool; or friction between 
workpiece, tool, chip breaker, breakage, and collision of the chip. There 
are also sources that appear from the functioning of the mechanical sub- 
system (bearings, gears) and high-frequency electrical sources. 

The propagation of acoustic emission is similar to radio waves. 
The source emits spherical wave packages that are influenced by the 
surfaces which are intersected; this creates reflections and surface waves 
(Fig. A3. 9). The heterogeneities of the propagation medium distort the 
front waves. Consequently, the mathematical relationships that describe 
the real propagation phenomena become very complicated when one 
needs to locate and measure the sources and effects of the phenomena. 
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Figure A3. 9 Wave packages influenced by surfaces they intersect. 




Further on, the area of utilizing these methods (based on acoustic emis- 
sion) on large steel structures where an uncertainty coefficient can be 
accepted, becomes very narrow. 

In industry acoustic emission is used for controlling wear and 
tool breakage and also for supervising the bearings and hydrodynamic 
bearings. 

A3. 2. 4. 3 Acoustic Emission Measurement 

As shown in the previous section, acoustic emission is constituted from 
surface waves. Special piezoelectric transducers connected to the struc- 
ture that is to be monitored can perceive them very well. The signal 
monitored by the transducer can be processed in two ways: for a contin- 
uous emission the most useful information is given by the voltameter; 
or for impulse type signals one can use the impulse analyzer. Another 
technique is based on counting the impulses using a peak indicator with 
an adjustable threshold level. 

A3. 2. 4. 4 Acoustic Parameters 

Propagation Velocity of Sonic Waves. Sound waves in solids follow 
the equation 



where E is the longitudinal elasticity module, G is the transversal one, 
and pi is the density of the medium, which is homogeneous and isotropic. 

The propagation of the sound waves in liquids follows a similar 
equation, where ko is the compression modulus: 



Sound wave propagation in gases is adiabatic and follows, therefore, the 
formula pV 7 = const. 2, and the propagation velocity is given by the 
relation: 





where 73 = cp/cv is the adiabatic exponent, while R is the constant of 
perfect gases. For the case of sound wave propagation through air, the 
relation remains 


c = 332Wl + 


t 

273 


where t is the air temperature (°C), which leads to determination of the 
sound wave propagation at a temperature of 20° C to be approximately 
340 m/s. 

Sonic pressure, is, by definition, the average of pressure variation, 
meaning: 


Ps = yf (A p) 2 

The acoustic level (level of sonic pressure) is determined as 



where po = 20 mPa is the threshold, considered to be the audibility limit 
for a sound having the frequency of 1 Hz. 

The intensity of the sound can be defined as the average value of 
the acoustic energy that travels through the surface unit, on a direction 
perpendicular to the propagation, in a unit of time. The equation of 
the sound intensity depends on the propagation area: for the progressive 
wave in a free field it is 



and for the progressive wave in a fuzzy field it is 



4 pc 


where p 2 fft is the square average value of the sonic pressure. 

The sound intensity level can be described similarly to the acoustic 
level; that is, 




where Io is the threshold of the acoustic intensity; that is, Io= 10 to 12 
W/m2. It can be pointed out that in normal conditions (po = 1 atm. and 
to = 22°C), the difference between the value of the acoustic pressure and 
acoustic intensity level is, for the same sound, 0.16 dB, and, therefore, 
it can practically be ignored. 

The level of acoustic power of a sonic source can be evaluated by 
using the relation 

L w = 10 log (JPj 

where Po is the threshold of the acoustic power level, Po = 10 to 12 W. 


A3. 3 CORRELATION BETWEEN SOUND 
ANALYSIS AND MACHINE TOOL 
PERFORMANCE: TECHNICAL 
DIAGNOSIS EQUIPMENT 

A3 .3.1 Time Analysis 

The acoustic signal can be considered to be created by a sum of pure 
harmonic pulsations having different intensities and frequencies, and 
therefore 
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As previously presented, it is possible to characterize a signal by using 
a series of parameters that define the progress of the signal in time 
(usually, for one period): absolute arithmetic mean, xA\ effective value 
(of the squared average), xef; and pick value, xv, xv = max [#(£)]. 

Correlation analysis is the most utilized tool for time analysis, be- 
ing recommended for linear systems that operate either with continuous 
or discrete signals, especially when the ratio between signal and noise has 
a small value. Acceptable input values are either stochastic signals or pe- 
riodic signals, leading, as a result, to the correlation functions and, as a 
particular case of these, the weighted function. The correlation functions 
define the level of similarity between two given signals in accordance with 
the time delay between them (intercorrelation function): 

T 

Rab{r) = lim 2- f a(t) ■ b(t + r)dt 
T — >oo ZJ J 



or, even, the level of similarity of a signal with itself, in accordance with 
the time delay (self-correlation function): 


T 

R aa {r) = lim A [ a(t) ■ a(t + r)dt 
T— >oo 21 J_ T 


For the case of a invariant stable and linear system, by noting with 
h(k) the weighted function or the reaction of the considered system at an 
impulse, the self-correlation function and the intercorrelation function 
are connected by a relation having the form: 

OO 

Rab(r) = h(k) ■ Raa{r - k) 
k = 0 

A correlation graph is a graphical representation of the self-correlation 
function in accordance with the time delay. Its form suggests the contents 
of the signal’s frequencies. The correlation graph gives a peak value for 
t = 0, which is sharper as the contents in high frequencies are richer; 
this value is equal to the square mean value of a random process, being, 
therefore: 

T 

Raa( o) = lim A [ a(t)a(t)dt = x 2 (t) 

T —>oo 21 J 


For a linear system, the integral is proportional to the energy along 
the interval considered; by dividing this value by the length of the in- 
terval one can obtain the average power along that interval, and this is 
the physical significance of the value of the square mean. 

Figure A3. 10 presents some significant types of correlation graphs: 
(a) the correlation graph of a random and stationary ideal process (no 
noise); (b) the correlation graph of a nonperiodic process for large values 
of the time displacement t (the self-correlation functions tend to attend 
the square mean value); (c) idem, when the square average value is 
zero; and (d) the correlation graph of a random process that conceals a 
periodic phenomenon (the correlation graph tends to become a periodic 
function in time t). 

By calculating the Fourier transform of the self-correlation function 
one can obtain the spectral density function of the square mean value 
(power spectral density): 

/ OO 

R{T)e~ l2 * fT dT = S(f) 

-oo 



Figure A3. 10 Significant types of correlation graphs: (a) correlation graph 
of random and stationary ideal process (no noise); (b) correlation graph of 
nonperiodic process for large values of time displacement t (self-correlation 
functions tend to attend the square mean value); (c) idem, when the square 
average value is zero; (d) correlation graph of random process that conceals a 
periodic phenomenon (correlation graph tends to become a periodic function 
in time t). 



This function depicts the manner in which the square mean value (and 
the average power, also) is distributed in the domain of frequencies. 
A periodic vibration can be represented in the frequency domain by a 
band formed by discrete lines, each line representing the square mean 
value of the respective harmonic component, while a random vibra- 
tion determines a continuous band in the frequency domain; so, the 
measured value of the square mean for a given frequency depends on 
the utilized bandwidth; that is why the spectral power density is use- 
ful. The self-correlation function and the spectral density power form a 
pair of Fourier transforms in such a manner that the first one can be 
rapidly obtained by applying an inverse Fourier transformation on the 
second one. 

A3. 3. 2 Analysis in the Frequency Domain 

The acoustic signals can not be analyzed by studying only the amplitude- 
time characteristic because this does not furnish sufficient data for a 
diagnostic interpretation. 

The separation of the vibrations in individual frequency compo- 
nents is called frequency analysis (Fig. A3. 11), this being a technique 
that can be considered as the fundament of the diagnosis, which is based 
on studying vibrations and acoustic signals. The curve that indicates the 
amplitude of the vibrations versus frequency is called a spectrogram. 

The most sophisticated and most precise techniques for establish- 
ing the diagnosis of the performance of machine tools and equipment are 
based on signal analysis in the frequency domain. This analysis can be 
accomplished using various protocols, depending on the element studied, 
diagnosis method, and apparatus controlled by the researcher or existing 
in a research laboratory. 

Gauge apparatuses for measuring vibrations and noise indicate a 
unique band of signals evaluated along the entire bandwidth (Fig. A3. 12a) . 
To evaluate the individual frequency components a filter is utilized, 
which will allow passing only those signals' components within a tide 
bandwidth. The filter's passband is successively displaced along the en- 
tire domain and as a result a reading of the level of vibration is obtained 
in the bandwidth. The filter can be formed either by a series of indi- 
vidual filters for fixed frequencies, which are adjacent and successively 
scanned (Fig. A3. 12b), or by a unique adjustable filter displaced within 
the studied frequency domain (Fig. A3. 12c). 


Time 


Frequency 





Frequency 


Figure A3. 11 Frequency analysis. 

A3. 3. 2.1 Frequency Analysis for Tight Bandwidth 

This kind of analysis is the most common vibration processing procedure 
and, due to its precision and the accuracy of the components of the 
processed signal, the most utilized method for monitoring and diagnosis. 

One can utilize two basic types of filters: constant bandwidth type, 
having an absolute bandwidth of 3 Hz, 10 Hz, and so on, and constant 
proportional bandwidth type, having a bandwidth expressed as a per- 
centage (3%, 10%, etc.) of the selected central frequency. 

In Figure A3. 13 the difference between these two types of filters 
is presented. It should be stressed that the constant proportional band- 
width filters were built in order to maintain a constant bandwidth on 
logarithmic scales of frequencies, which are ideal for large bandwidths. 
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Figure A3. 12 (a) Unique band of signals evaluated along entire band- 

width; (b) series of individual filters for fixed frequencies, adjacent and succes- 
sively scanned; (c) unique adjustable filter displaced within studied frequency 
domain. 
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Figure A3. 13 Difference between constant bandwidth and constant propor- 
tional bandwidth types of filters. 


However, if the frequency scale is linear, a constant bandwidth filter will 
give a constant representation, while the constant proportional band- 
width filter will show an enlarged bandwidth, which is not an advantage 
in practical cases. 

Constant proportional bandwidth filter analysis shows the natural 
response of the systems to mechanical vibrations and permits a com- 
pact representation of a large bandwidth; that is why this is the most 
commonly utilized method for measuring the vibrations. 

Constant bandwidth filter analysis is utilized for high frequencies, 
especially for a linear scale, in order to distinguish the harmonic com- 
ponents. 

The filter-pass bandwidth establishes the resolution of the fre- 
quency analysis to be obtained. The employment of a filter having a 
tight bandwidth offers numerous details and permits the isolation of in- 
dividual peaks within the band, but has, at the same time, the disadvan- 
tage of a increased processing time together with a narrowed bandwidth. 
Sometimes a reduction of the processing time is possible through time 
compression, that is, fast-forward playing of the recorded signal. In that 
way the filter’s bandwidth increases proportionally, and this leads to an 
increased scanning rate along the frequency domain. 

The ideal filter should let pass all the frequency components that 
appear in the bandwidth while eliminating all other frequency compo- 
nents. Practically, electronic filters have no perfectly vertical limits, and 
therefore do not totally eliminate the components outside the bandwidth 
domain. 

In practice two methods are utilized for measuring the filter’s band- 
width: as the bandwidth of an ideal filter that lets pass the same quantity 
of power coming from a white noise source such as the described filter 
(Fig. A3. 14a); and as the bandwidth of a filter that shows an altering of 
3 dB in rapport with the normal transmission level (Fig. A3. 14b). 

The 3 dB bandwidth will differ considerably from the bandwidth 
of effective noise only for low selectivity filters, and therefore one can 
consider that the two definitions lead to the same practical result. 

A3. 3. 2. 2 Frequency Analysis Using Fourier Series 

By using the Fourier theorem and for the given Dirichlet conditions, 
the periodic deterministic vibrations can be looked upon as a sum of 
harmonic pulsations, having frequencies equal to multiples of a funda- 
mental frequency. In this case, the transform in the frequency domain 
can be accomplished using Fourier series by decomposing; during the 
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Figure A3. 14 Two methods for measuring filter bandwidth used in practice. 


transforming process only one period of the signal is usable. The trans- 
formation relations (Fig. A3. 15a) convert the continuous periodic signal 
from the time domain into a discrete band in the frequency domain, a 
band that encloses all the harmonics of the signal. The inverse situation 
is also possible by which a discrete signal in time is transformed in a 
periodic frequency (Fig. A3. 15b). One can observe that due to the sym- 
metry and the periodicity of the frequency band a component having the 
frequency fc in the continuous signal will appear in the discrete signal 
at the frequencies fd = n * fs ± fc , where fs is the discrete frequency 
and n = 0,d=l,±2, ... . In this case, in order to avoid ambiguities in 
the frequency contents of the continuous signal, a pass filter is utilized, 
which is not as broad as the half value of the discrete frequency. 

The transform relations shown in Figure A3. 15 prove the basic 
symmetry of the Fourier transform between the time domain and the 
frequency domain. 

A3. 3. 2. 3 Frequency Analysis Using Discrete 

Fourier Transform 

The discrete Fourier transform (DFT) is applied to discrete and periodic 
signals in the time domain and the result is, also, a discrete and periodic 
signal, but in the frequency domain (Fig. A3. 16). Due to the periodicity 
in both domains, only a finite number of samples are employed, therefore 




Figure A3. 15 Transform in frequency domain can be accomplished using 
Fourier series. 

the transform can be calculated using digital processing, for instance, if 
a period of a signal is described in the time domain by N samples along 
one period. 

The determination of a discrete Fourier transform having N values 
within a sequence implies N * N = N2 multiplication and summation 
operations of complex numbers, which determine a rapid increase of 
the calculation time for a larger number of samples N. Nevertheless, 





Figure A3. 16 Discrete Fourier transform (DFT) is applied to discrete and 
periodic signals in the time domain. 


due to the symmetry of the frequency band, only TV/ 2 of them will be 
independent. In addition, by utilizing an anti-over imposing filter (having 
a bandwidth narrower than half of the sampling frequency) the number 
of significant frequency components is reduced even further. Usually, for 
a signal in time having 1024 samples only 400 frequency lines from the 
band will need to be processed. 

A3 .3.2.4 Frequency Analysis Utilizing Fast 

Fourier Transform 

The algorithm of the fast Fourier transform (FFT) is based on several 
features of the complex exponential function in order to shorten the cal- 
culation time of a regular Fourier transform. In this case, the number of 
operations is diminished from TV 2 to TV * log 27V . The signal processing 
procedure for the case of FFT analysis is presented in Figure A3. 17. 
One can observe that it is necessary to input the recorded signal into 
a sampling block having an analogue-to-digital conversion element be- 
cause the output signal obtained at the output of the transducer has a 
continuous variation. 

The FFT algorithm considers a record from the time domain as 
being a block of TV samples equally separated in time, and which can 
be transformed into another block of TV samples also equally separated 
but in the frequency domain (Fig. A3. 18). The lowest frequency to be 
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Figure A3. 17 Signal processing procedure for the case of FFT analysis. 
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Figure A3. 18 Block of N samples equally separated in time transformed 
into another block of N samples also equally separated but in the frequency 
domain. 




analyzed by FFT is determined by the length of the record in the time 
domain, while the maximum value of the frequency that can be separated 
is / max = (N/2) * (1/Ti?), TR being the record duration. 

The necessary calculations for using the FFT algorithm require a 
finite time, in terms of the number of samples (usually N = 1024), but 
also of the speed of the processor implied. If this calculation (TFFT) 
is less than the duration for the recorded signal (TR), the operation is 
called the analysis in real-time. In this case, the analyzer processes, along 
the duration of a record, a FFT analysis of the preceding record, so no 
information from the time domain is lost, a situation that can occur if 
TFFT is larger than TR. 

The basic relation of the FFT algorithm is still the one known from 
the discrete Fourier transform; that is, 

G W = jfT, 5We- j(2x ‘ n/ff) 

71 = 0 

a relation that can be written using matrix form as 

{Gk} = -jy * {Akn} * {g n } 

where {Gk} and {gn} are the column vectors that include N samples 
from the time domain, and {Akn} is an fV-square matrix that contains 
the complex unit vectors e — j2pkn/N. 

For example, the matrix equation is shown in Figure A3. 19 for 
N = 8. Each arrow in the square matrix depicts a complex unit vector, 
with respect to the attached coordinate system. One can observe that a 
direct calculation of this matrix implies N * N complex multiplication 
operations, which are time consuming. By using the FFT algorithm the 
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Figure A3. 19 Matrix equation for N = 8. 



number of multiplication operations reduces to TV* log 27V, with the con- 
dition that TV is a power of 2. Typically, when TV = 1024, the processing 
time is reduced approximately 100 times. 

Errors Introduced by FFT Analysis. The aliasing effect occurs due 
to signal sampling in the time field; it shows up by the appearance of 
some high frequencies within the field of low frequencies after sampling. 
This effect can be removed by purposely installing a low-pass filter before 
sampling, to eliminate frequencies larger than one half of the sampling 
frequency. 

The time window effect is the result of the finite length of record- 
ing in the time field; the FFT algorithm deals with this record as a 
periodic signal, with the period TR (Fig. A3. 20). This approach is good 
for transitory signals with a period less than TR; the effect is, however, 
harmful for signals with a period larger than TR. Since the signal is “cut 
off” through a rectangular window and then is introduced in a loop in 
order to apply the FFT algorithm, distortions and discontinuities occur 
in transitory areas (Fig. A3. 21a). Consequently, the frequency specter 
will have a few components that do not exist in the original signal. The 
solution is to use a “smooth” window, which has both its value and slope 
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Figure A3. 20 Result of finite recording length in time field treated by FFT 
algorithm as a periodic signal with period TR. 
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Figure A3. 21 (a) While applying the FFT algorithm, distortions and 

discontinuities occur in transitory areas; (b) Hanning window. 


equal to zero at both ends. Usually, this window is a Hanning window 
with period cos 2(2 pt/TR), as in Figure A3. 21b. The round peak of a 
Hanning window can lead to amplitude estimation errors up to 1.5 dB 
(16%). For other types of windows, errors can go down to 0.1 dB (1%). 

The picket fence effect is the result of the discrete sampling of 
the specter in the frequency field. It appears as if the specter is seen 
through the slits of a fence. Consequently, some values, such as peak 
values, can not be observed. The possible resultant error depends on 
how the characteristics of the adjacent filters overlap (Fig. A3. 22). This 
effect is much less as the overlap is bigger. For a Hanning window, the 
distortions introduced in this way do not exceed 1.4 dB, while for a 
rectangular window they reach 3.9 dB. The error can be compensated 
where a frequency component “fits” between two spectral lines. The 
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Figure A3. 2 2 Possible resultant error depends on how characteristics of 
adjacent filters overlap. 


picket fence effect also occurs in a TFD analysis and it is typical when 
third-octave filters are used. 

A3. 3. 2.5 Zoom FFT Analysis 

In a FFT analysis, the resolution of the result is determined by the 
Nyquist frequency (equal to one half of the sampling frequency), and by 
the number of lines in the specter up to the Nyquist frequency. When 
a resolution higher than the one offered by the 400 lines of the ba- 
sic specter is desired, a Zoom FFT analysis is used. The interest area, 
which is between / 1 and / 2, is selected by moving the origin of the 
frequency representation to / 1, simultaneously with the passing of the 
signal through a low-pass filter, in order to eliminate all components, 
except for the range between / 1 and / 2 (see Fig. A3. 23). Zoom FFT 
analysis is useful for processing both low-frequency modulated signals, 
and signals with a large number of harmonics, as well as for separating 
some vibratory phenomena that have very close frequencies. 

A3. 3. 2.6 Frequency Response Function 

The Fourier transform allows a theoretic frequency response of a me- 
chanical system for different types of excitations. For a linear system 



Figure A3. 2 3 When a resolution higher than that offered by 400 lines of the 
basic specter is desired, a Zoom FFT analysis is used. 


with one degree of freedom, the relationship between the excitation f(t) 
and the displacement x(t) of a mass m is given by a linear equation with 
constant coefficients: 


d n x A dx A „ d m f 


‘V 

dx 


+ Bof 


Applying the Fourier transform to both members of this relation, it 
yields: 


n m 

x(f) Mif) k = p(f) E B ^fy 

k = 1 j = 1 


where F(t) and X(t) are the Fourier transforms of the excitation and 
the response, respectively. The complex function defined through the 
relation 


m / n 

3=1 / k= 1 

is called the frequency response function. 

The relation X(f) = H(f) * F(f) shows that the response specter 
is easy to obtain when multiplying the excitation specter by the fre- 
quency response specter. The amplitude of the frequency response for 
each frequency band is the product of the excitation amplitude by the 
frequency response amplitude. The response phase is the sum of the ex- 
citation phase with the response phase. By squaring the relation between 
amplitudes, the power specter formula is derived: 

\ x ( f )\ 2 = \ F ( f )\ 2 *\ H { f )\ 2 

The frequency response function is most often represented through sep- 
arating the real member from the imaginary member. For a proportional 
damping, the real member is zero, and the imaginary member reaches 
its maximum (Fig. A3. 24a). Another way to represent the frequency re- 
sponse function is the Nyquist diagram, by representing the real member 
in terms of the imaginary member (Fig. A3. 24b). 

The graphical analysis of these diagrams, as proposed by Kennedy 
and Pancu as early as 1947 [82], proved to be the most accurate method 
to determine the dynamic parameters and the vibration mode shapes of 
a complex structure. The frequency response function is also useful for 
studying the effects of various excitations upon a mechanical system, for 
determining the mechanical impedance (if the excitation is a force and 
the response is a velocity) , as well as for removing the effects introduced 
by signal propagation in order to rebuild the excitation. 

A3. 3. 2. 7 Analysis in the Amplitude Field 

In general, defects caused by harmful low-speed phenomena (e.g., wear) 
have a typical way of evolution: initially they show up as local singular 
defects, which generate impulse excitations (shocks); the frequency of 
these excitations increases in time and determines the appearance of a 
few peaks in the frequency specter. In order to monitor and diagnose the 
signal in such situations, an analysis in the amplitude field is indicated, 



Figure A3. 24 Frequency response function: (a) for a proportional damping, the real member is zero, and the imaginary 
member reaches its maximum; (b) another way to represent the frequency response function is the Nyquist diagram. 



by means of two specific functions: the peak factor F v and the Kurtosis 
function /? 2 , defined by the relations: 

tv 

X ef 

1 

#2 = -4 / (X - X) 4 p(x)dx 

® J —00 

where X v is the amplitude of the peak signal, X e f is the effective value 
(mean square) of the signal with the amplitude X, X is the average 
value of the signal’s amplitude, p(x) is the probability density of the 
amplitude, and a is the amplitude dispersion. 

The peak factor is utilized especially in the amplitude analysis 
of determinist signals; the Kurtosis function can be used for random 
signals, too. For different situations, the Kurtosis function is given in 
tables in the literature. Note that, for machines and equipment in nor- 
mal operating status, this function is 3, which is considered a reference 
value. An increase of this value in time, determined by the appearance 
of impulse-type processes, shows the appearance and further evolution 
of a defect. 

A3 .3.3 Acoustic Emission Analysis 

More methods are available for acoustic emission study, depending on 
the type of the selected signal. Therefore, a general evaluation is first 
necessary. 

A3. 3. 3.1 Method of Counting Impulses 

This method is useful for impulse- type signals. An evaluation of impulses 
that pass over a previously set threshold value is done. For the evalu- 
ation, the measuring chain consists of a device that differentiates the 
amplitude level of impulses, followed by an impulse counter (per time 
unit or per total). 

The simple counting of impulses (Fig. A3. 25c) can be improved 
through an evaluation of the impulse area (Fig. A3. 25a), in order to take 
into consideration the duration of the impulse, possibly by introducing 
a combination of thresholds (Fig. A3. 25b). 

A3. 3. 3. 2 Method of Impulse Amplitude Mediation 

This method is utilized when the acoustic emission shows up through 
continuous signals. Calculation of the effective value (RMS) is significant 



Figure A3. 2 5 Evaluation of impulses: (a) Evaluation of impulse area; 
(b) combination of thresholds; (c) simple counting of impulses. 


because, as shown in Section A3. 2. 2.1, this value is proportional to the 
signal’s power/energy. This method is sometimes called the energetic 
analysis of impulses. 

A3 .3.3.3 Method of Source Localization 

This method provides information about characteristics and changes of 
the source and of the wave trajectory. Controlled elastic waves are cre- 



FlGURE A3. 26 Locating a source using two transducers. 



ated by means of ultrasound sources. The ultrasonic wavetrain is char- 
acterized by a factor of the simultaneous wave, which depends on the 
succession of impulse frequencies, impulse duration, and the number of 
peaks that exceed a certain threshold. In this way, fracture areas in 
composite materials can be localized in advance. 

In Figure A3. 26, the way in which a source is localized is presented, 
in two dimensions, by means of two transducers. The difference between 
the arrival time of the signal at two transducers determines a plane 
hyperbola, provided the propagation velocity of the signal is known. 
The intersection of the hyperbolas obtained from the transducer pairs 
(1,2), (1,3), and (2,3) defines the real location of the source. This method 
is extremely useful for diagnosing, for localizing the source of primary 
defects, and for reducing the time to restore the operational status. 

A3. 3. 4 CONCLUSIONS 

A3 .3.4.1 Requirements for Diagnostic Systems 

Usage of diagnostic systems should be related to the importance of the 
system monitored within a fabrication process, taking into account its 
complexity and performance. Usually, a monitoring process is conducted 
for complex or continuously working equipment. 

For efficiency in monitoring and detection of defects, the analysis 
and diagnostic have to respect a minimum requirement set: 

Must be adequate for the monitored equipment and have the ability 
and enough accuracy to detect defects. 

Notification of defects must be in due time; false alarms must be 
avoided. 

Must allow localization of defects in order to minimize the inter- 
vention time for repairs. 

Must ensure, as much as possible, a correlation of parameters 
that accompany the equipment work (vibration, noise, temper- 
ature, pressure, etc.), in order to gather complete and correct 
information. 

Must be easy to use and must not raise maintenance problems. 

Must be resistant to shipping and handling, to dust, moisture, and 
industrial liquids, to low and high temperatures, and sometimes 
to radiation. 

Must not be supplied from the same energy sources as the moni- 
tored system, but from special, stabilized, protected sources. 


The reliability of the diagnostic system must be clearly higher than 
the reliability of the monitored system. 

The cost of the diagnostic system and its installation must not 
exceed 10% of the cost of the monitored system. 

A3. 3. 4. 2 Implementation Stages 

Finding a technical diagnostic for complex systems is possible only after 
their dynamics and kinematics are well known. This will lead to deter- 
mination of factors to be monitored, types of transducers to be used, as 
well as locations in which those transducers would be installed. 

When processing signals, the normal and maximum admissible lev- 
els of the monitored parameters should be taken into consideration. The 
sensitive points of diagnostic systems are: the signal processor, which ex- 
tracts the necessary data from the raw signal, and the diagnostic proces- 
sor, which uses these data in order to identify the status of the monitored 
system and to localize and isolate the defect. 
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